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Abstract 


Laser  resonance  saturation  represents  a  powerful  new  form  of  interaction 
that  can  lead  at  high  densities  to  extremely  rapid  and  efficient  coupling  of 
laser  energy  into  both  gases  and  plasmas.  We  have  developed  both  a  comprehensive 
experimental  facility  and  a  theoretical  program  for  studying  these  interactions 
in  alkali-vapours.  Preliminary  calculations  have  suggested  that  laser  resonance 
saturation  could  lead  to  extensive  ground  level  burnout  of  certain  kinds  of 
atoms  or  ions  and  that  this  could  lead  to  the  creation  of  a  ground  level 
inversion  in  these  species.  We  believe  that  this  approach  could  lead  to  the 
development  of  short  wavelength  lasers  of  high  quantum  efficiency. 
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RESEARCH  OBJECTIVES 


Laser  resonance  saturation  represents  a  new  technique  for  rapidly  and 
efficiently  coupling  laser  energy  into  a  suitable  gas  or  plasma.  As  such  it 
could  have  considerable  application  potential.  Our  primary  research  objective 
is  to  study  this  new  kind  of  interaction  in  order  to  gain  a  better  understanding 
of  its  capabilities. 

We  feel  that  laser  resonance  saturation  offers  a  new  approach  towards 
the  development  of  relatively  efficient  short  wavelength  lasers.  This  approach 

is  based  on  rapidly  burning  out  the  ground  level  population  of  a  suitable  i 

ionic  species  within  a  plasma  to  the  point  where  it  becomes  possible  to  create 
a  ground  level  population  inversion  and  thereby  a  short  wavelength  laser  with 
a  high  quantum  efficiency.  Much  of  our  future  research  will  be  directed  at 
investigating  various  aspects  of  this  application  of  laser  interactions  based 
on  resonance  saturation ,  LIBORS. 

STATUS  OF  RESEARCH 

The  consequences  of  laser  saturation  of  a  resonance  transition  within 

an  atom  or  ion  depend  critically  upon  the  length  of  time  saturation  is  maintained,  j 

If  the  laser  locks  the  ground  and  resonance  level  populations  together  for  a 

period  that  is  short  compared  to  the  lifetime  of  the  resonance  level,  then  the 

primary  effect  is  a  burst  of  intensified  spontaneous  emission.  The  author  was 

the  first  to  recognize  the  diagnostic  potential  of  this  so  called  momentary 

(1  2) 

laser  saturation  *  and  subsequently  also  pointed  out  the  highly  perturbative 

(3) 

aspect  of  extended  laser  saturation.  The  word  extended  in  this  context 

relates  to  the  fact  that  the  laser  locks  the  resonance  and  ground  level  popula¬ 
tions  in  the  ratio  of  their  respective  degeneracies  for  a  time  that  is  very  long 

compared  to  the  lifetime  of  the  resonance  level.  Under  these  circumstances 

(A) 

considerable  ionization  and  heating  are  possible. 

During  the  past  year  we  have  continued  our  study  of  laser  resonance 
saturation,  including  both  the  diagnostic  and  ionizational  aspects.  In  this 
section  we  shall  present  only  the  highlights  of  our  recent  work  that  has  been 
(or  is  about  to  be)  published.  Copies  of  these  papers  are  appended  to  this  J 

report.  A  more  complete  description  of  ongoing  work  is  also  provided.  j 
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ABLATION  PLASMA  STUDIES 


As  we  have  indicated  laser  selective  excitation  spectroscopy  involves 

the  momentary  locking  of  the  population  of  two  atomic  levels  in  the  ratio 

of  their  respective  degeneracies  by  an  intense  pulse  of  suitably  tuned 

laser  radiation.  The  subsequent  burst  of  intensified  spontaneous  emission , 

ISE,  from  the  overpopulated  level  represents  a  spectroscopic  signal  that 

is  both  spatially  and  temporally  resolved.  This  powerful  diagnostic 

technique  has  been  used  in  a  variety  of  applications,  including:  atomic 

(5“7)  (8  9) 

lifetime  measurements,  trace  element  analysis,  *  single  atom  detec¬ 
tion,^^  combustion  and  flame  studies,  neutral  hydrogen  measure¬ 

ments  in  Tokamaki^’*^  and  fusion  plasma  impurity  diagnostics . (20,21) 
Ablation  of  a  solid  is  central  to  many  of  the  approaches  towards  attaining 
very  short  wavelength  laser  action.  It  is  to  be  expected  that  a  detailed 
knowledge  of  the  density  ratio  of  consecutive  ionization  stages  in  these 
rapidly  expanding  plasmas  will  be  of  value  if  these  methods  are  to  be 
optimized.  We  have  recently  shown  that  laser  selective  excitation  spectro¬ 
scopy  can  provide  this  information  with  both  spatial  and  temporal  resolution. 

In  a  proof  of  principle  experiment  we  have  used  two  short  pulse, 
appropriately  tuned,  dye  lasers  to  produce  bursts  of  intensified  spontaneous 
emission  (ISE)  on  one  resonance  line  from  each  of  the  two  stages  of  ioniza¬ 
tion  of  interest.  The  experimental  arrangement  is  portrayed  in  figures  1 
and  2.  The  results  (see  Appendix  A)  indicate  that  an  ionization  freeze-out 
at  about  0.1%  ionization  was  attained  in  strontium  2  psec  after  the  moment 
of  laser  ablation.  This  observation  was  made  about  1  cm  above  the  ablation 
target, 

LASER  IONIZATION  BASED  ON  RESONANCE  SATURATION 

Our  current  program  includes  both  a  theoretical  and  experimental  study 
of  the  consequences  of  extended  laser  resonance  saturation.  Our  recent 
theoretical  work  has  involved:  (i)  the  formulation  of  a  simple  model  of 
LIBORS,  (ii)  a  comparison  of  the  results  predicted  by  this  simple  model  with 
those  obtained  from  our  comprehensive  LIBORS  computer  code,*  and  (iii)  a 
study  of  the  influence  of  molecular  nitrogen  upon  the  ionization  and  heating 

*0ur  LIBORS  computer  code  treats  the  laser  pumped  species  as  a  20  energy 
level  system  and  takes  account  of  all  of  the  radiative-collisional  processes. 
It  solves  the  subsequent  set  of  rate  equations  in  conjunction  with  the  energy 
equations  of  the  free  electrons,  the  neutral  and  the  ion  species.  The  code 
provides  the  temporal  variation  in:  the  population  of  each  level,  the  free 
electron  density,  the  ion  and  electron  temperatures  and  the  absorbed  laser 
power  density. 
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associated  with  laser  resonance  saturation.  Our  experimental  program  has 
included:  (i)  some  preliminary  observations  of  the  plasma  emission  created 

through  LIBORS,  and  (ii)  the  measurement  of  an  anomalous  absorption  of  the 
laser  pulse  energy.  Currently,  we  are  assembling  a  new  and  much  improved 
facility  for  undertaking  a  more  definitive  series  of  experiments  aimed  at 

strengthening  our  understanding  of  this  new  form  of  laser  interaction 
in  order  that  we  may  explore  some  of  its  potential  applications. 

LIBORS  Simple  Model 

Although  we  have  developed  a  fairly  comprehensive  computer  code  for 
simulating  the  laser  interactions  resulting  from  resonance  saturation,  we 
have  felt  that  it  should  be  possible  to  formulate  a  relatively  simple  model 
to  describe  LIBORS.  During  the  past  year  we  have  successfully  completed 
this  task.  The  advantages  of  this  simple  model  are:  (i)  it  yields  closed 
form  solutions  for  the  temporal  growth  of  the  free  electron  density,  (ii)  it 
is  not  limited  to  one  alkali  metal,  but  can  readily  be  extended  to  study  the 
ionization  behaviour  of  most  elements  having  relatively  low  lying  energy 
states,  and  (iii)  it  provides  considerable  physical  insight  into  the  relative 
importance  of  the  major  interactions. 

In  this  simple  model  the  ionization  is  perceived  to  proceed  in  three 
stages,  once  the  resonance  transition  has  been  saturated:  (i)  seed  ionization 
creates  an  initial  pool  of  electrons  that  are  free  to  gain  energy  through 
superelastic  collisions ,  (ii)  these  hot  electrons  are  then  capable  of  directly 
ionizing  the  resonance  level  atoms  and  populating  the  intermediate  levels 
which  are  subsequently  photoionized ,  (iii)  finally  direct  electron  impact 
ionization  of  this  intermediate  level  population  leads  to  a  very  rapid  and 
near  complete  ionization  burnout  phase  -  see  figure  3. 

A  detailed  description  of  this  simple  LIBORS  model  is  provided  in  the 

(I 

form  of  two  papers  that  have  recently  been  published.  *  Copies  of  these 
papers  are  appended  as  Appendices  B  and  C.  Consequently,  we  shall  not 
describe  the  simple  model  or  its  solutions  in  this  section.  However,  we 
should  point  out  that  the  paper  presented  as  Appendix  B  deals  with  a  simpler 
version  of  the  model  that  does  not  take  continuity  into  account.  Although 
this  leads  to  a  singularity  in  the  solution  for  the  free  electron  density 
at  burnout,  it  nevertheless  provides  a  reasonable  prediction  of  the  time  to 
achieve  ionization  burnout  (>  95%  ionization)  as  seen  by  the  comparison  with 
the  full  LIBORS  code  for  sodium.  This  results  from  the  very  fast  rate  of 
ionization  in  this  last  phase  of  ionization,  see  figure  7  of  Appendix  B. 
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The  model  presented  as  Appendix  C  takes  account  of  continuity  but  yields 

somewhat  more  complex  solutions  for  the  temporal  growth  of  the  free  electron 

density.  These  solutions  are  found  to  be  in  excellent  agreement  with  the 

predictions  of  our  full  LIBORS  computer  code.  This  agreement  can  be  seen 

to  extend  over  a  dynamic  range  of  six  orders  of  magnitude  in  figure  4,  where 

16  -3 

the  initial  sodium  atom  density  was  assumed  to  be  10  cm  and  the  laser 
6  -2 

irradiance  was  10  W  cm  .  We  have  also  applied  this  simple  model  to  most 

of  the  alkali  metals  and  have  predicted  the  time  for  ionization  burnout 

(>  95%  ionization)  over  a  wide  range  of  densities  and  laser  irradiances, 

see  Table  1  of  Appendix  C.  Figure  5  clearly  reveals  that  this  ionization 

time  is  expected  to  be  close  to  100  nsec  (to  within  ±  15%)  for  most  of  the 

16  -3 

alkali-metal  vapours,  if  the  initial  atom  density  is  10  cm  and  the 
laser  irradiance  is  10^  W  cm 

We  should  also  mention  that  in  this  paper  (Appendix  C)  we  have  shown 
that  if  allowance  is  made  for  the  spatial  and  temporal  nature  of  the  laser 
pulse,  the  apparent  (i.e.9  observed)  ionization  time  is  found  to  be  approxi¬ 
mately  twice  that  for  the  step-like,  uniform  laser  pulse  assumed  in  both 
the  LIBORS  code  and  earlier  version  of  the  simple  model.  In  this  comparison 
the  peak  irradiance  of  the  realistic  laser  pulse  was  set  equal  to  that  of  the 
step-like  pulse.  It  should  be  noted  that  figure  6  and  Table  II  of  Appendix  C 
were  incorrect  and  the  Erratum  presented  here  as  Appendix  0  should  be  sub¬ 
stituted. 

Since  rapid  cooling  of  a  plasma  represents  one  approach  to  achieving 
a  population  inversion  we  have  considered  the  impact  of  a  background  density 
of  molecular  nitrogen  on  LIBORS.  Our  initial  study  has  indicated  that  the 
ionization  time  of  those  solutions  for  which  the  intermediate  phase  dominates 
the  ionization  time  (i.e.,  SI  solutions  in  the  nomenclature  of  Appendix  B) 
are  likely  to  be  quite  sensitive  to  the  presence  of  molecular  nitrogen,  under 
conditions  of  low  laser  irradiance.  This  can  physically  be  understood  by 
reference  to  figure  3  and  figure  6.  For  SI  solutions,  the  ionization  burnout 
time  can  be  approximately  expressed  in  the  form 

t*  *  ~  *n(N  /N  )  (l) 

B  I  e  e 


Under  conditions  of  low  laser  irradiance,  i.e., 


F  «  N.K_  /c(1) 
2  2c  max 


(2) 
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the  intermediate  ionization  rate  coefficient. 


1  ’  "Ac  +  kT  l,  K 

2c  m>3 


a(1) 
2m  me 


where  is  the  saturated  resonance  state  population  density  and  can  be 

1  3 

written  as  for  sodium  where  Nq  is  the  original  atom  density, 

is  the  electron  collisional  ionization  rate  coefficient  for  the 
2c 

resonance  state  |2>, 

K2m  is  the  electron  collisional  rate  coefficient  between  the  inter¬ 
mediate  state  |m>  and  | 2>, 

is  single  photon  ionization  cross  section  for  state  m,  and 
is  the  maximum  value  of  this  cross  section, 

max 

F  is  the  laser  photon  flux  density  (F  »  /  I^(v)dv/h\>) , 

**  .  . 

represents  the  critical  free  electron  density  for  which  the  burnout 

rate  of  ionization  equals  the  intermediate  rate  of  ionization  and 

in  a  like  manner, 

* 

represents  the  characteristic  electron  density  for  which  the  inter¬ 
mediate  ionization  rate  equals  the  seed  ionization  rate,  see 
Appendix  B. 

Jcfc 

In  general  for  SI  solutions  N£  /N  =  100  and  given  (2),  equations  (1) 

and  (3)  indicate  that  is  inversely  proportional  to  and  will  consequently 

depend  strongly  upon  the  free  electron  temperature.  The  low  laser  irradiance 

^  6  - 

criterion  for  sodium  indicates  that  this  situation  will  apply  for  I  <  10  W  cm 

For  much  higher  values  of  laser  irradiance  free  electron  creation  during  the 

intermediate  phase  will  predominantly  arise  through  collisional  excitation  to 

the  intermediate  levels  which  are  subsequently  photoionized.  Under  these 

circumstances  t_  will  be  less  sensitive  to  T  . 

B  e 


Influence  of  Molecular  Gas  ( N„)  on  LIBORS 

The  presence  of  an  appreciable  background  density  of  molecular  nitrogen 
will  not  only  aid  in  rapidly  cooling  the  plasma,  it  will  also  impede  the 
expansion  and  have  a  detrimental  effect  upon  the  degree  of  ground  level 
burnout  achieved.  In  the  case  of  resonance  saturation  of  alkali  metal 
atoms,  the  typical  energy  involved  is  low  enough  (few  eV)  that  only  vibra¬ 
tional  and  rotational  modes  of  need  to  be  considered.  Under  these 
conditions  a  set  of  vibrational-rate  equations  need  to  be  solved  in  conjunction 
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with  the  normal  set  of  LIBORS  equations  and  an  amended  electron  energy 
equation.  If  we  assume  that  only  a  small  fraction  of  the  nitrogen  molecules 
are  vibrational ly  excited  we  can  avoid  having  to  solve  the  set  of  vibrational 
rate  equations  and  can  write  the  electron  energy  equation  in  the  form 


1  dT 

y  k  ’T~  =  LIBORS  terms 
L  t  (see  Eq.  5  of 
Appendix  B) 


-  n  y  <ov>  AE 
o  r  ow  wo 

w>o 


<4) 


where  n  is  the  total  nitrogen  density,  <ov>  the  collision  rate  coefficient 

for  excitation  from  the  ground  vibrational  state  to  one  with  w-vibrational 

quanta.  AE  is  the  corresponding  energy  jump, 

wo 

Solution  of  this  equation  will  yield  a  lower  value  of  than  would  be 
expected  in  reality  and  will  consequently  lead  to  a  conservative  estimate 
for  the  lengthening  of  the  ionization  time.  It  should  also  be  noted  that 
<crv>ow  decreases  as  w  increases  and  so  in  our  preliminary  work  we  have  summed 

to  w  =  4  in  equation  (4).  These  rates  were  numerically  evaluated  from  the 

•  .  (23) 

cross-sectional  data  available  in  the  literature.  The  subsequent  vibra¬ 

tional  cooling  term  in  equation  (4)  is  dependent  upon  the  electron  temperature 

B 

and  is  fitted  by  an  expression  of  the  form  n^AT^  exp(CTe)  in  our  LIBORS  code, 
where  A  =  1 . 315xlO-35,  B  =  7.341  and  C  =  -7.32x10  . 

The  results  of  our  analysis  reveal  that  if  the  alkali  percentage 
seeding  drops  below  about  0.3%  substantial  cooling  results  which  in  turn 
lead  to  a  marked  increase  in  the  time  to  achieve  maximum  ionization  -  and 
furthermore,  full  ionization  (>  93%)  can  no  longer  be  attained;  see  figures 
7  and  8.  A  paper  describing  this  work  more  fully  has  just  been  accepted 
for  publication  in  the  Journal  of  Applied  Physics  and  a  copy  is  appended  to 
this  report  as  Appendix  E. 


Experimental  Facility  and  Preliminary  Results 

Our  first  LIBORS  experimental  facility  included:  a  Phase-R  Model 
DL-21G0A  flashlamp  pumped  dye  laser,  a  specially  designed  and  built  crossed 
heat  pipe  (Heat-cross) ,  an  energy  meter,  two  photodiodes,  two  cameras  and 
two  Heath  monochromators  with  accompanying  photomultipliers.  A  schematic  of 
the  experimental  arrangement  is  presented  as  figure  9.  The  Heat-cross  oven 
used  to  contain  the  sodium  vapour  has  four  windows.  These  provided  access 
for  the  laser  beam  and  enabled  us  to  evaluate  the  absorption  of  the  laser 
pulse  after  transmission  through  the  column  of  sodium  vapour,  the  length 


of  which  was  about  40  cm  and  determined  by  the  position  of  the  water  cooling 
coils  and  the  argon  buffer  gas  pressure.  Observations  of  any  laser  induced 
emission  or  laser  scattered  radiation  can  be  made  through  the  other  pair  of 
windows.  Khodamine  6G  in  ethanol  was  used  in  the  dye  laser  to  provide  output 
energy  pulses  that  ranged  between  100  and  400  mJ.  The  spectral  width  of 
the  dye  laser  was  narrowed  to  about  0.06  nm  by  the  inclusion  of  a  one  to 
two  intracavity  beam  expander  and  a  grating. 

A  50  cm  focal  length  lens  was  used  to  slightly  focus  the  laser  beam 
within  the  Heat-cross  oven,  attaining  a  diameter  of  about  0.50  cm  at  the 
centre.  The  resulting  luminous .plasma  column,  as  photographed  by  Camera 
Cl,  is  dramatically  illustrated  in  figure  10.  Both  t^e  intensity  and  diameter 
of  this  plasma  channel  was  observed  to  decrease  with  progressive  detuning  of 
the  laser,  disappearing  completely  when  the  laser  was  detuned  by  about  1  nn 
to  the  long  wavelength  side  of  the  589  nra  resonance  line. 

Figure  11  displays  both  the  incident  laser  pulse  (from  PD1)  and  a  4  nm 
band  of  emission  centred  about  406  nm  -  which  is  close  to  the  edge  of  the 
radiative  recombination  continuum  to  the  3p  resonance  levels.  It  is  evident 
that  the  continuum  radiation  has  a  much  shorter  rise  time  than  the  laser 
pulse,  but  a  very  much  slower  decay  time.  The  long  decay  time  for  the  emission 
is  consistent  with  recombination  of  the  plasma  channel.  The  delayed,  short 
rise  time  for  the  emission  is  in  general  agreement  with  the  LIBORS  concept 
where  extremely  rapid  ionization  burnout  is  predicted  after  some  onset 
period. 

Unfortunately,  substantial  absorption  of  the  laser  pulse  occurs  at  the 
vapour  pressures  of  interest.  This  means  that  our  present  experimental 
arrangement  is  unsatisfactory  because  the  energy  and  even  the  shape  (and 
more  important  the  onset  time)  of  the  laser  pulse  is  in  fact  unknown  at  the 
observation  region.  This  makes  it  very  difficult  to  relate  the  emission 
results  attained  to  date  to  our  theory.  Consequently,  we  have  designed  and 
built  a  new  kind  of  alkali  oven  that  will  provide  us  with  far  superior  optical 
access.  Indeed,  with  this  oven  we  should  be  able  to  observe  the  region  where 
the  laser  pulse  is  incident  upon  the  alkali  vapour. 

Anomalous  Absorption  of  Laser  Energy 

In  another  series  of  preliminary  experiments  the  laser  was  tuned  to  the 
589  nm  resonance  line  of  sodium  and  the  energy  absorbed  in  transmission 
through  a  40  cm  column  of  sodium  vapour  was  measured  by  a  pair  of  photodiodes. 


7 


These  were  calibrated  against  an  energy  meter  so  that  the  shape  and  energy 
of  both  the  incident  and  transmitted  laser  pulses  were  monitored.  The 
variation  in  the  absorbed  energy  as  a  function  of  the  incident  laser  energy 
was  thereby  evaluated. 

We  were  able  to  demonstrate  that  the  laser  energy  absorbed  can  be 
understood  in  terms  of  simple  radiation  transport  for  low  values  of  laser 
irradiance.  However,  anomalous  absorption  is  observed  to  arise  for  increasing 
values  of  incident  laser  energy.  We  have  been  able  to  show  that  this  addi¬ 
tional  absorption  of  energy  does  not  have  its  origin  in  either  two  photon 
ionization  or  laser  induced  Penning  ionization  of  the  large  population  of 
resonance  state  atoms. 

We  believe  that  this  anomalous  absorption  results  from  the  development 
of  an  appreciable  density  of  free  electrons  (primarily  through  associative 
ionization)  which  are  then  able  to  extract  a  significant  amount  of  energy 
from  the  laser  beam  through  superelastic  collision  quenching  of  the  laser 
maintained  resonance  state  population.  This  interaction,  in  effect,  heats 
the  free  electrons  and  enables  them  to  produce  further  ionization  and 
excitation.  We  view  these  observations  as  representing  our  first  (albeit 
indirect)  evidence  of  superelastic  electron  heating  induced  through  laser 
resonance  saturation. 

(24) 

This  work  was  recently  published  and  a  copy  of  the  paper  is  appended 

as  Appendix  F.  However,  in  this  analysis  we  had  assumed  that  the  resonance 
state  population  can  be  regarded  as  being  in  equilibrium  with  the  radiation 
field  of  the  laser  beam,  i.e.. 


N2(z,t) 


GNo<f>(z,t) 

1  +  <J>(z,t) 


(5) 


where  G  =  g/(l+g),  g  representing  the  resonance  to  ground  level  degeneracy 
ratio, 

represents  the  original  (prior  to  laser  excitation)  ground  level 
density  of  the  vapour, 

4>(z, t) {hF(z, t)/Fg}  represents  the  normalized  laser  photon  flux  density 

where  F(z,t)  represents  the  laser  photon  flux  density,  and 

(24) 

Fg  represents  the  saturated  photon  flux  density. 

This  led  to  a  simple  solution  of  the  radiative  transfer  equation  of  the 
form 

<J>(z,t)  +  £n{<j>(z,t)}  «  <!>(o,t)  +  £n{£(o,t)}  -  CN  o*z  (6) 

o 


I 

where  o  is  the  appropriate  absorption  cross  section  (Appendix  F) . 

Recently,  we  have  checked  the  validity  of  equation  (5)  by  solving 
the  pair  of  coupled  radiative-population  rate  equations  [(9)  and  (10)  of 
Appendix  F]  directly  and  comparing  the  results  with  that  obtained  using  the 
approximation  represented  by  equation  (5).  It  was  found  that  the  complete 
solution  yielded  a  form  of  pulse  sharpening  for  just  the  leading  edge  of 
the  laser  pulse  but  otherwise  made  little  difference  to  the  conclusions 
reached  in  Appendix  F. 

Development  of  New  LIBORS  Facility 

In  order  to  avoid  the  difficulties  encountered  with  our  earlier  experi¬ 
mental  setup,  we  have  constructed  a  new  facility  that  includes  a  more 
powerful  flashlamp  pumped  dye  laser,  a  new  design  of  oven  for  the  sodium 
vapour,  and  a  more  comprehensive  diagnostic  arrangement.  A  photograph  of 
this  facility  is  presented  as  figure  12  and  a  schematic  diagram  which  is 
helpful  for  understanding  the  photograph  is  presented  as  figure  13. 

An  important  feature  of  the  new  oven,  which  is  in  the  form  of  a  heat 
sandwich  rather  than  the  crossed  heat  pipe  used  previously,  is  its  360° 
optical  access.  This  will  enable  observations  to  be  made  of  the  entire 
alkali  vapour  region,  including  the  transition  zones  between  the  cool 
buffer  gas  and  the  alkali  vapour.  This  avoids  the  uncertainties  associated 
with  having  the  laser  pulse  traverse  an  appreciable  column  of  alkali  vapour 
before  reaching  the  observation  zone  -  a  serious  drawback  of  the  crossed 
heat  pipe  design. 

The  new  laser  has  been  found  to  have  at  least  five  times  the  energy  of 
our  earlier  system  and  should  enable  us  to  undertake  experiments  at  relatively 
high  densities  (10*^  -  10*^  cm  .  We  shortly  intend  to  commence  a  definitive 
series  of  experiments  that  will  be  designed  to  test  and  strengthen  our  theore¬ 
tical  understanding  of  the  complex  sequence  of  interactions  that  arise  when 
a  resonance  transition  is  saturated  for  an  extended  period  of  time.  These 
will  include: 

(i)  determining  the  ionization  burnout  time  and  evaluating  its  dependence 
upon:  the  initial  atom  density  of  the  laser  pumped  species,  the  peak 
laser  irradiance  (assuming  a  reproducible  pulse  shape)  and  the  wave¬ 
length  of  the  laser  radiation; 

(ii)  studying  the  influence  of  a  molecular  background  gas  (nitrogen)  on 
this  dependence; 
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(iii)  measuring  the  temporal,  spatial  and  spectral  characteristics  of  the 
laser  pulse  before  and  after  transmission  through  the  alkali  vapour. 

The  ionization  burnout  time  will  be  determined  from  the  temporal  behaviour 
of  the  plasma  emission  and  Stark  measurements  on  selected  emission  lines.  In 
particular  a  number  of  spectral  lines  and  a  segment  of  the  recombination 
continuum  will  be  monitored  with  a  dual  photodetection  system,  comprising 
two  monochromator-photomultiplier  assemblies,  and  a  common  collection  optics 
package,  see  figure  13.  The  temporal  structure  of  the  laser  pulse  before 
and  after  transmission  through  the  alkali-vapour  column  will  be  monitored, 
as  in  our  preliminary  experiment,  by  a  calibrated  pair  of  photodiodes.  In 
addition  we  plan  to  use  a  reticon  (photodiode  array)  to  measure  either  the 
spatial  structure  of  the  laser  beam  or  its  spectral  profile  (using  an  etalon 
arrangement,  see  figure  13)  before  and  after  transmission  through  the  vapour. 

Stark  Measurements 

The  peak  electron  density  and  the  time  of  its  occurrence  will  be  evaluated 

from  measurements  of  the  Stark  profile  of  selected  spectral  lines.  An  example 
.  (25) 

of  the  Stark  profiles,  for  one  possible  line  in  sodium,  corresponding  to 
several  electron  densities  is  presented  in  figure  14.  A  measurement  of  the 
ratio  of  the  intensity  within  the  two  spectral  intervals  shown  in  figure  14 
should  reveal  the  temporal  variation  of  the  electron  density  through  the 
burnout  phase.  The  sensitivity  of  such  measurements  to  uncertainty  in  the 
electron  temperature  can  be  gauged  by  reference  to  figure  15  where  this 
intensity  ratio  is  plotted  as  a  function  of  the  electron  density  for  three 
values  of  the  electron  temperature  that  encompass  the  likely  range  of  tempera¬ 
tures  expected  during  the  burnout  phase.  Indeed,  an  uncertainty  of  100%  in 

the  electron  temperature  is  expected  to  lead  to  only  a  20%  error  in  the 

16  -3 

electron  density  measurement  at  about  10  cm  .  Clearly,  operation  at 
higher  densities  will  provide  an  improvement  in  the  accuracy  of  such  measure¬ 
ments  . 

The  experimental  arrangement  for  this  undertaking  will  involve  imaging 
the  exit  slit  of  a  SPEX  monochromator  onto  the  face  of  a  three  channel 
fiber  optic  array,  see  figure  13.  An  appropriate  choice  of  tuning,  exit 
slit  opening  and  optics  will  then  ensure  that  the  two  outer  fiber  optic  bundles 
sample  the  spectral  profile  in  accordance  with  figure  14.  Use  of  all  three 
channels  should  lead  to  better  accuracy  and  might  enable  us  to  extend  the 
measurements  to  lower  values  of  electron  density. 


Ground  Level  Burnout  Calculation 


Of  considerable  importance  to  the  application  of  laser  resonance 
saturation  to  the  development  of  short  wavelength  lasers  is  the  degree 
of  ground  level  burnout  that  can  be  attained.  Clearly,  the  greater  the 
reduction  in  the  ground  level  population  the  easier  it  will  be  to  create  a 
ground  level  population  inversion.  In  order  to  gain  some  insight  into  what 
might  be  attainable  through  laser  resonance  saturation  we  have  undertaken 

an  extended  computer  run  of  our  L1B0RS  code.  The  result  is  presented  as 

#  .  17  -3 

figure  16  where  it  can  be  seen  that  sodium  vapour  at  10  cm  will  burn  to 

over  99.7%  ionization  within  100  nsec  when  subjected  to  laser  irradiance 
8  -2 

of  10  Wcm  .  Indeed,  after  150.7  nsec  the  percentage  of  ionization  was 


found  to  reach  99.78%  while  the  electron  temperature  attained  19,432  K. 

12 

Furthermore,  the  ground  level  population  burnt  out  to  a  value  of  9.49x10 

-3  - 

cm  corresponding  to  a  ground  level  burnout  factor ,  a  (=  N^/N^)  of  9.49x10 

which  is  very  close  to  the  value  expected  on  the  basis  of  LTE  (local  thermo¬ 
dynamic  equilibrium)  and  laser  saturation,  i.e.. 


N  xlO 
0 


-21e-3.04/kTe 


(7) 


We  assume  here  that 


°  "  N  gN 
e  e 

under  conditions  of  laser  saturation  and  if  LTE  prevails,  then 


N-  =  N  /S(kT  ) 
2  e  e 


where  is  the  free  electron  density  and 

Z+(kT  )  -E  ,/kT 

S(kT  )  s  6xl021  (kT  )3  2  - —  e  c2  e 

e  c  K„ 


if  kT  and  the  resonance  level  ionization  energy,  E  are  in  (eV) .  In  the 

•  4*  ^ 

case  of  sodium  the  ion  partition  function,  Z  (kT^)  =  1,  over  the  range  of 
temperatures  of  interest. 

It  is  quite  apparent  that  LTE  may  be  assumed  to  hold  (at  least  for 
16  —3 

N  >  10  cm  )  during  the  period  of  ground  level  burnout  and  so  we  have 
e 

undertaken  a  calculation  aimed  at  determining  the  limiting  value  for  the 
ground  level  burnout  factor  a. 


In  this  calculation  we  assume  that  a  plasma  comprising  almost  entirely 
of  ions  of  charge  q  is  irradiated  with  intense  laser  radiation  tuned  to 
saturate  the  ground  to  resonance  transition  of  the  ion  of  charge  q.  This 
leads  to  a  burnout  of  the  q-ions  and  the  creation  of  ions  of  charge  q+1.  If 
LTE  is  assumed  to  prevail  we  can  write  a  modified  Saha  equation  of  the  form 

XI  \  —  f  .Ar  Wl/T 


N  N 

S(T  )  a  — 


!  Nq+1  21  ,/2  Zq+1(kT  )  ~{Eq  -AE}/kT 

5 - 6.03xl021(kT  )3/2  - e  c2  q  € 

N2q  C  g2q 


where  Nq+^  represents  the  density  of  ions  of  charge  q+1, 

N2q  represents  the  density  of  ions  of  charge  q  in  the  laser  pumped 
resonance  state, 

Zq+*(kTe>  represents  the  partition  function  of  the  ions  of  charge  q+1, 

g^q  represents  the  degeneracy  of  the  resonance  level  of  the  ion  of 
charge  q, 

^*c2  rePresents  t^e  ionization  energy  of  the  resonance  level  of  the 
ion  of  charge  q, 

AE^  represents  the  ionization  depression  for  the  ions  of  charge  q, 
and 

kT^  represents  the  mean  thermal  energy  of  the  free  electrons  in  (eV)  * 
Continuity  of  mass  leads  to  the  relation 

Nq+1  +  Nq  s  N  q  (12) 

o 

where  N^q  represents  the  original  density  of  ions  of  charge  q.  Continuity 
of  charge  requires  that 

N  £  (q+l)Nq+1  +  qNq  (13) 


Furthermore, 


Nq  *  l  N  q  =  N-q  +  N  q  +  l  N  q 

nJl  "  1  2  nil  n 


Under  conditions  of  laser  saturation  of  the  resonance  transition, 


Nq=l2_N<i  ONq 
2  q  W1  S  N1 


I 


t 


(22) 


2  $ 

Equation  (20)  has  two  solutions,  one  of  which  is  negative  since  (b  +4d}'  >  b 
and  has  no  physical  interpretation.  The  other  can  be  expressed  in  the  form 


Ng  «  [{b2  +  4d>*  -  b]/2 


(23) 


It  is  evident  from  (21),  (22)  and  (23)  that  N  can  be  evaluated  once  the 

e 

electron  temperature  is  known. 

In  order  to  determine  the  limiting  value  of  the  electron  temperature 
we  consider  the  energy  balance  of  the  system.  If  we  neglect  gradient 
effects  then  an  energy  balance  will  arise  when  the  volume  rate  of  laser 
energy  deposition  (through  superelastic  heating)  just  equals  the  volume 
rate  of  radiation  loss  from  the  plasma.  The  primary  contributions  to  this 
radiation  energy  loss  will  arise  from  line  radiation  from  the  excited 
states  within  both  ions  and  recombination  radiation  in  the  ion  of  charge  q. 
We  shall  assume  that  bremsstrahlung  loss  is  negligible  and  that  inverse 
bremsstrahlung  heating  by  the  laser  radiation  is  also  negligible.  At  the 
high  densities  of  interest  in  this  work,  radiation  trapping  may  be  quite 
significant  -  this  would  in  fact  reduce  the  energy  loss  and  allow  the  tem¬ 
perature  to  rise  even  higher  than  expected  by  neglecting  this  effect.  Our 
radiation  loss  is  likely  to  yield  a  pessimistic  value  for  the  ground  level 
burnout  achievable.  On  the  other  hand,  oscillator  strength  information  for 
the  q+1  ion  tends  to  be  rather  sparse  so  that  our  radiation  loss  may  not  be 
as  complete  as  we  would  like. 

Finally,  there  is  the  matter  of  thermal  conduction  loss  from  a  plasma 
of  finite  dimensions  -  this  will  be  dominated  by  electron  conduction  in  most 
instances.  If  we  assume  that  the  plasma  is  cylindrical  in  nature  (being 
formed  along  the  laser  beam)  with  radius  a  and  that  the  lateral  gradient 
is  likely  to  be  much  more  severe  than  the  axial  gradient  then  the  energy 
balance  equation  can  be  expressed  in  the  form 


NeN  qKqEq  =  l  N  q  l  AqEq  +  N  N,q+1  T  0  qEq  +  ^  Nq+1  V  Aq+1Eq+1 

e  4  21  21  L  n  u  nn  nm  e>  ^  n  nr\  L  n  L  nm  nm 


,  n  “  nm  nm 
n>l  m<n 


n>l 


n  cn 


n>l 


m<n 


nm  nm 


9  f  9T  x 

;Ke{^}t 


where  the  electron  thermal  conductivity, 


♦  Vl  *12  Ell 


(26) 


(24) 


K  s  0.234  T5/2/*nA  (Wctn_1  T_1), 
e  e 

is  in  eV  and  f.nA  s  3  for  the  range  of  conditions  of  interest.  flq  represents 


14 


the  radiative  recombination  rate  coefficient  for  level  n  of  the  ion  of  charge 

q,  is  the  superelastic  collision  rate  coefficient,  is  the  q-ion 

energy  difference  between  levels  n  and  m  and  represents  the  Einstein 

nm 

spontaneous  transition  probability  for  the  nm  transition  in  the  ion  of 
charge  q.  If  we  assume  LIE  holds  then  we  can  use 


¥  q  =  N  q 


-Eq„/kT 
q  n2  e 
e 


(for  n  >  2) 


4-1  4.1  “Eqt1/kT 

Kq+1  q+le  m  e 

&n _ 

Zq+1(T  ) 
e 


(for  n  >  1) 


in  the  energy  balance  equation.  Under  these  circumstances  equation  (24) 
becomes  an  equation  in  terms  of  and  T^.  However,  we  saw  from  equation 
(23)  that  can  be  written  in  terms  of  T^.  Consequently,  the  elimination 
of  between  equations  (23)  and  (24)  leads  to  a  complex  transcendental 
equation  for  T£.  This  can  be  written  in  the  form 

■SH  (Te)  =  Lr(T£)  +  Lc(Te)  (27) 

where  (T^)  represents  the  volume  net  rate  of  heating  through  superelastic 

collisions ,  i.e.,  SH  sr  E^fl  ”  e  E^^Te],  under  LTE  conditions, 

L  (T  )  represents  the  volume  rate  of  cooling  through  radiation  loss 
K  e 

L  (T  )  represents  the  volume  rate  of  cooling  through  conduction  loss . 
c  e 

If  we  assume  that  the  temperature  of  the  plasma  column  drops  to  zero  at  the 

boundary,  then  we  can  approximate  the  radial  gradient  in  (27)  by  T^/a. 

SU(T  ),  L  (T  )  and  (L  (T  )  +  L  (T  )}  are  then  evaluated  and  plotted  as  a 
«eRe  Re  c  e  17  -3 

function  of  T  for  the  case  of  sodium  vapour  with  N  »  10  cm  in  figure  17* 
e  o 

Two  temperatures  can  be  determined  from  the  crossover  points  of  these  curves. 

R  c+R 

represents  the  radiative  balance  temperature  and  the  radiative- 

conduction  balance  temperature  for  the  superelastically  heated  plasma. 

C-fR 

T  probably  represents  a  more  likely  temperature  for  determining  the 

®  C+R  C+R  R 

limiting  value  of  the  ground  level  burnout  factor  a  .  Both  and 

are  plotted  as  a  function  of  sodium  density,  N^,  in  figure  18. 

The  general  definition  of  the  GLB  factor  takes  the  form 


a  -  Nlq/vq+1 
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Under  conditions  of  laser  saturation,  we  can  write 


(29) 


which  can  be  calculated  in  terms  of  N  ^  and  T  from  equations  (18)  and  (19). 

®  R  e  C+R 

In  the  case  of  sodium  we  have  used  both  T  and  T  to  determine  the  GLB 

e  e 

factor  and  the  results  are  also  plotted  in  figure  18.  If  only  radiation  losses 
are  considered  it  is  apparent  that  the  lower  the  density  the  lower  the  GLB 
factor  attained.  With  the  inclusion  of  thermal  conduction  losses  it  is  evident 
that  there  exist  an  optimum  density  for  minimizing  a  and  in  the  case  of  sodium 
we  see  that 


C+R  _  -  ,  - 

a  .  =  6.6  x  10 

mm 


17  -3 

N  =  2  x  10  '  cm  J 
o 


As  might  be  expected  from  equation  (24)  and  indicated  in  figure  17  the 
radiation  losses  from  the  species  being  pumped  decreases  with  increasing 
temperature  and  in  the  absence  of  other  losses  the  laser  might  well  be  expected 
to  literally  burn  the  species  out  of  existence.  In  reality,  radiation  loss 
from  the  next  stage  of  ionization  increases  with  increasing  temperature  above 
some  critical  temperature, witness  the  sudden  rise  in  the  radiation  loss  above 
2  eV  in  figure  17.  The  greater  the  ionization  energy  (or  more  properly  the 
resonance  state  energy)  of  the  next  ionization  stage  the  higher  will  be  the 
onset  temperature  for  significant  radiation  loss  from  this  ion  and  in  the 
absence  of  other  losses  (e.g.,  conduction  loss)  the  smaller  the  GLB  factor 
attained.  For  each  of  the  ions  (i.e.,  BeII,  Bill  or  CIV)  suggested  in  our 
proposed  approach  towards  the  development  of  a  short  wavelength  laser  there 
is  a  significant  jump  in  the  ionization  energy  for  the  next  ionization  stage 
and  so  in  all  probability  it  will  again  be  conduction  losses  that  limit  the 
temperature  and  degree  of  GLB  attained.  We  shall  obviously  have  to  study  this 
subject  in  much  greater  depth  and  attempt  to  confirm  experimentally  our 
predictions  on  GLB. 
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NEW  DISCOVERIES  STEMMING  FROM  RESEARCH 


Saturation  of  an  atomic  transition  by  the  intense  radiation  field  of  a 
suitably  tuned  laser  represents  an  important  kind  of  interaction  with  potential 
applications  in  many  areas.  The  consequences  of  laser  resonance  saturation  and 
the  subsequent  applications  stemming  from  this  interaction  depend  to  a  very 
large  extent  upon  the  period  of  saturation.  If  the  resonance  to  ground  level 
populations  are  only  momentarily  locked  (in  the  ratio  of  respective  degeneracies) 
the  principal  effect  is  a  burst  of  intensified  spontaneous  emission,  on  the 
other  hand  an  extended  period  of  saturation  (lasting  for  very  much  longer 
than  the  resonance  lifetime)  can  lead  to  extensive  perturbation  of  the  medium. 
Indeed,  almost  total  ionization  and  very  rapid  heating  can  arise  if  the 
superelastic  collision  time  is  very  short  compared  to  the  duration  of  the 
laser  pulse. 

The  author  was  the  first  to  recognize  the  importance  of  this  interaction 
(see  references  1,  2,  3)  and  many  of  its  possible  areas  of  application.  Sub¬ 
sequent  work  in  the  author’s  laboratory  and  elsewhere  have  proven  that  laser 
resonance  saturation  does  represent  a  significant  form  of  interaction  between 
laser  radiation  and  atomic  vapours  or  plasmas. 

The  combination  of  laser  ablation  and  selective  excitation  spectroscopy 
(EASES)  represents  a  new  approach  at  evaluating  fundamental  atomic  quantities, 
such  as:  radiative  lifetimes,  branching  ratios,  transition  probabilities  and 
selected  collision  cross-sections.  A  preliminary  paper  on  this  subject  was 
published  in  Physical  Review  in  1977;  see  publication  Nos.  7  and  15.  The 
LASES  technique,  as  well  as  being  convenient  and  accurate,  is  particularly 
well  suited  for  measurements  on  short  lived,  highly  ionized  species  created 
by  laser  ablation.  Furthermore,  the  LASES  approach  is  versatile  and  can  use 
multiphoton  or  stepwise  excitation  as  the  means  of  generating  the  bursts  of 
intensified  emission. 

As  a  spin-off  of  this  work  we  have  also  shown  that  the  LASES  concept  can 
also  form  the  basis  of  a  new  form  of  trace  element  laser  microprobe  called  a 
TABLASER;  see  publication  Nos.  10  and  11. 

More  recently,  we  have  demonstrated  in  a  proof  of  principle  experiment 
that  laser  selective  excitation  spectroscopy  can  be  used  to  directly  measure 
with  both  spatial  and  temporal  resolution  the  ion  to  neutral  atom  density  ratio 
in  a  rapidly  expanding  ablation  plasma,  see  publication  No.  2. 

Currently,  our  efforts  are  concentrating  upon  obtaining  a  better  under¬ 
standing  of  the  consequences  of  extended  laser  resonance  saturation  (see 
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Ion— to— neutral-atom  measurements  within  an  ablation  plasma 
through  laser  selective-excitation  spectroscopy 
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A  direct  measurement  of  the  ion-to-neutral-atom  density  ratio  has  been  made  in  a  rapidly  expanding  ablation 
plasma  by  means  of  laser  selective-excitation  spectroscopy.  Evidence  of  an  ionization  freeze  out  at  a  degree  of 
ionization  of  about  0.  Wo  is  presented. 


INTRODUCTION 


THEORETICAL  BASIS  OF  MEASUREMENTS 


Laser  selective -excitation  spectroscopy  (often 
referred  to  as  laser  saturation  spectroscopy  in 
the  literature)  involves  momentarily  locking  the 
population  of  two  atomic  levels  in  the  ratio  of 
their  degeneracies  by  an  intense  pulse  of  laser 
radiation,  as  first  suggested  by  Measures. *“3  The 
subsequent  burst  of  intensified  spontaneous  emis¬ 
sion  from  the  overpopulated  level  represents  a 
spectroscopic  signal  that  is  both  spatially  and 
temporally  resolved.  This  powerful  diagnostic 
technique  has  been  used  in  a  wide  range  of  appli¬ 
cations,  that  include  (i)  atomic  lifetime  measure¬ 
ments  ,3"*  (ii)  trace  element  analysis,7'*  i.iii)  sin¬ 
gle-atom  detection,*"12  (iv)  combustion  and  flame 
studies, l3“ia  (v)  neutral -hydrogen  measurements 
in  Tokamaks,17*48  and  (vi)  fusion-plasma-impurity 
diagnostics,1*20  A  broad  review  of  time-resolved 
laser -fluorescence  spectroscopy  has  been  written 
by  Gauthier  and  Delpeck.24 

Laser  ablation  of  a  solid  is  central  to  many 
diverse  areas  of  endeavour  including  laser  fu¬ 
sion,®  x-ray-laser  development,23'25  trace- ele¬ 
ment  analysis,7'8*27  ion-beam  generation,”  and 
laser  sputtering.29  In  addition,  some  Interest  has 
recently  been  shown  in  using  laser  vaporization 
of  thin  metal  films  for  chemical  kinetic  stud¬ 
ies.29'30  Although  there  have  been  several  investi¬ 
gations  of  the  ablated  material,  some  concerned 
with  the  electrons  and  ions,31, 32  others  with  the 
neutrals,33,34  these  measurements  were  either 
spatially  or  temporally  integrated. 

In  the  present  paper  we  demonstrate  that  laser 
selective-excitation  spectroscopy  can  provide  a 
direct  means  of  evaluating,  with  spatial  and  tem¬ 
poral  resolution,  the  density  ratio  of  two  consecu¬ 
tive  stages  of  ionization  in  a  rapidly  expanding 
ablation  plasma.  This  diagnostic  approach  uses 
two  short  pulse,  appropriately  tuned  lasers  to 
produce  a  burst  of  intensified  spontaneous  emis¬ 
sion  (LSE)  on  one  resonance  line  from  each  of  the 
two  stages  of  ionization  of  the  species  of  interest. 


In  accordance  with  the  basic  Ideas  of  selective - 
excitation  spectroscopy1  the  resonance -to- ground  - 
level  population  ratio,  for  an  ensemble  of  atoms 
excitated  by  a  steplike  pulse  of  laser  radiation 
tuned  to  the  resonance  transition,  approaches  the 
corresponding  degeneracy  ratio  (equivalent  to  an 
effective  Infinite -temperature  distribution)  once 
the  rate  of  stimulated  emission  plus  absorption 
dominates  the  decay  rate  of  the  resonance  level 
(in  the  absence  of  laser  radiation).  This  require¬ 
ment  for  saturation  of  the  transition  can  be  ex¬ 
pressed  in  the  form 

(l  g)f2jl  >>  l/^j  ,  (1) 

where  the  rate  of  stimulated  emission 

Rn  =  Bu  f  l,{v)£u{v)dv/ 4» ,  (2) 


ll{v)  is  the  laser  spectral  irradiance,  £tl(v)  is  the 
resonance  transition  profile  function,  Btl  Is  the 
Milne  stimulated  emission  coefficient  for  the  reso¬ 
nance  transition,  g  is  the  ratio  of  the  resonance  to 
ground  level  degeneracies,  and  rs  is  the  resonance 
level  lifetime  in  the  absence  of  laser  radiation. 

In  the  case  of  broad-band  laser  radiation,  this 
criterion  for  saturation  can  be  expressed  in  the 
form 


where  the  saturated  spectral  irradiance 


(3) 


,  8»hi/»  /r"d\ 


(4) 


Here,  h  is  Planck’s  constant,  v  is  the  laser  (i.e., 
resonance  transition)  frequency,  c  is  the  velocity 
of  light,  and  t™*  is  the  radiative  lifetime  of  the 
resonance  transition. 

For  a  steplike  laser  pulse  the  resonance -to- 
ground-level  population  ratio  Ns /Sx  approaches 
the  degeneracy  ratio  g  with  an  exponential  time 
constant1 
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EXPERIMENTAL  FACILITY  AND  MEASUREMENTS 

Our  facility  for  undertaking  these  measurements 
comprise  two  nitrogen-laser-pumped  dye  lasers 
of  the  Hansch  design37  (excluding  the  etalon),  a 
dual -wavelength  photodetection  system,  and  a 
low-pressure  ablation  chamber  that  has  six  opti¬ 
cal  ports  and  a  rotatable  multisample  target  (see 
Fig.  1)  that  allows  us  to  change  the  target  during 
an  experiment  without  breaking  the  system’s 
vacuum.  A  <?-switched  ruby  laser  is  used  to 
create  the  ablation  Dtasma.  An  RCA  8575  PMT 
and  a  SPEX  1700  monochromator  are  used  in  one 
channel,  while  an  RCA  4526  PMT  and  a  Jobin- 
Yvon  H20  monochromator  are  used  in  the  other 
channel.  The  background  pressure  in  the  chamber 
during  the  experiments  was  less  than  10“*  torr. 

Strontium  was  selected  for  this  work  since  both 
tr.e  neutral  and  singly  ionized  form  possesses 
resonance  line  wavelengths  that  conveniently  fell 
within  the  operating  range  of  our  dye  lasers.  An 
example  of  the  two  intensified  spontaneous  emis¬ 
sion  signals  generated  about  1  cm  out  from  the 
target  surface  and  about  1  psec  from  the  moment 
of  ablation  is  presented  as  Fig.  2.  The  upper 
trace  corresponds  to  the  ISE  signal  frctn  the 
strontium  atoms  at  460.7  run  corresponding  to 
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FIG.  2.  Laser- Intensified  spontaneous-emission 
signals  arising  from  an  ablation  plasma  (1.2  psec  after 
ablation).  The  upper  trace  corresponds  to  the  ISE 
signal  from  strontium  atoms  (at  460.7  nm)  while  the 
lower  trace  corresponds  to  the  ISE  signal  from  strontium 
ions  (at  421.6  nm).  The  first  pulse  on  the  lower  trace 
arises  from  a  photodiode  that  monitors  one  of  the  dye 
lasers.  Time  scales  are  20  n3ec/division. 


FIG.  3.  Laser-intensified  spontaneous -emission 
signals  arising  from  an  ablation  plasma  (1.2  (isec  after 
ablation).  The  upper  trace  corresponds  to  the  ISE 
signal  from  strontium  atoms  (at  460.7  nm)  while  the 
lower  trace  corresponds  to  the  ISE  signal  from  strontium 
ions  (at  421.6  nm).  The  first  pulse  on  the  lower  trace 
arises  from  a  photodiode  that  monitors  one  of  the  dye 
lasers.  Time  scales  are  20  nsec/division.  N.B.  The 
longer  decay  of  the  atom  ISE  signal  relative  to  Fig.  2 
arises  from  radiation  trapping  as  density  of  strontium 
is  higher  in  this  experiment. 

Sr  I  {lPl  -  lS0)  while  the  lower  trace  corresponds 
to  the  ISE  signal  from  the  strontium  ions  at  421.6 
nm  corresponding  to  Sr  n  2Sin).  The  first 

pulse  (not  clearly  seen)  on  the  lower  trace  arises 
from  a  photodiode  that  monitors  one  of  the  dye 
lasers. 

In  our  preliminary  experiments  radiation  trap¬ 
ping  of  the  atom  ISE  signal  was  quite  apparent 
since  the  radiative  lifetime  of  the  lPt  level  of  Sr  I 
was  known38' 39  to  be  about  5  nsec  which  was  close 
to  an  order  of  magnitude  smaller  than  the  ob¬ 
served  decay  time  for  the  460.7-nm  ISE  signal  as 
seen  in  Fig.  3.  In  order  to  eliminate  the  problem 
of  radiation  trapping,  the  strontium  doping  of  the 
target  material  was  reduced  from  10 %  to  below 
1%.  This  had  the  unfortunate  effect  of  also  reduc¬ 
ing  the  already  weak  ion  ISE  signal.  To  compen¬ 
sate  for  this  the  diameter  of  the  421. 5-nm  laser 
beam  was  increased  from  0.25  to  0.8  mm,  an  ef¬ 
fective  increase  in  the  excitation  volume  of  close 
to  a  factor  of  10.  (The  PMT  used  on  the  ion  ISE 
signal  was  also  changed  from  a  degraded  RCA 
31034  to  a  new  RCA  8575.)  A  comparison  of  Fig. 

3  with  Fig,  2  clearly  reveals  the  elimination  of 
radiation  trapping  on  the  atom  460.7-nm  ISE  signal 
in  the  latter  work  (Fig.  2).  The  residual  differ¬ 
ence  between  the  observed  decay  time  and  the 
known  radiative  lifetime  can  be  accounted  for  in 
terms  of  the  finite  duration  of  the  laser  pulse 
(about  5  nsec)  and  the  combined  temporal  response 
time  of  the  RCA  4526  PMT  and  the  Tektronix  8744 
oscilloscope. 

The  temporal  variation  in  the  degree  of  ioniza¬ 
tion  of  the  strontium  in  the  ablation  plasma  (at 
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FIG-  4.  The  percentage  of  singly  ionized  strontium 
(to  neutral  strontium)  as  a  function  of  time  from  the 
moment  of  ablation  and  at  about  X  cm  from  the  target 
as  evaluated  from  the  ratio  of  peak  ISE  signals  from 
SR  n  and  Sri. 

about  1  cm  from  the  target  surface)  as  determined 
from  the  ion  to  neutral  ISE  signal  ratio  ts  pre¬ 
sented  as  Fig.  4.  The  lower  limit  of  the  time 
scale  (about  1  psec)  was  determined  by  the  fLxed 
delay  In  firing  the  nitrogen  laser  with  respect  to 
a  photodiode  signal  from  the  ruby  laser  pulse. 

The  upper  limit  was  given  by  the  weakness  of  the 
ion  ISE  signal.  The  relationship  between  the  ISE 
peak  signal  ratio  S\l/Slp  and  the  ion  to  neutral 
density  ratio  (prior  to  laser  irradiation)  is  given 
in  our  experiments  by  the  expression 

V11  cti 

^■  =  2.3xio->|v.  (9) 

In  the  case  of  our  experiments,  the  saturated 
spectral  irradiance  equation  (4)  is  estimated  to  be 
about  10**  W cm’1  Hz*1  for  each  of  the  transitions 
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(assuming  negligible  collision  quenching),  while 
the  laser  bandwidth  is  close  to  7xlOl°  Hz.  Con- 
sequently,  saturation  requires  that  the  peak  laser 
irradiance  greatly  exceeds  70  W cm**.  Indeed, 
since  the  duration  of  the  laser  pulses  was  only  5 
nsec,  we  would  expect  that  a  peak  laser  irradiance 
of  closer  to  10s  Wcm*1  is  probably  required  to  en¬ 
sure  that  saturation  is  achieved  within  the  period 
of  laser  excitation. 

Unfortunately,  the  dye-laser  power  available  to 
us  was  limited  (by  the  output  of  the  N,  laser)  and 
it  was  found  that  when  the  ion-dye-laser  area  was 
increased  it  led  to  only  marginal  saturation..  The 
consequences  of  this  were  that  the  ion  ISE  signal 
had  a  stronger  dependence  on  the  laser  irradiance 
than  it  would  if  true  saturation  occurred.  Thus 
not  only  did  the  ion  ISE  signal  vary  as  a  result  of 
the  shot -to -shot  Ir  reproducibility  of  the  dye  laser, 
it  was  also  vulnerable  to  some  extent  to  the  fur¬ 
ther  shot-to-shot  variation  associated  with  attenu¬ 
ation  of  the  dye -laser  pulse  in  propagating  through 
the  ablation  plasma.  These  effects  contributed  to 
the  spread  in  the  data  presented  in  Fig.  4.  An 
additional  source  of  shot- to- shot  variation  was 
structure  in  the  expanding  plasma. 

Previous  work  on  this  facility*  has  revealed  that 
laser  scattering  from  ablation  generated  particu¬ 
lates  should  not  constitute  a  source  of  signal  in 
these  experiments  due  to  the  judicious  choice  of 
the  space-time  window  of  observation.  It  is  quite 
clear  that  if  a  more  powerful  Na  laser  were  avail¬ 
able  and  a  few  other  improvements  made  on  our 
facility  (such  as  better  collection  optics  and  better 
photomultipliers)  substantial  improvements  in  ac¬ 
curacy  could  be  achieved.  Nevertheless',  it  is 
apparent  from  our  measurements  that  ionization 
freeze -out*'  occurs  in  our  system  at  about  0.1% 
ionization  in  strontium.  This  represents  the  first 
such  measurement  of  its  kind  and  constitutes  a 
reasonable  proof  of  principle  that  laser  selective- 
excitation  spectroscopy  can  be  applied  to  such 
transient  Ionization  measurements.  If  the  more 
powerful  N,  laser  could  also  be  fired  closer  In 
time  to  the  moment  of  ablation  a  more  compre¬ 
hensive  time  history  should  be  attainable. 
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Laser  ionization  based  on  resonance  saturation — a  simple  model  description 
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We  have  developed  a  simple  physical  model  of  laser  ionization  based  on  resonance  saturation  that  involves  the  most 
important  collisional  and  radiative  interactions,  yet  lends  itself  to  analytical  solutions  that  enable  the  time  history  of 
the  free  electron  density  to  be  evaluated.  We  have  been  able  to  demonstrate  that  in  the  case  of  stniium  the 
predictions  of  this  simple  model  are  within  15%  of  the  values  calculated  by  our  extensive  i  ihors  computer  code  The 
model  has  also  been  used  to  estimate  the  ionization  time  for  each  of  the  alkali  metals  over  a  wide  range  of  conditions. 

These  results  are  found  to  be  consistent  with  several  experimental  observations. 


INTRODUCTION 

Laser-induced  ionization  is  of  conside*  in  - 
terest  both  from  the  standpoint  of  understanding 
the  mechanisms  involved  and  from  attempts  to  ex¬ 
ploit  these  interactions  in  a  diverse  range  of  ap¬ 
plications.  During  the  last  few  vears  there  has 
grown  a  considerable  bodv  of  experimental  evidence 
to  indicate  that  intense  illumination  of  a  resonance 
transition  of  an  atomic  vapor  in  the  density  range 
(lO^-lO17  cm"*)  leads  to  substantial  excitation 
and  ionization.1'”  Indeed,  the  degree  of  ioniza¬ 
tion  achieved  has  been  found  to  approach  100'V 
if  the  duration  of  the  laser  pulse  greatly  exceeds 
the  lifetime  of  the  resonance  state. 

Although  at  low  densities  (below  1C3?  cm*1)  mui- 
tiphoton  processes  tend  to  dominate  the  ioniza¬ 
tion,  for  densities  above  10M  cm'1  collisional  ef¬ 
fects  start  to  plav  an  important  role  and  signifi¬ 
cant  levels  of  excitation  and  ionization  can  be 
achieved  for  quite  low  values  of  laser  irradi - 
ance.'1"”  Lucatorto  and  Me II rath1  2  were  the  first 
to  demonstrate  that  if  laser  radiation  of  relatively 
modest  power  density  (~10,;  Wcm*2)  were  tuned  to 
the  resonance  line  of  an  alkali  metal  vapor  (den¬ 
sity  lO^-lO1”  cm"’),  almost  complete  ionization 
resulted  in  a  remarkably  short  period — far  short¬ 
er  than  could  be  accounted  for  in  terms  of  anv 
multiphoton  or  atom -atom  collisional  process.9 
The  validity  of  these  /'esults  has  been  substantiated 
by  similar  experimental  wc  k  of  Skinner^  and 
Young  ct  al.  1  The  explanation  for  this  efficient 
and  rapid  new  method  of  ionization  can  be  found 
in  a  much  earlier  paper  devoted  to  laser  enhance¬ 
ment  of  ionization  for  MHD  power  generation. 10 
In  that  paper  one  of  the  present  authors  had  pro¬ 
posed  that  under  conditions  of  laser  saturation  of 
a  resonance  transition,  the  energy  of  any  free  el¬ 
ectrons  would  be  rapidly  increased  by  superelas¬ 
tic  collision  quenching  of  the  overpopulated  reson¬ 
ance  level.  Subsequently,  the  rate  of  ionization 
would  be  dramatically  enhanced,  not  only  as  a 
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result  of  the  increased  electron  temperature,  but 
also  by  virtue  of  the  fact  that  the  laser -maintained 
resonance -state  population  constitutes  a  substan¬ 
tial  pool  of  atoms  (roughly  *.)  having  an  ionization 
energy  that  is  reduced  by  the  laser-photon  ener¬ 
gy. 

In  experiments  where  there  is  essentially  no 
initial  ionization,  multiphoton  ionization,9, 11  laser- 
induced  Penning  ionization,'*  or  associative  ioniza¬ 
tion12  could  liberate  the  “seed”  electrons  required 
for  the  above  mechanism.  Although  associative 
ionization  has  been  found  to  have  very  large  cross 
sections  for  high-lying  Rydberg  states, n  an  ef¬ 
fective  value  of  around  2  0‘1T  env  seems  to  he  more 
appropriate12, 1 1  for  alkali  resonance  transitions. 

In  spite  of  this  small  cross  section,  Cardinal1' 
has  shown  that  associative  ionization  would  dom¬ 
inate  the  initial-rate  of  ionization  in  the  experi¬ 
mental  work  with  sodium  of  Lucatorto  and  Me  li¬ 
ra  th. 1,2 

We  have  developed  a  fairly  comprehensive  code 
for  studying  laser  ionization  based  on  resonance 
saturation  (LiBORS).  This  code  treats  the  laser- 
excited  atom  as  a  20-energy -level  system  and 
solves  the  subsequent  set  of  energy  and  population 
rate  equations  by  a  fourth-order  Runge-Kutta 
technique.  A  detailed  description  of  the  code  is 
presented  elsewhere1^"1”  and  will  not  be  repeated 
here.  It  will  suffice  to  mention  that  the  code 
evaluates  the  temporal  behavior  of  (i)  each  of  the 
20-atoniic -level  populations,  (ii)  the  free  electron 
density,  (Hi)  the  free  electron  temperature,  (ivt 
the  ion  and  neutral  temperatures,  (v)  the  laser 
power  absorbed  per  unit  volume,  and  (vi)  the  pow¬ 
er  radiated  per  second  per  unit  volume  as  a  re¬ 
sult  of  radiative  recombination  to  the  resonance 
level,  assuming  a  steplike  laser  pulse. 

The  results  of  our  computer  simulations  for 
sodium  are  consistent  with  the  experimental 
observations  of  Lucatorto  and  Mcllrath1  and  sug¬ 
gest  that  LIBORS  might  be  ideal  for  creating  the 
plasma  guide  channels  that  are  deemed  to  be 
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necessary  for  transportation  of  electron  or  ion 
beams  to  the  fuel  pellet  of  future  particle -beam 
fusion  reactors.17  We  have  also  shown  that  the 
concepts  of  superelastic  laser -energy  conversion 
(SELEC)  could,  under  certain  conditions,  have 
several  advantages  over  inverse  bremsstrahlung 
as  a  means  of  coupling  laser  energy  into  a  plasma 
or  an  un-ionized  gaseous  medium,  fn  the  event 
that  SELEC  is  applied  to  a  strong  transition  with¬ 
in  a  suitable  ionic  species,  extremely  rapid  rates 
of  plasma  heating  can  be  achieved19 — a  feature 
that  is  quite  conducive  to  the  development  of  short  - 
wavelength  lasers.  Indeed,  in  the  case  of  a  boron 
in  plasma,  we  predict19  that  a  heating  rate  in 
excess  of  1013  Ksec'1  could  be  attained  with  a 
laser  irradiance  of  close  to  108  Wcm'2. 

The  purpose  of  the  present  paper  is  to  reveal  a 
relatively  simple  physical  model  of  LIBORS  that 
lends  itself  to  an  analytical  expression  for  the 
temporal  growth  of  ionization,  yet  retains  almost 
all  of  the  important  interactions.  This  model 
allows  us  to  predict  the  ionization  time  under  a 
wide  range  of  conditions  and  is  very  helpful  in 
providing  a  physical  interpretation  of  the  more 
comprehensive  LIBORS  code  results.  Of  parti¬ 
cular  significance  is  the  fact  that  this  simple 
model  enables  us  to  quickly  and  inexpensively 
extend  our  ionization  time  estimates  to  other 
elements.  Several  examples  are  provided  which 
are  found  to  be  consistent  with  experimental 
observations.  In  the  case  of  solium  our  simple - 
model  ionization -time  predictions  have  been  found 
to  agree  to  within  1 5r,  of  our  full  LIBORS  code 
results. 

In  passing,  we  feel  obliged  to  comment  on  the 
term  saturation,  as  used  in  the  above  context. 
Recently,  some  confusion  has  been  displayed  in 
the  literature20'2*  regarding  the  attainment  and 
meaning  of  saturation.  In  particular,  several 
workers  have  suggested  that  their  experiments  have 
failed  to  demonstrate  saturation  of  the  laser  -induced 
fluorescence  (or  intensified  spontaneous  emission) 
even  at  high  valuesof  laser  irradiance.  In  the  context 
of  the  present  pape  r ,  saturation  is  taken  to  mean  that 
the  resonance  to  ground-level  population  densities 
are  locked  in  the  ratio  of  their  degeneracies  bv  the 
laser  fields.  This  concept  was  first  enunciated 
bv  Measures,  71  Recently,  Salter  ct  nl . 71  have 
demonstrated  that  this  condition  can  in  fact  ap¬ 
ply,  even  though  the  laser-induced  fluorescence 
fails  to  saturate  at  the  corresponding  laser  ir¬ 
radiance.  Salter22  has  attributed  this  failure  to 
observe  saturation  to  radiation  trapping.  How¬ 
ever,  Rodrigo  and  Measures25  and  Daily2*  have 
provided  an  alternative  explanation  of  this  effect 
in  terms  of  the  laser  beam’s  radial  profile  and 
the  spatially  averaged  signal.  Some  combination 


of  these  mechanisms  might  in  fact  be  needed  to 
adequately  account  for  this  anomalous  observation. 


LIBORS  SIMPLE-MODEL  FORMULATION 


We  shall  consider  the  case  of  a  homogeneous, 
single -constituent,  atomic  gaseous  medium  sud¬ 
denly  irradiated  with  an  intense  pulse  of  laser 
radiation  that  is  tuned  to  one  of  the  electronic 
resonance  transitions.  Although  the  interaction 
of  this  radiation  field  with  the  atoms  is  most  ac¬ 
curately  treated  through  use  of  the  density-matrix 
approach,  a  rate  equation  analysis  is  reasonably 
reliable27  if  the  bandwidth  of  the  laser  radiation 
is  large  compared  to  the  Rabi  frequency  I  Slr  l.  In 
the  case  of  a  strong  transition,  !ftr  1  =  1 0R v7? 
rad  sec"1,  where  P  Wcm"2  is  the  frequency-aver¬ 
aged  laser  irradiance,  and  a  rate  equation  analy¬ 
sis  will  be  reasonable  for  most  laser-ionization 
experiments  (where  the  typical  laser  line  width 
-0.05  nm)  up  to  a  laser  irradiance  of  about  10a 
W  cm'2. 

Rapid  collisional  dephasing  also  improves  the 
reliability  of  the  rate  equation  analysis.  Under 
conditions  of  resonance  collisional  broadening,  the 
dephasing  time  can  be  less  than,  or  of  the  order 
of,  1  nsec  for  the  densities  of  interest  according 
to  Mcllrath  and  Carlsten.20 

The  model  we  shall  develop  is  particularly  ap¬ 
propriate  to  atoms  with  low-lying  resonance 
states,  such  as  the  alkali  metals.  For  this 
class  of  atoms,  2£81<£el  *  3£21,  where  £21  repre¬ 
sents  the  energy  difference  between  the  ground  and 
the  laser -excited  resonance  state  (and  is  con¬ 
sequently  equal  to  the  laser-photon  energy  hv)  and 
£cl  represents  the  ionization  energy  of  the  ground 
state.  If  the  laser  radiation  is  assumed  to  be 
suddenly  applied  at  /  =  0,  then  the  subsequent  chain 
of  events  can  be  viewed  as  proceeding  in  the  four 
stages  indicated  in  Fig,  1. 

According  to  Measures23  the  resonance  to  ground- 
state  population  densities  are  locked  in  the  ratio 
of  their  degeneracies,  viz. ,  g2/gx*g,  in 

a  characteristic  time  rR  *| (l  +  g)R2l\~l ,  provided 
P(v)>' ls(v)  in  the  case  of  broadband  laser  radia¬ 
tion.  £2i(3^2i  J  l'(v)£{v)dv /4ff]  represents  the 
stimulated  emission-rate  coefficient  for  the  reson¬ 
ance  transition.  fl2l  represents  the  appropriate 
Milne  coefficient  and  £(p)  the  corresponding  line- 
profile  function,  /'(v)  represents  the  laser  spec¬ 
tral  irradiance  and  /s(^)  is  the  saturated  spectral 
irradiance,  given  by 


Birins' 
cz(\  +g) 


where  M'J)  and  r2  is  the  resonance -state 
lifetime  •  'l  2.  is  the  resonance  transition  Einstein 
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probability*  Invariably,  is  very  short  (much 
less  than  1  nsec  for  the  typical  laser  fields  under 
consideration)  compared  to  any  other  process  and 
we  assign  laser  saturation  to  stage  1. 

Once  the  large  pool  of  resonance-state  atoms  is 
created,  several  kinds  of  interactions  give  rise 
to  a  linear  growth  of  the  free-electron  density 
associated  with  stage  2.  The  most  important  of 
these  are 

(i)  two-photon  ionization  of  the  resonance  level: 

A*  +  2hu  —  /I  +  e  +  KE  , 

(ii)  laser -induced  Penning  ionization: 

A *  +A*  +  hv  -  A  4  A  ’  +  c  +  KE  , 

and 

(iii)  associative  ionization: 

A*  +  A*  -  +  <»  4  KE  . 

.-I*  represents  a  laser-excited  resonance-state 
atom,  A  *  the  ionic  species,  and  A  a  ground-state 
atom  of  the  same  species.  A£  represents  the 
dimer  ion.  In  the  case  of  the  alkali  metals,  Gelt- 
man,9  de  Jong  and  van  der  Valk,12  and  Klucharev 
ct  al.u  have  provided  the  relevant  cross  sections. 
The  appropriate  rate  of  ionization  during  stage  2 
can  be  expressed  in  the  form 

1  =  A’c fS'F2  +  t-F  +  jSlo  .r  .  ( 1 ) 

<11  I,**,  2  2c  '  2  L  -  ' 


where  a^2’  (cm'1  sec)  is  the  two-photon,  resonance- 
state,  ionization-rate  coefficient,  vL  (cm1  sec)  is 
the  laser-induced  Penning  ionization-rate  coef¬ 
ficient,  aA  (cm2)  is  the  effective  associative  ion¬ 
ization  cross  section  (allowing  for  the  energy 
defect  of  the  interaction1,1),  r  (cm  sec"1)  is  the 
mean  atom  velocity,  and  F  (photons  cm  “a  sec*1) 
represents  the  laser  photon  flux  density,  i.e.  , 

F—  J  Il(v)dv/hv. 

The  electrons  created  bv  these  initial  processes 
rapidly  gain  energy  through  superelastic  collision 
quenching  of  the  laser  sustained  resonance-state 
population,  as  first  suggested  by  Measures.1” 

In  stage  3,  the  free-electron  temperature29  is 
viewed  as  stabilizing  as  a  result  of  a  balance  be¬ 
tween  the  rate  of  superelastic  heating  and  excita¬ 
tion  cooling.  This  will  be  considered  in  more 
detail  later.  These  electrons  now  give  rise  to  an 
exponential  growth  in  the  free-electron  density 
due  to  both  direct  collisional  ionization  of  the 
resonance  level  and  single  photon-laser  ioniza¬ 
tion  of  the  collisionallv  populated  intermediate 
levels. 

Finally,  in  stage  4,  runaway  collisional  ioniza¬ 
tion  of  the  intermediate  levels  occurs  once  a  cri¬ 
tical  electron  density  is  achieved.  This  process 
leads  to  the  ionization  burnout  that  results  in  al¬ 
most  complete  ionization  of  the  laser  pumped  spe¬ 
cies.  Our  full  LIBORS  computer  code1"'-18  indicates 
that  during  this  rapid  burnout  phase,  superelastic 
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FIG.  1.  Four  stages  of  U BOR S .  Stage  1:  Laser  rapidly  locks  ground  and  resonance-level  imputations  in  ratio  of 
degeneracies.  Stage  2:  <  i  >  Rapid  growth  of  free  electrons  due  to  two-photon  ioni/ation  of  resonance  level  and  laser- 
induced  Penning  Ionization.  Associative  ionization  is  important  for  some  elements,  (li)  Free  electrons  rapidly 
gain  energy  through  superelastic  collisions.  Stage  3:  ti)  Direct  electron- impact  ioni/ation  of  resonance  level  and  sin¬ 
gle-photon  ionization  of  eollisionally  j>opuIated  upper  levels  dominate  the  rate  of  ioni/ation.  iii)  Electron  temperature 
stabilizes  ns  the  rate  of  superelastie  heating  balances  rate  of  eollisional  eooling  through  excitation.  Stage  -4:  (i)  Runa¬ 
way  collisional  ionization  of  upper  levels  occurs  once  a  critical  electron  density  is  achieved,  (ii)  Superelastie  heating 
can  no  longer  balance  eollisional  eooling  and  the  electron  temperature  falls. 
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heating  is  unable  to  maintain  the  free-eiectron 
temperature  and  a  sudden  drop  in  the  temperature 
is  expected. 

In  order  to  obtain  an  analytical  solution  for  the 
ionization  time  from  our  UBORS  simple  model, 
we  shall  assume  that  the  free-electron  tempera¬ 
ture  stabilizes  and  remains  constant  throughout 
both  stages  3  and  4.  Some  justification  for  this 
assumption  can  be  found  in  the  following  argument. 


During  the  linear  growth  period  (stage  2),  the 
dominant  ionization  processes  are  independent 
of  the  free-electron  temperature.  This  is  for¬ 
tunate  because  it  is  questionable  whether  the  elec¬ 
tron  energy  distribution  could  be  characterized 
by  a  temperature  during  this  period.29 

The  appropriate  energy  equation  for  the  free 
electrons  takes  the  form 


+  <2 E2i-Ec2).\X;'F2+  E 

ai 

+  -EjN\oLrF+  ‘.<2£.,  - £  S^tm 

l 

-A*  C  -  A"2//  .  -  A’  A  //  +  A*  y*  K  £  -  ArZf  3{»?)  , 

e  *  e  et  e  <r  e a  e  i—u  cm  cm  e  e  '  * 


where  *e(J)  represents  the  mean  translation  energy 
of  a  free  electron,  Km  (cm^sec'1)  represents  the 
rate  coefficient  for  an  electron  collision  induced 
transition  from  level  »i  to  >/,*''  KmC  (cm1  sec'1) 
represents  the  rate  coefficient  for  an  electron  col¬ 
lision  induced  ionization  transition  from  level 
m,1*  KCm  (cm'  sec"1)  represents  the  three-bodv 
recombination  coefficient, <?^  (cm2)  represents 
the  single-photon  ionization  cross  section™  for 
level  w,  and  the  sum  extends  over  all  levels  w 
> m *  for  which  ionization  by  a  single  laser  photon 
can  be  achieved,  Ar  (cm'1)  represents  the  num¬ 
ber  density  of  atoms  in  the  m  energy  state,  C 
represents  the  loss  of  energy  per  electron  due 
to  the  net  collisionally  induced  upward  movement 
of  bound  electrons  (exclusive  of  the  resonance 
superelastic  term) 


<-vA,  -  XJi -  AVv'3l£21  , 


Emn  (J)  represents  the  energy  separation  between 
levels  w  and  >j,  E^  (J)  represents  the  ionization 
energy  of  level  n,  E  A  (J)  represents  the  energy 
required  for  the  associative  ionization  of  two 
laser-excited  atoms,  Hpi  (Jcm7sec‘l)  represents 
the  rate  of  elastic  energy  transfer  to  ions  through 
Coulomb  scattering  collisions,17  A'a  (cm"1)  repre¬ 
sents  the  total  neutral  atom  density,  H ea 
(Jem7  sec”1)  represents  the  rate  of  elastic  energy 
coupling  between  the  free  electrons  and  the  neu¬ 
tral  atoms.17  {S  (cm7 sec"1)  represents  the  rate  of 
radiative  recombination  to  level  m.  Equation  (2) 
can  be  divided  into  two  equations,  one  for  the  rate 
of  growth  of  the  free  electron  density 


- l*Tthito£'F*  +  »(<7A  +  atf)  +  £  ^a"'r  +  A-,£  JV„Kw-\=£[N^  +  d(m)l),  «> 

Ul  \  *  m*  1  m*  1  f 

where  Te  (K)  represents  the  free-electron  temperature  (viz, ,  ),  and  the  other  for  the  variation  of 

the  free-electron  temperature, 

■v,  ^ k r,>  -  SJiJtnEn  +  (2£„  -  Ee2  -  Ife T'){Nto« ’F  +  |jV|at «•  F) 

+  i(2 En  - -  lkTe)N\c,v  +  Wj£  KcJECm  + § hT,) 

171*1 

A'Jf.  (£„  +  l*r,)+  E,  (*2,  -  -  IkTJS^F  -  NeC  -  -  NX,//,, .  (5) 

1  m>  m* 


Considerable  simplification  of  these  equations 
can  be  achieved  if  we  restrict  our  time  of  inter¬ 
est  to  the  interval  for  which  the  degree  of  ioniza¬ 
tion  is  less  than  a  few  percent.  This  will  incur 
only  a  small  loss  of  accuracy  in  predicting  the 


ionization  time  as  the  ionization  proceeds  from  a 
few  percent  to  almost  total  ionization  in  a  very 
small  fraction  of  the  total  time  for  full  ioniza¬ 
tion,  due  to  the  very  fast  rate  of  ionization  achieved 
in  stage  4— the  runaway  region— a  result  confirmed 
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by  our  LEBORS  code  which  does  not  make  these 
simplifications. 

Under  these  conditions  we  can  write  for  the  al¬ 
kali  metal  vapors,  where  2E2i  ^  Ec{  3E21,  the  rate 
of  ionization  in  the  form 


sion 


Ts* 


E.i2  k 
In  ^1  +  .cS,, 


(9) 


^  -  SJloZ'F  4  i-Vl®,  +  cL  F)  +  .Y,K.J 

4  £  (6) 

The  first  term  on  the  right-hand  side  of  Eq. 

(5)  represents  the  superelastic  heating  rate  and 
it  is  possible  to  show  that  for  times  ^E 
(SjK2XEn)<  where  ±E  ~2En  -  Ec2  -  \k'Te,  this 
superelastic  heating  term  will  dominate  the  heat¬ 
ing  terms  associated  with  the  seed  creation  pro¬ 
cesses.  For  the  alkali  metals,  this  characteris¬ 
tic  time  is  very  short  compared  to  the  ionization 
time  (at  10'  cm*',  it  is  typically  less  than  about 
l  nsec).  The  gain  in  energy  due  to  single-photon 
ionization  and  the  loss  in  energy  from  direct -im¬ 
pact  ionization  of  the  resonance  level  are  also  neg¬ 
lected  in  comparison  to  the  superelastic  heating 
and  excitation  cooling  terms  in  Fq.  (5).  Under 
these  circumstances  the  free-elect ron  energy 
equation,  during  stage  3  can  be  reasonably  well 
expressed  in  the  form 

j-{\ k r ,i * .Y.A,, A',,  _ .v,a,,a\,  - .v,  £  .  <?> 

The  first  term  on  the  right-hand  side  of  Eq.  (7) 
represents  the  superelastic  heating  rate  per  elec¬ 
tron,  while  the  subsequent  terms  represent  the 
cooling  rate  per  electron  associated  with  the  dom¬ 
inant  collisional  excitation  processes.  The  sum 
over  tit  in  the  last  term  extends  over  the  group  of 
intermediate  levels  that  are  rapidly  populated  by 
hot  electrons  and  depopulated  by  single-photon 
laser  ionization. 

Using  electron  collision-rate  coefficients  that 
are  based  upon  Seaton  cross  sections,1  Eq.  (7) 
yields  under  conditions  of  laser  saturation,  a 
transcendental  equation  for  the  slcntly-stntc  s/*/>cr- 
i'lnslic  tt'wfH'nitHrc , 

7-. =  _ _ !^Ji _ 

■’  Infcxpl -(i  k‘J*)  1  x  A*i  ;)>  * 

(8) 

where  fm  ~  E m2  -  /: , .  and  Jim'  \fim$ ./^repre¬ 
sents  the  absorption  oscillator  strength  for  the 
im  transition  and  the  relevant  Gaunt  factor 

averaged  over  a  Maxwellian  velocity  distribu¬ 
tion.  This  equation  can  be  solved  bv  an  iterative 
procedure,  or,  in  light  of  the  logarithmic  nature 
of  the  denominator  and  the  fact  that  Z:^  - /  ,,  (in 
the  model),  it  can  be  approximated  by  the  exp  res - 


where  we  have  also  assumed  L  (njv  1.  This  value 
of  the  free -elect ron  temperature  is  then  used  to 
compute  the  various  collision-rate  coeffic ients. 

The  relevant  intermediate -level  population -rate 
equation  can  be  expressed  in  the  form 

dXm  <it~  X2SJ<2m-Sm(u^F  SJ\J  , 

(10) 

where  .1^- -A  ■>  rfis  represents  the  radiation-trapped 
Einstein  transition  probability  for  the  ms  transi¬ 
tion,  s  can  correspond  to  either  the  ground  or 
resonance  level,  >  is  the  Holstein  escape  fac¬ 
tor/"  and 

A  A'  4  Y\  A'  .  (11) 

m  **tc  j  rm 

nt 

Wo  have  assumed  that  in  light  of  the  large  popula¬ 
tion  of  the  resonance -state  atoms  tdue  to  laser 
saturation).  .V  .V  J\ ,  will  dominate  any  other  rate 
of  population  of  level  m  from  other  levels.  This 
may  not  be  strictly  true  for  closely  spaced  ad¬ 
jacent  levels,  but  such  terms  only  serve  to  rapid¬ 
ly  equilibrate  the  population  within  the  manifold 
of  informed  'to  states,  ff  wo  introduce 

l\  cyp  (12) 

and  set 

p  P  4  .1*  ,  (13) 

then  over  a  significant  fraction  of  the  ionization 
time 

-  .V  A  .  (14) 

r  rr\  .  »*l 

Under  these  circumstances,  Fq.  (10)  has  a  solu¬ 
tion  of  the  form 

SJt)  A\A  >n.oxp(-/\/)  j  S\Sf''expi()J'\ir  (15) 

if  the  vapor  is  assumed  to  be  unexcited  prior  to 
laser  irradiation,  i.o.,  .V^(0)  0.  /•',  A*.,  and 

h2m  are  also  taken  to  be  constant  over  the  period 
of  interest. 

The  constancy  of  A/  over  an  appreciable  fraction 
of  the  time  to  ionize  is  quite  reasonable  on  the 
grounds  that  the  upper  levels  [m  3)  only  become 
significantly  populated  just  prior  to  ionization  burn- 
burnout  (stage  4).  In  fact,  if  we  use  the  saturation 
condition,  we  can  write  A’, ^  (*A„,  where  (»  g(l 
-‘'.g)"1  and  An  represents  the  initial  atom  density. 

The  time  independence  of  A  1|fi  is  based  upon  our 
previous  assumption  of  a  stabilized  temperature 
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during  the  intermediate  phase  (stage  3)  and  F  has 
been  set  to  be  constant  by  our  assertion  of  a  step¬ 
like  laser  pulse.  The  viability  of  these  assump¬ 
tions  will  become  apparent  later  when  we  review 
the  LIBORS  computer-code  results. 

It  is  apparent  that  substitution  of  Eq.  (15)  into 
the  ionization  equation  yields  a  polynomial  expan¬ 
sion  in  Ng  on  the  right-hand  side  of  Eq.  (6).  At 
the  earliest  times  from  the  moment  of  laser  sat¬ 
uration,  i.e.,  stage  2,  the  free-electron-density 
growth  rate  is  expected  to  be  independent  of  Ne , 
while  an  exponential  growth  is  to  be  expected 
during  the  subsequent  period,  stage  3.  This  sug¬ 
gests  that  a  first  approximation  to  the  free- 
electron  density  can  be  expressed  in  the  form 


Xe(t)  *  iVe0[exp(/3/)  -  1J .  (16) 

This  gives  a  linear  growth  for  1  and  an  ex¬ 
ponential  growth  for  /3/  r  1.  Arp(,  can  be  taken  to  be 
the  seed  electron  density  for  which  Xe  departs 
significantly  from  a  linear  growth. 

Substituting  Eq.  (16)  into  Eq.  (15)  yields 


Nm(t)  = 


0  + 


1  -  exp(-p„<)\1 
exp(/J/)  -  1  /J  ' 


(17) 


For  times  sufficiently  great  for  relaxation  of  the 
intermediate-level  populations,  we  obtain  an  ap¬ 
proximation  relation  of  the  form 


S  (/)*  ASA2m-Ve(/) 


provided 

/  >  (1//3)  in(l  +&/pm). 


(18) 

(19) 


In  many  instances  this  is  more  restrictive  a  con¬ 
dition  than  necessary.  Further  discussion  of  this 
will  be  provided  later.  If  Eq.  (18)  is  used  in  Eq. 
(6)  to  eliminate  S m  we  arrive  at  a  simple  form  of 
the  ionization  equation, 


(20) 


where  we  have  introduced 


as  the  critical  electron  density  at  which  the  burn¬ 
out  rate  of  ionization  equals  the  intermediate  rate 
of  ionization.  In  alike  manner  we  can  also  define 


JV*  sS // 


(25) 


as  the  characteristic  electron  density  for  which 
the  intermediate  ionization  rate  equals  the  seed- 
ionization  rate.  It  is  apparent  from  a  comparison 
of  Eq.  (16)  with  Eq.  (20)  that  we  can  identify 
with  /  and  X  n  with  N+ .  In  which  case  we  can  write 

e u  e 


I  —  X 


z 

m*  3 


(26) 


and 


B  =  A’, 


tf+pj 


(27) 


1  can  be  evaluated  in  an  iterative  manner  by  sub¬ 
stituting  pm-  pmix  =  p  into  Eq.  (26)  and  obtaining  a 
first  approximation  of  the  form 


I  ~  'Ak  - p)  +  j |u  -  p)2  +  4 ^<p  +  \2  X) 


(28) 


where  k  ~  X^K.^.  If  this  first  approximation  is  then 
reintroduced  into  the  right-hand  side  of  Eq.  (26), 
we  have  found  that  in  almost  all  cases  /  is  within 
lfvr  of  the  correct  (limiting)  value. 

The  simple -model  LIBORS  ionization  (Eq. 

(20)  |  has  two  forms  of  analytical  solution  de¬ 
pending  upon  the  relative  values  of  S,  /,  and  B. 

If  P<  4 SB }  vhich  is  equivalent  to  X**  <  X^ ,  then 
ionization  proceeds  directly  from  stage  2  to  stage 
4— ^hai  is  to  say,  seed  ionization  leads  directly 
to  runaway  ionization.  The  solution  to  Eq.  (20) 
under  these  circumstances  takes  the  form 

(29) 

where 

a*(4SB -l2)1"2  (30) 


S  +  i.V2r(tr,  +  aL  P)|  (21 ) 

as  the  seed-ionization  rate, 


I  ? 


K^P. 


(a  +  p 


j 


(22) 


as  the  intermediate  ionization-rate  coefficient, 
and 


w+pj 


as  the  burnout-ionization-rate  coefficient. 
We  see  that  we  can  now  define 


(23) 


This  class  of  solutionis  indicative  of  a  very  high 
rate  of  seed  ionization  and  we  shall  refer  to  this 
as  an  S3  solution. 

If  /2>  4 SB,  which  is  equivalent  to  X**  >  A**,  stage 
3  ionization  will  play  some  role  in  separating  stage 
2  from  stage  4— that  is  to  say  intermediate  ioniza¬ 
tion  is  significant  in  increasing  the  free-electron 
density  to  the  point  where  runaway  ionization  oc¬ 
curs.  The  solution  of  the  ionization  Eq.  (20)  under 
these  circumstances  takes  the  form 


Arp  =/  B 


(24)  where 


/ 


k 


/ 


H 10 
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A  :(/2-4S£)I/2,  (32) 


If,  however,  /*»4SB  (i.e.  ,  N**  >iV*),  then  the 
intermediate  processes  play  a  major  role  in  de¬ 
termining  the  ionization  burnout  time  rB.  We  shall 
discriminate  between  solutions  that  arise  when 
100 SB  l2  4 SB  and  those  that  apply  when  V 
>  100SB,  by  designating  the  former  as  S2  solu¬ 
tions  and  the  latter  as  SI  solutions. 

The  ionization  burnout  time  for  l2  >  4SB  solutions 
takes  the  form 


However,  for  SI  solutions,  where  l2  lOOSf*  (or 
N**  100A’*)  we  can  approximate  this  by  the  ex¬ 

pression 


This  relation  lends  itself  to  the  pictorial  inter¬ 
pretation  presented  as  Fig.  2. 

For  intermediate  levels  having  zero  single-phot¬ 
on  ionization  cross  section  at  the  laser  wavelength, 
we  have  used  the  minimum  photoionization  cross 
section  in  Eq.  (18)  since  collision  coupling  be¬ 
tween  adjacent  levels  plays  a  role,  as  stated  earl¬ 
ier,  in  equilibrating  adjacent  levels.  We  have  also 
found  that  in  some  instances,  particularly  S3  solu¬ 
tions  ,  Eq.  (18)  grossly  overestimates  the  inter- 
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FIG.  2.  Tht-  three  distinct  phases  (stages  2  to  4  of  Fig. 
1)  of  laser  ionization  based  on  resonance  saturation  arc 
indicated  in  this  graph  of  electron  density  against  time. 


mediate-level  populations  at  times  corresponding 
to  rfl/3,  as  seen  from  a  comparison  of  Sm  ob¬ 
tained  by  Eqs.  (17)  and  (18).  Under  these  cir¬ 
cumstances  we  have  assumed  that  the  intermediate- 
level  populations  do  not  contribute  to  the  inter¬ 
mediate  ionization  stage  and  reassign  1  = 

In  the  limit  of  J2^  4 SB,  the  ionization  burnout 
time 


(37) 


LIBORS  SIMPLE-MODEL  RESULTS  FOR  ALKALI 
METALS 

In  this  section  we  present  some  representative 
results  of  the  LIBORS  simple-model  calculations. 
The  relevant  cross  sections  and  parameters  re¬ 
quired  for  this  analysis  are  presented  for  each 
of  the  alkali  metals  in  Table  I.  The  type  of  solu¬ 
tion  obtained  and  the  appropriate  ionization  burn¬ 
out  time  for  each  element  are  also  displayed  under 
a  wide  range  of  densities  and  values  of  laser  ir- 
radianee.  A  clear  pattern  is  readily  discernible. 
In  the  case  of  lithium  all  solutions  are  of  the  SI 
type,  which  means  that  the  intermediate  process¬ 
es  are  always  important— at  least  for  10ls  A',, 

'  1017  (cm_<)  and  10  / 1  10K  (Worn"’)  ranges 

under  consideration.  SI  solutions  are  interpreted 
to  mean  that  the  intermediate  processes  plav  a 
dominant  role  and  this  appears  to  be  particularly 
true  where  associative  ionization  is  absent.  The 
predicted  high  value  of  T*  is  also  important. 

SI  solutions  are  also  found  for  sodium  under 
most  of  these  conditions.  Reference  to  Table  I 
shows  that  for  sodium,  T*  is  slightly  higher  and 
<ja  is  slightly  lower  than  for  the  remaining  alkali 
metals  (n>. ,  K,  Rb,  and  Cs). 

At  the  other  end  of  the  periodic  table,  cesium 
has  predominantly  S3  solutions.  This  is  not  sur¬ 
prising  in  light  of  the  predicted  low  value  of  7* 
and  relatively  high  value  of  oA.  It  is  also  appar¬ 
ent  that  potassium  and  rubidium  have  S3  solutions 
at  the  highest  laser  powers.  In  the  case  of  rub¬ 
idium  this  can  be  understood  in  terms  of  the  high 
cross  sections  for  both  two-photon  resonance  ion¬ 
ization  and  laser-induced  Penning  ionization. 

A  representative  set  of  ionization  curves  is  pre¬ 
sented  in  Fig.  3.  The  original  density  was  as¬ 
sumed  to  be  10l”  (cm'1)  for  each  of  the  elements 
and  the  laser  irradiance  was  taken  as  10’  (Wcm*!). 
It  is  interesting  to  note  that  the  predicted  ioniza¬ 
tion  burnout  times  for  both  lithium  and  sodium 
are  shorter  than  for  potassium  and  cesium,  even 
though  the  seed-ionization  rates  for  potassium 
and  cesium  are  greater.  This  illustrates  that  a 
strong  intermediate  ionization  process  can  com¬ 
pensate  for  a  lower  seed-ionization  rate. 
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TABLE  I.  Basic  parameters  ,  ionization  times,  and  types  of  libors  solutions  for  alkali 
metals.*  SI:  /2>100*S,  S2:  /2^4RS,  and  S3:  72<4BS. 


Li 

Na 

K 

Rb 

Cs 

(CV) 

1.85  2.10 

1.62 

1.58 

1.43 

rj  (K) 

11  442  8433 

7757 

8036 

4163 

aA  (cm*l 

0  3.26  x  10"18 

1.27  xlO'17 

1.22  xio"17 

2.19  xlO"17 

(at  .Y0=  1016) 
aL  (cm*  sea 

9.1  xlO"44  2. 

,0  xlO"44 

1,2  xlO"43 

9.6  xlO"43 

4.5  xlO'42 

(7  2c  (cm4  sec) 

2.2x1  O'48  1. 

,7  x  10~49 

2.9  xlO"49 

7.7  xlO"48 

5.8  xlO'48 

1.0  *10"17  2, 

,9  xlO"18 

2.3  <10",s 

8.2  xlO"19 

4.5  xlO"1' 

Am2)b  (sec’1) 

3.5  xlO7  2.4  xlO7 

2.9  xlO7 

9.8  xio6 

9.8  xlO6 

(cm3sec"1> 

2.5  xio'7  1, 

.78  x  10“7 

1.5  xlO'7 

1.6  xio”7 

3.4  xlO'8 

K2m)  (cm3  sec"1) 

1.73  xlO"8  3.02  v  10“R 

8.3  xlO"8 

1.26  xlO’7 

1.42  xlO'7 

V0  !l 

(cm'3)  (H'cm'2) 

Ionization  time, 

rB  (nsec),  and  class  of  solution 

10t& 

105 

1060  (SI) 

837  (SI) 

994  (S2) 

976  (S2> 

2740  (S3) 

1015 

10* 

479  (SI) 

435  (Sli 

905  (S2) 

762  (S2) 

1940  (S2) 

1015 

107 

251  (SI) 

204  (SI) 

595  (S2) 

222  (S3) 

294  (S3) 

1015 

108 

124  (SI) 

93.0  (SI) 

131  (S3) 

27.2  (S3) 

31.8  (S3) 

10u 

10* 

180  (SI) 

124  (SI) 

92.3  (S2) 

87.3  (S2) 

171  (S3) 

1016 

io‘; 

84.7  (SI) 

80.3  (SI) 

91.3  (S2) 

81.4  (S2) 

138  (S3) 

1016 

107 

38.4  (SI) 

41.3  (SI) 

82.1  (S2) 

44.4  (S3) 

55.7  (S3) 

10hl 

108 

18.8  (SI) 

17.0  (SI) 

32.5  (S3) 

8.04  (S3) 

9.29  (S3) 

10" 

105 

23.2  (SI) 

13.9  (S2) 

8.56  (S2) 

7,63  (S2) 

10.6  (S3) 

1017 

106 

13.6  (SI) 

11.8  <S2) 

8.46  (S2) 

7.30  (S2) 

9.42  (S3) 

1017 

107 

6.39  (SI) 

7.63  (SI) 

7.83  (S2) 

5.46  (S2) 

5.50  (S3) 

1017 

108 

2.84  (SI) 

3.60  (SI) 

4.88  (S3) 

1.89  (S3) 

1.78  (S3) 

aThe  “intermediate  levels*’  used  correspond  to  the  next  four  levels  that  are  optically  con¬ 
nected  to  the  resonance  level . 

bThc  results  presented  in  this  table  assume  radiation  trapping  with  a  scale  length  of  0.1 
cm. 

In  Fig.  4  we  plot  the  variation  of  the  ionization  consistent.  This  is  reasonable  in  light  of  the 

burnout  time  against  the  laser  irradiance  for  an  fact  that  both  these  elements  have  lower  values  of 

initial  atom  density  of  101*  cm-*  for  all  of  the  al-  oA  and  vL  and  are  controlled  much  more  by  the 

kali  metals  considered  (Li  through  to  Cs).  It  is  intermediate  stage  of  ionization.  In  Fig.  5  we 

clear  that  in  the  case  of  potassium,  rubidium,  and  have  plotted  the  variation  of  the  ionization  burnout 
cesium,  the  laser  dependence  of  the  ionization  time  against  the  initial  atom  density  at  a  given 

burnout  time  rB  is  weak  at  the  lower  values  of  laser  irradiance  of  10'  (Wcm"2)  for  each  of  the 

laser  irradiance.  This  can  be  understood  in  terms  alkali  metals.  From  this  figure  it  is  clear  that 
of  the  dominance  of  associative  ionization.  Then  full  ionization  is  expected  for  both  lithium  and 

at  the  higher  values  of  laser  irradiance  the  in-  sodium  under  the  experimental  conditions  of  the 

fluence  of  laser-induced  Penning  ionization  (and  to  work  by  Lucatorto  and  Mcllrath.1  2  They  operated 

a  lesser  extent  two-photon  ionization)  becomes  with  a  laser-pulse  duration  of  about  500  nsec,  a 

apparent  by  the  rapid  fall  in  rB  with  peak  irradiance  of  about  10,;  (Wcm*2)  and  densities 

On  the  other  hand,  the  variation  of  tb  with  I1  for  6  x  101*  (cm"1)  and  10l,i  (cm-3)  for  lithium  and  sod  - 

lithium  and  sodium  can  be  seen  to  be  much  more  ium,  respectively.  On  the  other  hand,  appreciable 
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FIG.  3.  Ionization  time  histories  for  each  of  the  alkali 
metals  as  predicted  by  the  LI0ORS  simple  model.  The 
initial  number  density,  Af0  lO^om'3,  and  the  laser  ir- 
radianee  I1  10b  VV  cm"2. 


ionization  would  not  be  expected  under  the  experi¬ 
mental  conditions  of  Salter  of  al . 

The  ionization  times  predicted  in  Table  I  assume 
that  radiation  trapping  is  present  and  corresponds 
to  a  saturating  laser-beam  radius  of  0. 1  cm. 


FIG.  4.  The  variation  of  the  ionization  burnout  time 
r  with  laser  irradiance  for  each  of  the  alkali  metals 
at  an  Initial  density  of  1GU  cm*3  as  predicted  by  the 
MROKK  simple  model. 
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FKi.  f».  The  variation  of  the  ioniz.it  km  hurnoui  lime 
tb,  with  initial  density  at  a  laser  irradiance  of  10*  W 
cm*2  for  each  of  the  alkali  metals  as  predicted  by  the 
MBORS  simple  model. 

Although  this  is  quite  reasonable,  wo  have  con¬ 
sidered  the  influence  of  neglecting  radiation  trap¬ 
ping.  As  might  be  expected  the  importance  of  ra¬ 
diation  trapping  is  less  pronounced  at  high  values 
of  density  and  laser  irradiance.  This  can  he  seen 
by  reference  to  Table  II,  where  the  ionization 
times  are  recalculated  assuming  “no  radiation 
trapping”.  The  differences  in  the  ionization 
curves  trifh  and  trillion f  radio  lion  trapping  are 
illustrated  for  sodium  at  an  initial  density  of  101 
(cm"1)  in  Fig.  6. 

Although  we  have  attempted  to  take  account  of 
all  of  the  important  collisional-radiative  process¬ 
es,  ionization  resulting  from  the  collision  of  two 
excited  atoms,  where  one  is  in  the  resonance  state 
and  the  other  is  in  an  intermediate  state  has  not 
been  included  in  the  work  presented  in  this  paper. 
Such  collisions  could  have  a  large  cross  sec¬ 
tion3*’*4  and  could  be  fitted  into  our  model,  in¬ 
creasing  the  intermediate  ionization-rate  coef¬ 
ficient  /,  and  slightly  decreasing  the  burnout  rate 
coefficient  B.  Our  reason  for  not  including  them 
at  this  time  lies  in  a  lack  of  reliable  cross-section 
data. 

Nevertheless,  we  have  undertaken  a  preliminary 
parametric  study  using  a  range  of  possible  cross 
sections  in  the  case  of  cesium.  Our  results  sug¬ 
gest  that  the  ionization  time  could  be  considerably 
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TABLE  It.  Ionization  times  (nsec)  anti  types  of  LXBQRS  solution  for  each  of  the  alkali  metals 
assuming  no  radiation  trapping. 


N0  (cm*3) 

/*  (Wcm*2) 

Li 

Na 

K 

Rb 

Cs 

1015 

10s 

2430  (SI) 

1970  (SI) 

1720  (SI) 

1240  (SI) 

6440  (S3) 

1015 

10* 

809  (SI) 

755  (SI) 

1670  (SI) 

985  (SI) 

3340  (S2) 

1015 

lo7 

288  (SI) 

250  (SI) 

1160  (S2) 

342  (S2) 

764  (S3) 

1015 

108 

131  (SI) 

97.2  (SI) 

300  (S3) 

50.3  (S3) 

82.8  (S3) 

10’6 

10s 

229  (SI) 

159  (SI) 

112  (S2) 

92.0  (S2) 

255  (S3) 

10J* 

10* 

101  (SI) 

94.3  (SI) 

110  (S2) 

85.5  (S2) 

196  (S3) 

10}6 

107 

40.7  (SI) 

44.0  (SI) 

94.3  (S2) 

48.3  (S2) 

75.9  (S3) 

10“' 

I08 

18.9  (SI) 

17.5  (SI) 

41.5  (S3) 

9.00  (S3) 

12.7  (S3) 

1017 

10s 

24.0  (SI) 

14.6  (S2) 

8.77  (S2) 

7.65  (S2) 

12.2  (S3) 

10n 

10" 

14.0  (S2) 

12.2  (S2) 

8.70  (S2) 

7.35  <S2) 

10.8  (S3) 

1017 

107 

6.49  (SI) 

7.75  (SI) 

8.09  (S2) 

5.52  (S2) 

6.11  (S3) 

1017 

108 

2.86  (SI) 

3.64  (SI) 

5.05  (S3) 

1.92  (S3) 

1.96  (S3) 

shorter  than  we  currently  predict  for  cross  sec¬ 
tions  of  about  10~14  (cm2)  and  under  these  circum¬ 
stances  most  of  the  ionization  curves  switch  from 
S3  to  SI  solutions. 

SIMPLE-MODEL  COMPARISON  WITH  LIBORS 
COMPUTER-CODE  RESULTS 

As  a  further  check  on  the  reliability  of  our  simple 
model  of  LIBORS ,  we  have  compared  the  ioniza¬ 
tion  burnout  times  predicted  by  the  simple  model 
with  that  of  our  extensive  computer  code.1*'18 
This  comparison  has,  at  the  moment,  been  re¬ 
stricted  to  sodium  due  to  the  cost  and  time  in¬ 
volved  in  running  the  code  for  a  different  alkali. 

The  results,  as  evident  in  Fig.  7,  are  very  en¬ 
couraging  for  they  indicate  that  over  the  range  of 
densities  and  laser  irradiances  of  interest  in  the 
present  work  excellent  agreement  (^15%)  is  ob- 


FIG.  Pi.  Ionization  curves  for  sodium  with  and  without 
radiation  trapping,  as  predicted  by  the  LIBORS  simple 
model. 


tained.  Furthermore,  there  is  also  excellent 
agreement  in  the  detailed  features  of  the  ioniza¬ 
tion  time  histories.  This  can  be  gauged  by  com¬ 
paring  the  curve  in  Fig.  8  with  the  correspond 
ing  ionization  curve  for  sodium  in  Fig.  6.  Both 
display  the  same  shape  and  have  changes  of  slope 
at  about  the  same  values  of  Ne. 

It  is  also  worth  pointing  out  that  the  UBOR8 
computer-code  results,  as  shown  in  Fig.  8,  con¬ 
firm  that  our  basic  assumptions  regarding  the 
constancy  of  N2  and  Tg  over  the  major  portion  of 


FIG.  7.  A  comparison  of  the  ionization  burnout  time 
(s-95%)  variation  with  density  as  predicted  by  the  LIBORS 
simple  model  and  the  full  computer  code  for  sodium. 
Radiation  trapping  assumed  for  both  with  scale  length 
of  0.1  cm. 
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FIG.  8.  Time  histories  of  Se,  Te,  A2«  Q*  and  T{  as 
predicted  by  the  full  LlBORS  computer  code  (Refs.  15  to 
18)  for  sodium  with  an  original  density  of  101**  cm'3  and 
a  laser  irradiance  /*  10*  W  cm“2. 

the  ionization  time  are  reasonable.  Not  shown, 
but  in  fact  observed  in  our  UBORS  code  results  is 
the  proportionality  of  the  intermediate-state  pop¬ 
ulation  densities  and  the  electron  density  over 
much  of  the  ionization  period — corresponding  to 
Eq.  (18). 

CONCLUSIONS 

We  have  developed  a  relatively  simple  model  of 
laser  ionization  based  on  resonance  saturation. 
This  model  includes  essentially  all  of  the  import¬ 
ant  collisional  and  radiative  interactions  that  come 
into  play  when  an  alkali  metal  vapor  (density  range 
101  r'  to  1017  cm"1)  is  subject  to  a  sudden  pulse  of 
intense  laser  radiation  that  is  tuned  to  one  of  the 
resonance  transitions  of  the  alkali  atom.  This 
model  reveals  that  once  laser  saturation  of  the 
relevant  resonance  transition  has  been  achieved, 
ionization  proceeds  in  three  stages:  First,  seed 
ionization  (including  associative,  laser-induced 
Penning  and  multiphoton  ionization)  creates  an 
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initial  pool  of  seed  electrons.  Subsequently,  in 
the  intermediate  ionization  phase  these  free  elec¬ 
trons  are  heated  through  superelastic  collisions 
and  then  populate  intermediate  energy  states  which 
in  turn  are  photoionized  by  the  laser  radiation. 
Finally,  runaway  ionization,  through  direct  elec¬ 
tron-impact  ionization  of  the  intermediate-level 
population,  leads  to  a  very  rapid  and  near  com¬ 
plete  ionization  burnout  of  the  laser -irradiated 
species. 

We  have  also  been  able  to  show  that  this  model 
leads  to  two  types  of  analytical  solutions  for  the 
ionization  time  history  of  the  irradiated  vapor. 
Physical  interpretation  of  the  conditions  that  med¬ 
iate  which  type  of  analytical  solution  characterizes 
the  ionization  time  history  is  also  provided.  The 
ionization  time  for  each  of  the  alkali  metals  (Li 
through  to  Cs)  has  been  evaluated  under  a  wide 
range  of  conditions:  initial  number  densities  lO1^ 
to  1017  cm*1  and  laser  irradiances  from  10s  to  10R 
W  cm-2.  Although  to  date  there  has  been  very 
limited  experimental  work  on  this  new  class  of 
dense  plasma  production,  where  results  have  been 
available  they  are  consistent  with  our  predictions. 
We  have  also  been  able  to  show  that  the  predicted 
ionization  times  of  this  model  are  in  excellent 
agreement  with  the  results  of  our  comprehensive 
LEBORS  computer  code,  over  a  wide  range  of  den¬ 
sities  and  laser  irradiances ,  in  the  case  of  sod¬ 
ium.  In  general,  the  predictions  of  our  simple 
UBORS  model  should  only  be  regarded  as  a  guide 
to  the  expected  time  to  achieve  ionization  burnout 
since  the  spatial  and  temporal  characteristics  of 
the  laser  beam  will  not  in  reality  be  rectangular 
and  because  many  of  the  cross  sections  used  are 
not  known  with  great  precision.  Furthermore,  our 
neglect  of  continuity  is  not  valid  for  electron  den¬ 
sities  10f/<i  ,  and  this  leads  to  an  underestimate 
of  the  full  95'^ )  ionization  time  by  about  25'£ 
for  all  but  S3  solutions.  A  more  detailed  analysis 
of  this  will  shortly  be  published  by  us.35 
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We  have  developed  a  theoretical  model  of  the  ionization  that  results  from  extended  laser 
saturation  of  an  atomic  resonance  transition  within  a  dense  alkali  vapor.  Although  this  model 
includes  all  of  the  essential  collisional-radiative  interactions,  nevertheless,  it  leads  to  relatively 
simple  analytical  relations  that  predict  the  ionization  time  history  to  within  20%  of  that 
calculated  by  our  extensive  laser  ionization  based  on  resonance  saturation  (LIBORS)  computer 
code  in  the  case  of  sodium.  When  allowance  is  made  for  the  temporal  and  spatial  nature  of  a 
realistic  laser  pulse,  the  apparent  ionization  time  has  been  shown  to  be  approximately  double  that 
of  the  value  predicted  on  the  basis  of  a  steplike,  flat  laser  pulse. 

PACS  numbers:  52.50.Jm,  34.80. Dp,  32.80.Kf,  32.80.Fb 


INTRODUCTION 

The  concept  of  locking  the  population  of  two  atomic 
levels  in  the  ratio  of  their  degeneracies  by  an  intense  pulse  of 
laser  radiation  for  both  diagnostic  and  atomic  measurement 
applications  was  first  suggested  by  Measures. 1-3  Since  then 
the  field  of  laser  saturation  spectroscopy ,  as  it  has  become 
known  in  the  literature,  has  become  firmly  established  and 
its  range  of  applications  has  expanded  to  include  (i)  atomic 
lifetime  measurements,4**  (ii)  trace  element  analysis,7,8  (iii) 
laser  isotope  spearation,*  (iv)  single  atom  detection,10  (v) 
combustion  studies,1112  (vi)  flame  measurements,1'  (vii)  neu¬ 
tral  hydrogen  density  measurements  in  Tokamaks/4  M  and 
(viii)  fusion  plasma  impurity  diagnostics.16 

In  some  experiments 17~'*  the  laser  intensified  spontane¬ 
ous  emission  (referred  to  as  laser  induced  fluorescence  in 
much  of  the  literature)  was  not  observed  to  saturate  even  for 
values  of  laser  spectral  irradiancethat  far  exceed  the  saturat¬ 
ed  spectral  irradiance.  Recent  interferometric  measure¬ 
ments,20  however,  have  revealed  that  even  in  situations 
where  the  intensified  emission  does  not  appear  to  saturate, 
the  atomic  level  populations  are  close  to  the  infinite  tempera¬ 
ture  distribution  associated  with  a  population  ratio  equal  to 
the  degeneracy  ratio.  Rodrigo  and  Measures21  offered  one 
possible  explanation  of  this  paradox  in  terms  of  the  laser 
beam’s  radial  profile  and  the  observed  spatially  averaged 
emission  signal.  Recently,  Salter22  has  proposed  that  radi¬ 
ation  trapping  could  account  for  the  apparent  lack  of  satura¬ 
tion.  Some  combination  of  these  mechanisms  might,  in  fact, 
be  needed  to  adequately  account  for  these  anomalous 
observations. 

In  another  group  of  laser  saturating  experiments  per¬ 
formed  at  relatively  high  atom  densities  ( >  10'5  cm" 3)  sub¬ 
stantial  ionization  has  been  observed2’"26  when  the  laser 
pulse  is  tuned  to  a  resonance  transition  and  its  duration  is 
much  greater  than  the  resonance  state  lifetime.  Measures27 
was  the  first  to  recognize  that  such  extended  laser  saturation 
of  a  resonance  transition  represents  a  very  effective  means  of 
coupling  laser  energy  into  a  gaseous  medium.  In  essence  the 
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dense  population  of  resonance  state  atoms  serve  to  funnel 
laser  energy  directly  to  the  free  electrons  through  superelas¬ 
tic  collisions.  These  superelastically  heated  electrons  can 
then  readily  excite  and  ionize  resonance  state  atoms  since,  in 
effect,  these  atoms  have  their  ionization  energy  reduced  by 
the  laser  photon  energy.27,28  These  same  processes  have  re¬ 
cently  been  shown2*  to  account  for  the  fluorescence  ob¬ 
served  from  levels  that  lie  well  above  the  laser  pumped  level 
in  another  series  of  experiments.  '0-32 

The  initial  pool  of  free  electrons  in  all  of  these  experi¬ 
ments  is  expected  to  be  created  by  the  various  laser  induced 
collisional  processes  summarized  by  Geltman." 

We  have  developed  a  fairly  comprehensive  code  for 
studying  laser  ionizaton  based  on  resonance  saturation  (LI¬ 
BORS).  This  LIBORS  code  treats  the  laser  excited  atom  as  a 
20-energy  level  system  and  solves  the  subsequent  set  of  ener¬ 
gy  and  population  rate  equations  by  a  fourth-order  Runge- 
Kutta  technique.  A  detailed  description  of  the  code  is  pre¬ 
sented  elsewhere34'36  and  will  not  be  repeated  here.  It  will 
suffice  to  mention  that  the  code  evaluates  the  temporal  be¬ 
havior  of  (i)  each  of  the  20-atomic  level  populations,  (ii)  the 
free-electron  density,  (iii)  the  free-electron  temperature,  (iv) 
the  ion  temperature,  (v)  the  laser  power  absorbed  per  unit 
volume,  and  (vi)  the  power  radiated  per  second  per  unit  vol¬ 
ume  as  a  result  of  radiative  recombination  to  the  resonance 
level,  assuming  a  steplike  laser  pulse. 

The  results  of  our  computer  simulations  for  sodium  are 
consistent  with  the  experimental  observations  of  Lucatorto 
and  Mcllrath23  and  suggest  that  LIBORS  might  be  ideal  for 
creating  the  plasma  guide  channels  that  are  deemed  to  be 
necessary  for  transportation  of  electron  or  ion  beams  to  the 
fuel  pellet  of  future  particle  beam  fusion  reactors.34  We  have 
also  shown  that  the  concept  of  superelastic  laser  energy  con¬ 
version  (SELEC)  could,  under  certain  conditions,  have  sever¬ 
al  advantages  over  inverse  bremsstrahlung  as  a  means  of 
coupling  laser  energy  into  a  plasma  or  an  un-ionized  gaseous 
medium.  In  the  event  that  SELEC  is  applied  to  a  strong 
transition  within  a  suitable  ionic  species,  extremely  rapid 
rates  of  plasma  heating  can  be  achieved37 — a  feature  that  is 
quite  conducive  to  the  development  of  short  wavelength  la- 
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broadband  laser  radiation. 


sers.  Indeed,  in  the  case  of  a  boron  III  plasma,  we  predict37 
that  a  heating  rate  in  excess  of  10’ 3  °K  s" 1  could  be  attained 
with  a  laser  irradiance  of  close  to  10“  W  cm"2. 

In  the  present  paper  we  have  developed  a  more  accurate 
theoretical  model  than  has  been  previously  published38  and 
have  shown  that  it  can  be  used  to  predict  the  ionization  time 
history  for  alkali  vapors  exposed  to  a  realistic  laser  pulse, 
where  allowance  has  to  be  made  for  both  the  spatial  and 
temporal  characteristics  of  the  laser  pulse. 

LIBORS  SIMPLE  MODEL  FORMULATION 

We  shall  consider  the  case  of  a  homogeneous,  single 
constituent,  atomic  gaseous  medium  suddenly  irradiated 
with  an  intense  pulse  of  laser  radiation  that  is  tuned  to  one  of 
the  electronic  resonance  transitions.  Although  the  interac¬ 
tion  of  this  radiation  field  with  the  atoms  is  most  accurately 
treated  through  use  of  the  density  matrix  aproach,  a  rate 
equation  analysis  is  reasonably  reliable39  if  the  bandwidth  of 
the  laser  radiation  is  large  compared  to  the  Rabi  frequency, 
and  the  collisional  dephasing  time  is  short  compared  to  the 
times  of  interest.40  For  the  typical  conditions  found  in  LI¬ 
BORS  work  a  rate  equation  analysis  should  be  quite 
accurate. 

The  model  we  shall  develop  is  particularly  appropriate 
to  atoms  with  low-lying  resonance  states,  such  as  the  alkali 
metals.  For  this  class  of  atoms:  2  E2  \<EC  ,  5  3  E1U  where  E2  r 
(  =  hv ,  the  laser  photon  energy)  represents  the  energy  differ¬ 
ence  between  the  ground  and  the  laser  excited  resonance 
state  and  Ec ,  represents  the  ionization  energy  of  the  ground 
state.  If  the  laser  radiation  is  assumed  to  be  suddenly  applied 
at  t  =  0,  then  the  subsequent  chain  of  events  can  be  viewed  as 
proceeding  in  the  four  stages  indicated  in  Fig.  1 . 

According  to  Measures'  the  resonance  to  ground-state 
population  densities  are  locked  in  the  ratio  of  their  degener¬ 
acies,  viz.,  N2/N\  ^g2/g{^g\  in  a  characteristic  time, 
t*ss[(1  +  g)R21]  \  provided  / 1  (v)> Is (v)  in  the  case  of 


represents  the  stimulated  emission  rate  coefficient  for  the 
resonance  transition.  B2i  represents  the  appropriate  Milne 
coefficient  and  J’{  v)  the  corresponding  line  profile  function. 
Jl(v)  represents  the  laser  spectral  irradiance  and  7s(v)  is  the 
saturated  spectral  irradiance ,  given  by 

/»  =  [&irhvVc2(\  +g)]r?AD/r2, 

where  rf ad(=M  a,  '),  and  r2  is  the  resonance  state  lifetime. 
A 2 ,  is  the  resonance  transition  Einstein  probability.  Invari¬ 
ably,  tr  is  very  short  (much  less  than  1  ns  for  the  typical  laser 
fields  under  consideration)  compared  to  any  other  process 
and  we  assign  laser  saturation  to  Stage  1. 

Once  the  large  pool  of  resonance  state  atoms  is  created, 
several  kinds  of  interactions  give  rise  to  a  linear  growth  of  the 
free-electron  density  associated  with  Stage  2  The  most  im¬ 
portant  of  these  are  two-photon  ionization  of  the  resonance 
level ,  laser  induced  Penning  ionization ,  and  associated  ioniza¬ 
tion .  In  the  case  of  the  alkali  metals,  Geltman,33  de  Jong  and 
van  der  Valk,41  and  Klucharev42  have  provided  an  estimate 
of  the  magnitude  of  the  relevant  cross  sections.  The  appro¬ 
priate  rate  of  ionization  during^  Stage  2  can  be  expressed  in 
the  form 


dN{. 

”5T 


~N2</*F:  +  \N\(<tlvF  +  oAv), 

StaRc2 


(1) 


where  (cm4  s)  is  the  two-photon  resonance  state  ioniza¬ 
tion  rate  coefficient,  oL  (cm4  s)  is  the  laser  induced  Penning 
ionization  rate  coefficient,  v  (cm  s  1 )  is  the  mean  atom  veloc¬ 
ity,  F  (photons  cm  2  s  ')  represents  the  laser  photon  flux 
density,  i.e.,  F=  $Il(v)dv/hv%  and  aA  (cm2)  is  the  appropri¬ 
ate  associative  ionization  cross  section. 


FIG.  ]  Four  stages  of laser  ionization  based  on 
resonance  saturation.  STAGE  1.  Laser  rapidly 
locks  ground  and  resonance  level  populations  in 
ratio  of  dengeracies.  STAGE  2.  (i)  Rapid  growth 
of  free  electrons  due  to  two-photon  ionization  of 
resonance  level  and  laser-induced  Penning  ion¬ 
ization.  Associative  ionization  is  important  for 
some  elements,  (ii)  Free  electrons  rapidly  gain 
energy  through  superelastic  collisions.  STAGE 
3.  (i)  Direct  electron  impact  ionization  of  reso¬ 
nance  level  and  single-photon  ionization  of 
collisionally  populated  upper  levels  dominate 
the  rate  of  ionization,  (ii)  Electron  temperature 
stabilizes  as  the  rate  of  superelastic  heating  bal¬ 
ances  rate  of  collisional  cooling  through  excita¬ 
tion.  STAGE  4.  (i>  Runaway  collisional  ioniza¬ 
tion  of  upper  levels  occurs  once  a  critical 
electron  density  achieved,  (ii)  Superelastic  heat¬ 
ing  can  no  longer  balance  collisional  cooling 
and  electron  temperature  falls. 
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The  electrons  created  by  these  initial  processes  rapidly 
gain  energy  through  superelastic  collision  quenching  of  the 
laser  sustained  resonance  state  population,  as  first  suggested 
by  Measures.27 

In  Stage  3,  the  free-electron  temperature  is  viewed  as 
stabilizing  as  a  result  of  a  balance  between  the  rate  of  supere¬ 
lastic  heating  and  excitation  cooling.  This  will  be  considered 
in  more  detail  later.  These  electrons  now  give  rise  to  an  expo¬ 
nential  growth  in  the  free-electron  density  due  to  both  direct 
collisional  ionization  of  the  resonance  level  and  single-pho¬ 
ton  laser  ionization  of  the  collisionally  populated  intermedi¬ 
ate  levels. 

Finally,  in  Stage  4,  runaway  collisional  ionization  of  the 
intermediate  levels  occurs  once  a  critical  electron  density  is 
achieved.  This  process  leads  to  the  ionization  burnout  that 
results  in  almost  complete  ionization  of  the  laser  pumped 
species.  Our  full  LIBORS  computer  code34'36  indicates  that 
during  this  rapid  burnout  phase,  superelastic  heating  is  un¬ 
able  to  maintain  the  free-electron  temperature  and  a  sudden 
drop  in  the  temperature  is  expected. 

In  order  to  obtain  an  analytical  solution  for  the  ioniza¬ 
tion  time  history  from  our  LIBORS  simple  model  we  shall 
assume  that  the  free-electron  temperature  stabilizes  and  re¬ 
mains  constant  throughout  both  Stages  3  and  4.  Some  justifi¬ 
cation  for  this  assumption  can  be  found  in  the  following  ar¬ 
gument.  During  the  linear  growth  period  (Stage  2)  the 
dominant  ionization  processes  are  independent  of  the  free- 
electron  temperature.  This  is  fortunate  because  it  is  ques¬ 
tionable  whether  the  electron  energy  distribution  could  be 
characterized  by  a  temperature  during  this  period.  During 
Stage  3,  the  free-electron  energy  equation  can  be  reasonably 
well  expressed  in  the  form 

kTrXN2K2,E2t  -  NtKxiE2l  -  JV2  V  K2mEml,  (2) 

dt  rn.2 

where  Tr  is  the  free-electron  temperature,  Knm  is  the  appro¬ 
priate  electron  collision  rate  coefficient,  and  Em2  the  energy 
separation  between  the  intermediate  state  |  m  >  and  the  reso¬ 
nance  state  1 2  >  .  The  first  term  on  the  right-hand  side  repre¬ 
sents  the  superelastic  heating  rate  per  electron  while  the  sub¬ 
sequent  terms  represent  the  cooling  rate  per  electron 
associated  with  the  dominant  collisional  excitation  pro¬ 
cesses.  The  sum  over  m  in  the  last  term  extends  over  the 
group  of  intermediate  levels  that  are  rapidly  populated  by 
hot  electrons  and  depopulated  by  single-photon  laser  ioniza¬ 
tion.  Electron  heating  arising  from  the  seed  ionization  pro¬ 
cesses  is  assumed  to  be  of  secondary  importance  during 
Stage  3  and  4,  Also  the  gain  in  energy  due  to  single-photon 
ionization  and  the  loss  in  energy  from  direct  impact  ioniza¬ 
tion  of  the  resonance  level  are  assumed  negligible  in  com¬ 
parison  to  the  superelastic  heating  and  excitation  cooling 
terms  in  Eq.  (2). 

Using  electron  collision  rate  coefficients  that  are  based 
upon  Seaton  cross  sections,43  Eq.  (2)  yields  under  conditions 
of  laser  saturation  a  transcendental  equation  for  the  steady - 
state  superelastic  temperature 
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In  exp(  -  e ,/kT’t )  +  — S—  £  (f2m  )exp(  -  em/kT’t ) 

V 12/  3 


where -  Eit and  <fim)=fim «/  represents 

the  absorption  oscillator  strength  for  the  im  transition  and 
(G)im  the  relevant  Gaunt  factor  averaged  over  a  Maxwel¬ 
lian  velocity  distribution.  This  equation  can  be  solved  by  an 
iterative  procedure,  or  in  light  of  the  logarithmic  nature  of 
the  denominator  and  the  fact  that  Eml  zzEn  (in  the  model),  it 
can  be  approximated  by  the  expression 

t*s  _ _ Eyi/k _  ,a\ 


lnfl+gX/2.,//12) 

\  / 


The  ionization  equation  appropriate  to  this  model  takes 
the  form 

+  ^^+^1  +  NeK2c  ] 

at 

+  £A rm«F+V'Kmc),  (5) 

►n>3 

where  NmNeKmc  represents  the  volume  rate  of  direct  colli¬ 
sional  ionization  of  the  \m  >  state,  Nmoi^)eF  represents  the 
volume  rate  of  single-photon  laser  ionization  of  the  \m  > 
state,  Kmc  is  the  appropriate  electron  collisional  ionization 
rate  coefficient  of  |m  >  state,  and  o ^  is  the  single-photon 
ionization  cross  section44  for  the  |  m  >  state. 

The  relevant  intermediate  level  population  rate  equa¬ 
tion  can  be  expressed  in  the  form 

dN 

-  N2NcK2m  -  Nm  (<7^  +  A  l2  +  N,Cm  ),  (6) 

at 

where  A  fm2  =  Am2  ym2  represents  the  radiation- trapped  Ein¬ 
stein  transition  probability  for  the  m2  transition,  ym2  is  the 
Holstein  escape  factor,4*  and 

Cm^Kmc  4-  £Kmn  (7) 

n  #  m 

If  we  introduce 

(8) 

and  set 

p„=Pm  +Al2,  (9) 

then  over  a  significant  fraction  of  the  time  to  ionize 

pm>NeCm,  (10) 

as  at  low  densities  and  weak  fields  A  *,2  Cm  until  burnout 

phase,  while  under  intense  illumination  Pm  >NeCm .  Under 
these  circumstances  Eq.  (6)  has  a  solution 

NJt)  =  N2K2m  exp(  —  ^mt)J^iV,(»*)exp(^„t,)dt*,  (11) 

if  the  vapor  is  assumed  to  be  unexcited  prior  to  laser  irradia¬ 
tion,  i.e.,  Nm  (0)  =  0.  F ,  iV2,  and  K2m  are  also  taken  to  be 
constant  over  the  period  of  interest. 
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The  constancy  of  N2  over  an  appreciable  fraction  of  the 
time  to  ionize  is  quite  reasonable  on  the  grounds  that  the 
upper  levels  (m>3)  only  become  significantly  populated  just 
prior  to  ionization  burnout  (Stage  4}.  In  fact,  if  we  use  the 
saturation  condition  we  can  write  N2zzGNw  where 
G  =  g\  1  +  g)  1  and  Ao  represents  the  initial  atom  density. 
The  time  independence  of  Klm  is  based  upon  our  previous 
assumption  of  a  stabilized  temperature  during  the  interme¬ 
diate  phase  (Stage  3)  and  Fhas  been  set  to  be  constant  by  our 
assertion  of  a  steplike  laser  pulse.  The  viability  of  these  as¬ 
sumptions  is  apparent  from  our  LIBORS  computer  code 
results.'8 

Since  the  ionization  burnout  phase  under  the  condi¬ 
tions  of  interest  is  short  compared  to  the  time  to  attain  ion¬ 
ization  burnout  we  can  deduce  that  the  solution  prior  to 
burnout  will  take  the  form 

)~  A;  [exp09r)  -1].  (12) 

This  gives  a  linear  growth  for  0t  <  1  corresponding  to  the  seed 
ionization  phase  and  an  exponential  growth  for  0t  >  1  which 
takes  account  of  the  influence  of  the  intermediate  processes, 
i.e.,  Stage  3.  A  *  can  be  taken  to  be  the  seed  electron  density 
for  which  Ne  departs  significantly  from  a  linear  growth. 

Substituting  this  form  of  solution,  i.e.,  Eq.  (12)  in  Eq. 
(11),  yields  an  approximate  relation  for  the  intermediate  level 
population 

Nm[t)*N2K2mNe[t)/\p+pm),  (13) 

that  is  valid  fort  >(1/0  )ln(  1  •+•  0  /pm ).  This  constraint  is  dis¬ 
cussed  elsewhere.38  Equation  (13)  can  be  used  in  Eq.  (5)  to 
eliminate  Nm  and  yield  a  simple  form  of  ionization  equation 

=  S  +  INt+B\ l,  (14) 

dt 

where  we  have  introduced 

S  asATjogF2  +  +  oLF)\  (15) 

as  the  seed  ionization  rate. 


as  the  intermediate  ionization  rate  coefficient ,  and 


K  K 

m>10  +pm 

as  the  burnout  ionization  rate  coefficient. 
We  see  that  we  can  now  define 

a;w/s 


(17) 

(18) 


as  the  critical  electron  density  at  which  the  burnout  rate  of 
ionization  equals  the  intermediate  rate  of  ionization.  In  a  like 
manner  we  can  define 


a;  =s// 


(19) 


as  the  characteristic  electron  density  for  which  the  interme¬ 
diate  ionization  rate  equals  the  seed  ionization  rate. 

There  are  two  forms  of  analytical  solution  for  Eq.  (14) 
depending  on  the  relative  values  of  S,  /,  and  B.  If  / 2  <  4SB, 
then  this  is  equivalent  to  A  **  <  4A  *  and  so  the  ionization 
proceeds  directly  from  Stage  2  to  Stage  4.  That  is  to  say  seed 


ionization  leads  directly  to  runaway  ionization.  The  solution 
to  Eq.  (14)  under  these  circumstances  takes  the  form 

2/  ,2  BNC  +  I  t\ 


a  =  (4SB-l2)1'2.  (21) 

This  class  of  solutions  is  in  general  indicative  of  a  high  rate  of 
seed  ionization. 

If  / 2  >  4 SB,  then  A  *#  >  4 A  *  and  Stage  3  ionization  will 
play  some  role  in  separating  Stage  2  from  Stage  4.  The  solu¬ 
tion  to  Eq.  ( 14)  under  these  circumstances  takes  the  form 


=  1  infill!  *4 

A  UN,  +  R2)  Rj 


A  =  (/2  —  4SB  )I/J, 
Rt  =  U-A)/2B, 


R2  =  (/  +  /!  )/2 B.  (25) 

From  Eq.  ( 14)  it  is  apparent  that  /  is  associated  with  the 
exponential  growth  phase  and  we  therefore  identify  0  with  I. 
Under  these  circumstances  Eq.  (16)  represents  a  transcen¬ 
dental  equation  for  /.  We  have  found  that  /  can  be  deter¬ 
mined  to  about  1  %  accuracy  by  setting/?  m  equal  to  the  maxi¬ 
mum  value  of  pm  in  Eq.  (16),  then  solving  the  resulting 
quadratic  equation  for  /  and  using  this  approximate  value  of 
I  in  the  denominators  of  the  sum  on  the  right-hand  side  of 
Eq.  (16).  The  more  accurate  value  of  /is  then  substituted  for 
0  in  Eq.  (17)  to  ascertain  B. 

A  pair  of  representative  solutions  for  sodium  and  potas¬ 
sium  are  presented  (broken  curves)  in  Fig.  2.  Although  the 
ionization  time  histories  and  burnout  times  predicted  from 
this  analysis  are  expected  to  provide  a  reasonable  compari- 


0  75  150 


TIME  (nsec) 

FIG.  2.  Representative  ionization  time  histories  with  and  without  the  bur¬ 
nout  solution  for  sodium  and  potassium  with  initial  atom  density  ty,  =  10,ft 
cm  '  and  a  steplike  laser  irradiance,  / '  =  10*  W  cm 
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son  between  the  elements  and  to  indicate  the  dependence  of 
these  features  on  the  laser  irradiance  and  initial  atom  densi¬ 
ty — the  form  of  the  ionization  equation  ( 14)  has  one  draw¬ 
back.  It  does  not  take  into  account  continuity  and  therefore 
predicts  runaway  ionization  leading  to  a  singularity  in  jV  .  In 
reality  this  is  not  too  serious  as  the  burnout  phase  is  very 
rapid,  in  all  cases,  so  that  the  loss  in  accuracy  in  predicting 
the  ionization  burnout  time  from  Eqs.  (20)  or  (22)  by  setting 
ymax  _  ^  the  jnitiai  atom  density,  is  minimal. 

A  more  reliable  total  solution  can  be  obtained  by  treat¬ 
ing  the  burnout  phase  separately.  In  essence  the  burnout 
phase  can  be  regarded  to  arise  when  either  Nt.  /B  (for 
1 2  >  )  or  A.  ~  (5  /B  ) ' 1 1 2  (for  1 :  <  455 ).  At  these  high  elec¬ 

tron  densities,  the  rate  of  ionization  is  dominated  by  electron 
collision  ionization  and  single-photon  ionization  of  the  inter¬ 
mediate  level  populations.  In  some  cases  electron  collisional 
ionization  of  the  resonance  level  may  also  be  significant. 
Thus  the  burnout  ionization  equation  can  be  expressed  in  the 
form 

dN 

—j—  =  A r2N,K2l  +  £  Nm  ).  (26) 

Ut  m.,11 

In  the  burnout  phase  electron  collisions  will  tend  to  bring  the 
intermediate  level  populations  into  local  thermodynamic 
equilibrium  (LTE)  with  the  resonance  level  population. 
Hence  we  can  write  for  the  intermediate  level  population 

Nm  =  N2{gm/g2)  exp(  -  Em2/kTe).  (27) 

Continuity  can  be  taken  into  account  by  assuming  that 


N'  +*2  +  JV,+  ^Nm=Na  (28) 

m  s  3 

applies  in  the  burnout  phase.  If  we  also  assume  laser  satura¬ 
tion  of  the  resonance  transition,  then  the  continuity  relation 
can  be  rewritten  in  the  form 

^  =  ii-AU  (29) 

where 


1  +  SlSm/Sslex P(  -  Em2/kTc 


(30) 


If  Eq.  (29)  is  used  with  Eq.  (27)  in  Eq.  (26),  the  burnout  phase 
ionization  equation  can  be  written 


dN. 

dt 


where 


and 


H  = 


r  V 


\F  exp- 


(31) 


+  (32) 

V  kT  ) 


(33) 


~*gz  kT 

Equation  (31)  has  a  simple  solution  of  the  form 

,  .  ,  (  1  ,  KNt  +Hy. 

t  =  f,  4-  I - In — - )  ,  (34) 

\{xNa  +  H)  Nn  —  Ne  )s[  '  ’ 

where  t ,  is  the  time  at  which  the  electron  density  reaches  the 


TABLE  I:  Basic  parameters  and  95%  ionization  times  calculated  from  LIBORS  model  for  several  alkali  metal  vapors. 


Li 

Na 

K. 

Rb 

E2{ 

(eV) 

1.85 

2.10 

1.62 

1.58 

n 

CK) 

11442 

8433 

7757 

8036 

*4 

(cm2) 

(at  N0  =  10lft) 

0 

3.26X10  ,B 

1.27x10  ” 

1.22x10  17 

<*L 

(cm4  s( 

9.1X10  44 

2.0x10  44 

1.2X10' 43 

9.6X10  41 

< 

(cm4  s) 

2.2X10  4* 

1.7x10  4* 

2.9  x  10  ** 

7.7X104" 

(cm2) 

1.0X10  17 

2.9X10  ,H 

2.3X10' 

8.2X10'"’ 

ml  > 

(S  '( 

3.5X107 

2.4  X107 

2.9  X107 

9.8  X  10* 

(cm’s  ') 

2.5X10  7 

1.78x10  7 

1.5X10  7 

1.6X10-’ 

<Kim) 

(cm*  s  ‘) 

1.73X10  " 

3.62X10  * 

8.3X10  K 

I.26X10-’ 

V„ 

r 

95%  Ionization  Time  (ns) 

(cm  ’) 

[Wcm  J) 

10'5 

10* 

1120 

972 

1200 

1140 

!0'5 

to* 

521 

539 

1080 

914 

10” 

107 

264 

262 

707 

343 

10'5 

10* 

125 

102 

169 

84.4 

10,ft 

10* 

188 

140 

113 

104 

10,ft 

10* 

90.9 

93.8 

112 

98.0 

10lh 

107 

42.3 

51.7 

100 

60.2 

10” 

10" 

20.0 

22.8 

44.3 

20.2 

10'7 

10* 

24.2 

15.7 

10.6 

9.28 

I0’7 

106 

14.3 

13.4 

10.6 

8.98 

10'7 

107 

7.00 

8.98 

9.88 

7.14 

10'7 

10" 

3.26 

4.65 

6.68 

3.46 

'The  results  presented  in  this  table  assume  radiation  trappping  with  a  scale  length  of  0.1  cm. 

"The  intermediate  levels  chosen  for  these  calculations  correspond  to  the  next  four  levels  that  are  optically  connected  to  the  resonance  level. 
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FIG.  3.  Variation  of  the  95%  ionization  time  with  laser  irradiance  for  lith¬ 
ium,  sodium,  potassium,  and  rubidium  assuming  a  steplike  laser  pulse  of 
extended  duration. 

critical  *ue  A  J,  for  switching  solutions.  A  'e  is  chosen  to  be 
the  electron  density  at  which  the  slope  of  Ne(t ),  correspond¬ 
ing  to  Eq.  ( 1 4),  becomes  steeper  than  the  slope  corresponding 
to  Eq.  (31).  This  procedure  results  in  a  smooth  transition  as 
seen  by  reference  to  Fig.  2  where  two  representative  electron 
density  time  histories  are  presented,  with  and  without  bur¬ 
nout  solutions ,  for  both  sodium  and  potassium.  Similar 
curves  have  been  obtained  for  lithium  and  rubidium  under 
the  same  conditions  of  density  and  laser  irradiance. 

The  basic  parameters  and  the  time  to  achieve  95%  ion- 
izatijn ,  evaluated  from  the  above  analysis  is  displayed  in 
Table  I  for  lithium,  sodium,  potassium,  and  rubidium.  As 
might  be  expected,  these  95%  ionization  times  are  somewhat 
longer  than  the  ionization  times  evaluated™  solely  on  the 
basis  of  Eq.  (14),  but  for  the  most  part  the  difference  is  within 
30%.  The  variation  in  the  95%  ionization  time  as  a  function 
of  laser  irradiance  at  an  initial  atom  density  of  1016  cm-3  is 
presented  in  Fig.  3  for  Li,  Na,  K,  and  Rb.  In  the  case  of 
cesium,  the  combination  of  a  predicted  low  electron  tem¬ 
perature  and  a  large  associative  ionization  cross  section 
leads  to  a  situation  where  the  slopes  cross  at  an  electron 
density  that  is  too  small  to  justify  the  use  of  the  burnout 
solution. 

SIMPLE  MODEL  COMPARISON  WITH  LIBORS 
COMPUTER  CODE  RESULTS  FOR  SODIUM 

In  order  to  test  the  accuracy  of  our  simple  model  of 
LIBORS  we  have  made  a  direct  comparison  between  the 
electron  density  time  histories  predicted  by  our  LIBORS 
computer  code34'™  and  that  predicted  by  the  model  de¬ 
scribed  herein.  This  comparison  has,  at  the  moment,  been 
restricted  to  sodium  due  to  the  cost  and  time  involved  in 
running  the  full  code  for  a  different  alkali.  The  results  for 
Art  =  I015  and  I0,6cm  3  and  /'  =  10*  and  10*  W  cm'  2  are 
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presented  in  Fig,  4.  Both  linear  and  logarithmic  plots  are 
used  in  order  to  display  the  degree  of  agreement  over  a  range 
of  electron  density  spanning  six  orders  of  magnitude. 

It  is  clear  that  the  simple  LIBORS  model  gives  excellent 
agreement  with  the  full .  jBORS  computer  code  over  a  wide 
range  of  conditions.  1  he  e;fferences  in  the  95%  ionization 
times  are,  in  general,  within  20%,  while  the  times  to  achieve 
only  10%  ionization  are  much  closer.  What  is  particularly 
worthy  of  note  is  the  close  overlap  of  the  ionization  histories 
over  six  orders  of  magnitude  variation  in  the  electron  densi¬ 
ty.  Clearly  the  LIBORS  simple  model  can  be  viewed  as  con¬ 
taining  the  essence  of  the  complex  array  of  collisional-radia- 
tive  processes. 

ALLOWANCE  FOR  SPATIAL  AND  TEMPORAL 
CHARACTERISTICS  OF  THE  LASER  PULSE 

In  order  to  be  able  to  make  a  comparison  between  the 
simple  model  description  of  LIBORS  and  experimental  ob¬ 
servations,  it  is  immediately  apparent  that  we  have  to  take 
account  of  both  the  radial  profile  of  the  laser  pulse  and  its 
temporal  history. 

The  above  analysis  has  been  predicated  on  the  basis  of  a 
step  function  form  of  laser  pulse,  that  is  to  say  the  laser  irra¬ 
diance  was  assumed  to  be  instantly  switched  to  its  maximum 
value  and  remain  at  that  value  indefinitely.  In  reality  we 
have  found  in  a  preliminary  experiment46  that  a  typical 
flashlamp  dye  laser  pulse  can  be  well  represented  by  the 
expression 

/'(*)» U  X  107£'(t  /tt)  (1  —  r  /r,)exp(  -  t  /r2),  (35) 

where  I*(t)  represents  the  laser  irradiance  in  W  cm  " 2  and  E 1 
the  total  energy  in  the  laser  pulse  in  J.  ry  ( ~  160  ns),  r, 

( ~  1154  ns),  and  r2  ( »  8 10  ns)  represent  the  three  character¬ 
istic  times  needed  to  define  the  shape  of  the  pulse. 

The  temporal  variation  of  the  free-electron  density  cor¬ 
responding  to  the  laser  pulse  described  by  Eq.  (35)  can  be 
evaluated  by  replacing  this  laser  pulse  by  a  series  of  represen¬ 
tative  rectangular  pulses  and  integrating  the  appropriate 
ionization  equation  over  each  time  interval.  The  lower  limit 
for  the  electron  density  in  each  step  is  set  equal  to  the  maxi¬ 
mum  value  achieved  in  the  previous  integration  interval. 
The  initial  time  for  the  integration  is  obtained  by  displacing 
the  time  origin  by  an  amount  equal  to  the  difference  in  the 
time  required  to  reach  the  starting  electron  density  along  the 
ionization  curve  corresponding  to  the  previous  value  and  the 
new  value  of  laser  irradiance. 

In  order  to  predict  the  temporal  variation  of  the  radia¬ 
tive  recombination  emission  the  free-electron  line  density , 

.  *  „(/)=  f  Ne\r,t)2nrdr ,  (36) 

Jo 

within  the  effective  field  of  view  of  our  photodetection  sys¬ 
tem46  (radius  R  )  has  to  be  evaluated  as  a  function  of  time. 
This  is  accomplished  by  using  cylindrical  geometry  and  di¬ 
viding  the  observed  alkali  vapor  column  segment  into  a  se¬ 
ries  of  equal-area  annular  zones  colinear  with  the  laser  beam 
axis.  Each  zone  is  assumed  to  be  exposed  to  a  laser  pulse 
having  a  temporal  variation  given  by  Eq.  (35).  The  laser  pulse 
is  also  assumed  to  have  a  Gaussian  radial  profile  of  the  form 
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FIG  4.  Comparison  of  the  ionization  time  histories  (on  both  linear  and  logarithmic  scales)  as  calculated  from  the  simple  LIBORS  model  and  as  evaluated 
from  the  full  LIBORS  code  for  a  range  of  conditions  in  the  case  of  sodium. 
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FIG  5  Si  hematic  rept osculation  of  the  lom/ation  in  both  space  and  time 
indicating  both  the  temporal  and  spatial  nature  of  the  laser  pulse 

J  Appl  Phys  Vol.  52.  No.  3.  March  1981 


/  I  =  /'(*  )exp[  -  (r/r0):],  (37) 

where  1 '(/ )  is  given  by  Eq.  (35)  and  r0  represents  the  expo¬ 
nential  radius  for  the  laser  pulse.  Equation  (36)  is  evaluated 
using  a  Simpson’s  integration  technique  where 
/  '{rj )  and  )  are  determined  at  the  mean  radius  and  at 
the  two  boundary  radii  of  each  annular  zone.  In  order  to 
facilitate  computation  we  assumed  that  the  electron  density 
froze  at  either  the  burnout  value  or  the  value  ihal  cauv.v’ 
breakdown  of  the  laser  saturation  condition.  A  schema.  • 
representation  of  the  variation  of  ionization  in  both  space 
and  time  allowing  for  the  spatial  and  temporal  nature  of  the 
laser  pulse  is  presented  in  Fig.  5 

In  Fig.  6,  we  present  a  comparison  between  the  tempo¬ 
ral  variation  of  the  free-electron  line  density,  1  t  evaluated 
to  a  radius  of  0.6  cm)  for  the  three  eases  (i).a  steplike  laser 
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FIG  6.  Comparison  of  the  temporal  variation  of  the  free  electron  line  densi¬ 
ty.  ,  in  a  sodium  column  evaluated  to  a  radius  of  0.6  cm  for  the  three  cases, 
tit  Steplike,  flat  laser  pulse  for  / 1  ~  10*  W  cm  \  (id  Laser  pulse  is  assumed 
to  have  realistic  temporal  variation.  Eq.  (35)  of  text,  but  uniform  radial 
distribution.  Energy  in  the  pulse  is  set  to  be  0.227  J  corresponding  to  a  peak 
irradiance  of  10*  W  cm  :.  (lit)  Laser  pulse  is  assumed  to  have  both  a  realistic 
temporal  variation  and  a  Gaussian  radial  distribution  with  an  energy  of 
0.227  J 

pulse  that  also  has  a  flat  radial  distribution  so  that  /  ;  =  10* 
W  cm  2  for  t>0  and  0  <  r<0.6  cm.  (ii)  A  laser  pulse  having  a 
realistic  time  history,  given  by  Eq.  (35),  but  a  flat  radial  dis¬ 
tribution.  The  energy  of  the  laser  pulse  is  assumed  to  be 
0.227  J  which  corresponds  to  a  peak  laser  irradiance  of  10* 
W  cm  (iii)  A  realistic  laser  pulse  having  a  time  history, 
given  by  Eq.  (35),  and  a  Gaussian  radial  distribution,  given 
by  Eq.  (37),  with  an  exponential  radius  of0.2  cm.  The  energy 
of  the  laser  pulse  is  again  assumed  to  be  0.227  J. 

The  maximum  free-eleclron  line  density ,  /  *)a'  and  corre¬ 
sponding  maximum  free-electron  density,  averaged  over  the 
column  radius,  (A^)max  appropriate  to  each  of  these  three 
cases  is  presented  in  Table  II. 

It  is  clear  from  Fig.  6  that  in  an  experiment  the  apparent 
time  to  achieve  95%  ionization  would  be  about  double  that 
predicted  on  the  basis  of  our  simple  model. 

CONCLUSIONS 

We  have  developed  a  theoretical  model  of  the  ionization 
that  arises  when  an  alkali  metal  vapor  is  suddenly  irradiated 
by  an  intense  pulse  of  laser  radiation  that  is  tuned  to  one  of 
the  alkali  resonance  transitions.  This  model  includes  essen¬ 
tially  all  of  the  important  collisional  and  radiative  interac¬ 
tions.  Once  saturation  of  the  resonance  transition  is  achieved 
ionization  is  assumed  to  proceed  in  three  stages:  (i)  seed  ion¬ 
ization  creates  an  initial  pool  of  free  electrons  that  are  free  to 
gain  energy  through  superelastic  collisions ,  (ii)  these  hot  elec¬ 
trons  are  then  capable  of  directly  ionizing  the  resonance  level 
atoms  and  populating  the  intermediate  atomic  levels,  which 
are  subsequently  photoion  ized ,  and  (iii)  direct  electron  impact 
ionization  of  the  intermediate  level  population  subsequently 
leads  to  a  very  rapid  and  near  complete  ionization  burnout  of 
the  laser  irradiated  species. 

Analytical  solutions  for  the  ionization,  derived  from 
this  relatively  simple  model,  are  found  to  lead  to  free-elec¬ 
tron  density  time  histories  that  are  in  very  close  agreement 


TABLE  II.  Sodium  .V()  =  10'*  cm  \ 


Case 

Intern  ') 

(cm  ') 

(II 

1.13  x  10,f’ 

10‘* 

HI) 

1.13  x  lO1* 

10u‘ 

(III) 

6x10** 

5.3x10" 

with  those  obtained  on  the  basis  of  our  detailed  LIBORS 
computer  code  in  the  case  of  sodium.  The  model  has  also 
been  applied  to  other  alkali  metal  vapors  and  the  predicted 
time  to  achieve  95%  ionization  has  been  tabulated  over  a 
wide  range  of  conditions. 

In  an  experiment,  the  laser  pulse  will  not  necessarily 
have  a  steplike  time  history  and  a  flat  spatial  distribution,  as 
assumed  in  the  simple  LIBORS  model.  We  have  shown  that 
the  apparent  time  to  achieve  95%  ionization  may  be  nearly 
twice  that  predicted  by  our  simple  model,  if  allowance  is 
made  for  the  temporal  and  radial  nature  of  an  experimental 
laser  pulse. 
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Erratum:  Theoretical  model  of  laser  ionization  of  alkali  vapors  based  on 

resonance  saturation 

[J.  Appl.  Phys.  52, 1296  (1981)] 

R.  M.  Measures,  P.  G.  Cardinal,  and  G.  W.  Schinn 

Institute  for  Aerospace  Studies.  University  of  Toronto.  4925  Dufferin  Street,  Downsview,  Ontario .  Canada 
PA CS  numbers:  52.50.Jm,  34.80.Dp,  32.80.Kf,  32.80.Fb,  99.10.  +  g 


In  the  above  paper.  Fig.  6  and  Table  II  were  incorrect. 
The  correct  Fig.  6  and  Table  II  are  provided  herein.  A  small 


TABLE  II.  Sodium  N0  =  10'*  cirT3. 
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FIG.  6.  Comparison  ofihe  temporal  variation  of  ih tfree  electron  line  densi¬ 
ty.  I  \  in  a  sodium  column  evaluated  to  a  radiuso!T).6cm  for  the  three  cases. 
|t)  Steplike,  flat  laser  pulse  for  / '  =  10*  W  cm“2.  ill)  Laser  pulse  is  assumed 
to  have  realistic  temporal  variation.  Eq.  (35i  of  test,  but  uniform  radial 
distribution  out  to  0.6  cm.  Energy  yf  the  pulse  is  set  to  be  0.808  J  correspond¬ 
ing  to  a  peak  irradiance  of  10*  W  cm*2.  (Ill)  Laser  pulse  is  assumed  to  have 
both  a  realistic  temporal  >  ariation,  Eq.  (35),  and  a  Gaussian  radial  distribu¬ 
tion  \*iih  an  exponential  radius  of  0.20cm  and  an  energy  of  0.091  J,  giving  a 
peak  irradiance  of  10*  W  cm  , 


Case 

(I) 

1.13x10"“ 

ur 

(II) 

U3X10"* 

10" 

(HI) 

4.96 X  10“’ 

4.39  X  10'5 

section  of  the  text  relating  to  Fig.  6  and  Table  II  was  also 
incorrect.  The  correction  reads: 

In  Fig.  6,  we  present  a  comparison  between  the  tempo 
ral  variation  of  the  free-electron  line  density,  (evaluated 

to  a  radius  of  0.6  cm)  for  the  three  cases:  (I)  A  steplike  laser 
pulse  that  also  has  a  flat  radial  distribution,  so  that  I1 
=  106  W  cm'2  for  f>0  and  0  < r^0.6  cm.  (II)  A  laser  pulse 
having  a  realistic  time  history,  given  by  Eq.  (35),  but  a  flat 
radial  distribution  out  to  0.6  cm.  The  energy  pulse  of  the 
laser  is  assumed  to  be  0.808  J,  which  corresponds  to  a  peak 
laser  irradiance  of  106  W  cm'2.  (Ill)  A  realistic  laser  pulse 
having  a  time  history,  given  by  Eq.  (35),  and  a  Gaussian 
radial  distribution,  given  by  Eq.  (37),  with  an  exponential 
radius  of  0.2  cm.  The  energy  of  the  laser  pulse  is  assumed  to 
be  0.09  iJ. 

It  might  be  pointed  out  that  these  errors  do  not  in  any 
way  modify  our  conclusions  as  the  change  in  the  effective 
ionization  time  for  the  case  of  the  realistic  laser  pulse  is  still 
close  to  double  that  obtained  with  a  steplike  pulse. 
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Abstract 

Laser  resonance  saturation  of  alkali  vapours  represents  a  very  attrac¬ 
tive  net  hod  of  creating  the  extended  plasma  channels  of  interest  in  light 
ion  he  a n  fusion.  We  have  undertaken  a  preliminary  study  of  the  influence 
of  molecular  nitrogen  upon  this  laser  ionization  technique.  Our  results 
indicate  that  the  electron  cooling,  resonance  quenching  and  the  increase 
in  the  laser  energy  requirement  are  acceptable,  providing  that  the  alkali 
seeding  exceeds  0,5u» 


*Work  supported  by  the  U.S.  Sandia  National  Laboratory  under  Contract  No 
U9-2859A  and  a  grant  from  the  Natural  Science  and  Engineering  Research 
Council  of  Canada, 
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rHTROPtCTION 


A  fusion  reactor  baaed  upon  light  ion'  bean  drivsrs  prcaises 

to  be  less  expensive  and  more  efficient  tt  x  other  inertial  confinement 
reactor  schemes.  The  light  ion  bean  app;  cw  ch  to  fusion  also  lends  itaelf 
to  nethoia  of  minimizing  the  vail  loading  of  aucb  a  reactor  by  ezpl  crying  a 
large  reactor  cavity  and  using  an  atmosphere  that  will  absorb  ao«t  of  the 
energy  liberated  (other  than  that  released  In  the  for*  of  fast  neutrem) 
when  the  fuel  pellet  undergoes  fusion.  A  a  (kind  off  distance  of  several 
ns  ter s  nay  be  required  la  order  to  protect  the  array  of  loo  diode a  free 
the  ndcroexploalon  and  this  vlll  necessitate  tbs  development  of  an  efficient 
mode  of  ion  beam  transport  across  the  reactor  chamber. 

There  is  a  reasonable  body  of  experience  that  has  demonstrated  that 
light  Ion  beams  can  be  efficiently  transported  for  distances  of  at  least 
one  mater  by  means  of  current  carrying  plasm*  channels. ^~3^  These  plasma 
channels  provide  both  charge  aai  beam  current  neutralization  and  sufficient 
azimuthal  magnetic  field  to  confine  and  guide  the  loos.  Recent  experiments 

have  also  ahovn  that  these  discharge  channels  can  be  initiated  and  gulled 
fit 

by  laser  heating  of  a  molecular  gas.'  '  '  In  this  work  a  COg  laser  la  used 
to  vib rationally  excite  ammonia  molecules  thereby  creating  a  low  density, 
very  weakly  ionized  channel  for  gulling  the  discharge.  Although  this 
approach  may  be  adequate  for  testing  the  concept  it  is  unlikely  to  be  * 
candidate  for  a  fusion  pever  reacto-  where  the  antolent  tenperature  could 
be  in  excess  of  1000 ‘R. 

We  have  proposed^  that  saturating  a  strong  transition  within  one  of 
the  gaseous  constituents  (such  as  lithium  vapour)  of  a  fusion  reactor  with 
an  appropriately  tuned  laser  represents  an  almost  ideal  method  of  creating 
the  multiple  electrical  discharge  guide  paths.  In  regard  to  this  task  the 
advantages  inherent  in  our  proposed  Icser  ionisation  tassd  on  vsscnionaa 
saturation  -  LBCRS  technique  tret 

(1)  Bigh  EfficitKOii  due  to  the  large  cross  sections  Involved  in  directly 
converting  the  laser  energy  Into  Ionization  and  electron  energy. 

(2)  Complete  Ionisation  {>  95"0  of  laser  pua^>ed  species  along  the  path  of 
the  laser  bean,  even  for  relatively  modest  values  of  laser  irrtliance. 

(3)  Uniform  Plasm  Charms l»  s'-o-li  be  formed  due  to  the  linear  rate  of 
energy  deposition  along  the  path  of  the  laser  beam. 

(4)  Super* l as etc  Hiatinj  of  t*t  free  electrons  to  provide  high  conduc¬ 
tivity  channels. 


^  to  the  formation  of  a  plasm*  channel,  with  an  electron  density  of  close 
to  1015  cm  3,  In  *  time  ^  le,,  than  1  nsec  in  the  presence  of  mlecular  nitro¬ 
gen  at  a  density  of  at  least  2  x  1017  cm’3,  in  audition  we  have  evaluated  the 
lncrmaae  la  the  mini**  laser  energy  required  for  the  creation  of  such  channels 
by  this  background  gas.  Where  possible  we  have  made  slxpUryirg  as  suctions 
that  are  strongly  conservative  sad  we  therefore  expect  that  our  results  will, 

If  anything,  probably  overpredict  the  laser  energy  required  to  create  these 
P1***  channels.  Even  so,  our  analysis  is  encouraging  for  it  appears  that 
UB<*S  still  represents  an  efficient  and  low  energy  method  of  produclj* 
the  discharge  guide  path  necessary  for  ion  beam  transport . 
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(5)  R’lpu!  Formation  of  the  plazra  channels  through  exponential  grewth  of 
ionization. 

(6)  Itiniizc  Fellst  Frsiisat  by  the  laser  radiation  due  to  the  modest  demand 
of  laser  energy. 

(T)  Cbrjtzttbiiiijf  of  the  constituent  (such  as  lithium  vapour),  that  is 
•unable  to  LI30'riS,  with  the  p-oposed  operating  conditions  within  a 
fusion  reactor , 

(8)  Resonor.-e  tins  laesro  currently  under  development'''  could  lead  to  a 
alnplification  of  the  design  and  operation  of  the  lasers  required  for 
LIBORS. 

In  our  initial  etudy^  we  have  ahcnm  that  a  laser  pulse  of  less  than 
1  Joule  would  be  required  to  Creite  a  plasma  channel  with  an  electron  density 
of  close  to  10*^  cm'3  over  a  5m  path  in  the  case  of  a  0.1  torr  aodlua  atmos¬ 
phere,  Furthermore,  we  can  say,  baaed  on  our  recent  development  of  a  aisg>le 
rndel  of  LIB0RS, that  this  estimate  would  also  apply  to  a  similar  lithium 
atmosphere.  This  can  be  seen  by  reference  to  figure  1  where  the  variation  in 
the  95<(  ionization  burnout  time  with  laser  i  r  radiance  is  presented  for  lit  hi  urn, 
sodium,  potassium  and  rubidium.  Although  these  ionisation  times  were  estimated 
on  the  basis  of  a  spatially  uniform,  step- like  laser  pulse,  we  have  recently 
sbovn^  that  in  the  case  of  sodium,  allcwance  for  a  realistic  laser  pulse 
(both  spatially  and  t amorally)  should  only  lead  to  a  doi&llng  of  the  obssrvsd 
Ionization  time.  By  way  of  illustration,  in  figure  2,  we  present  the  tlms 
history  of  the  growth  in  the  frss  slsatron  line  density,  <>V^(t)  for  three 
kinds  of  laser  pulse. 

An  analyst  s^0^  of  the  radiant  heat  pulse,  from  the  fusion  ersatsd 
fireball  within  the  atmosphere  of  s  light  Ion  beam  reactor,  has  suggested 
that  it  might  lead  to  excessive  loadlrg  of  the  reactor  walls.  A  more  racant 
analysis has  found  that  seeding  this  argon  atmosphere  with  as  little  as 
0.24  of  sodium  vapour  could  degrade  this  heat  pulse  to  a  more  satisfactory 
level.  Thus  It  would  appear  that  the  alkali  vapour  could  serve  two  useful 
purposes  within  the  reactor,  'ftie  possibility  of  using  a  molecular  gas  as 
the  main  constituent  of  the  reactor  cavity  atmosphere  has  also  been  nggested/11 
In  this  paper  we  report  or  an  initial  study  of  the  Influence  of  molecular 
gas  upon  the  plasms  channel  format 1 on  potential  of  LXBCFS.  Wit rag on  was  choeeo 
as  a  representative  molecular  gas  with  no  obvious  disadvantage,  in  particular, 
w*  set  out  to  see  If  laser  resonance  saturation  of  aodlua  vapour  could  still 
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2.  USER  IQNIfcATIOK  WITH  NITROGEN  BACKGROUND  GAS 


The  presence  or  an  appreciable  background  density  of  a  molecular  gas 
such  as  nitrogen  will  Invariably  have  a  detrimental  effect  on  the  laser 
ionization  process.  Clearly  the  large  nustoer  of  Internal  energy  states 
associated  with  a  molecule  will  represent  an  effective  sink  for  the  energy 
of  both  the  free  electrons  and  the  laser  excited  gieciea.  The  most  serious 
of  these  are  considered  below: 


Botmtial  Influmos  of  jHtrogmn 

(1)  Tree  electron  cooling  through  excitation,  dissociation  and  Imita¬ 
tion  of  the  nitrogen  molecules,  vis., 

•  +Hg-»l^+2e  -  (15.58  eV)  Ionisation 

■"  ♦  Kg  **«+  ♦  W  +  2#’*  (lAtJfe  mV)  Kseoetmtiw  Ionisation 

•  ♦  1^  ♦#"  *  (9*78  eV)  Oissooiation 

The  group  of  Interactions  shown  show*  era  all  expected  to  be  wljfxir- 
tant  iue  to  the  relatively  low  electron  teqperetur*  found  In  LIBORS 
of  the  alkali  metal  vapours  (typically  £  I  *v) ,  Ih  the  case  of 

•“  *  ^  ♦  •’  *  (>  6.2  eV)  SUetnmic  Sanitation 

th«  energy  spacing  between  the  ground  state  of  ^  nd  its 

first  excited  electronic  state  (A3!*)  is  about  6.2  aV  tod  so  elec¬ 
tronic  excitation  of  should  also  ba  relatively  unitor  tent  unless 
the  alkali  seeding  drops  below  0 .14 . 

♦  Kp  -♦  *2  +  •”  “  (■•  0 .01  eV)  Rotational  Sroxtation 

has  also  been  neglected  in  the  present  snalyais  due  to  the  anil 
cross  eectlcr  end  —all  saount  of  energy  exchanged  in  this  inter¬ 
action.  Thus  the  primary  cooling  mechanism  for  the  free  electrons 
has  been  ass—md  to  result  from: 

•  +  ^g(v)  R^( u)  ♦  e’  -  (t^  -  s^)  fibrattonal  Emattation 

of  the  nitrogen  molecules,  whore  *«u-«v“»(u-w)*  o.?9  (eV) . 
^  r< /resents  the  energy  of  vibrational  state  described  by  the  w- 
rJbrttionel  qu&ntt*  madbtr  and  u  >  w. 

(ii)  Quenching  of  the  lmeer  excited  teeonance  state  through  vibrational 
excitation  of  nitrogen  molecules 

♦  «a(w)  -*le(3«)  ♦  Hj(u)  ♦  [t^  - 
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•uperelastlc  term)  end  la  given  by 


where  E?1  represents  the  iodJ.ua  resonance  energy  (equal  to  the 
laser  photon  energy). 

We  see  that  for  each  resonance  state  atom  quenched  by  collision  with 
a  nitrogen  molecule,  the  excess  energy  between  the  resonance  energy  end 
that  required  to  vlb rationally  excite  the  molecule  from  vtou,  goes  into 
translational  energy  of  the  colliding  species,  thereby  heating  the  gaseous 
mixture .  (13) 

(ill)  Depletion  of  the  alkali  seed  atme  through  the  formation  of  nitrides 
has  been  assumed  to  be  negligible  due  to  the  low  4eccwg>o*ltioc  tem¬ 
perature  of  these  easground*  (e.g.,  SaH^  decc^>osei  f or  T  >  300  K) . 

(iv)  Removal  of  low  energy  electrons  through  the  foraetioq  of  j£  appears 
to  be  negligible  due  to  the  negative  electron  affinity  of  KjA1  ^  , 

It  is  apparent  that  the  two  most  lap  art ant  interactions  to  be  considered 
are,  *L*otn?n  oojling  and  direct  ratonano#  Itvfl  quenching  through  vibrational 
excitation  of  the  nitrogen  molecules.  We  shall  now  consider  each  in  turn. 
gZsgtron  Coaling  Through  Vibrational  Excitation  of  Nitrogen 

The  total  electron- nitrogen  scattering  cross  section  la  found 
to  peak.  Just  in  the  energy  range  of  interest The  magnitude  and 
structure  of  this  cross  section  is  associated  with  a  reacnance  interaction 
that  leads  to  the  vibrational  excitation  of  Hg  through  the  momentary  forma¬ 
tion  of  Rg  according  to  Herzerherg'  and  Birtwistle  sad  Herzeifcerg. 

The  absence  of  an  electric  dipole  coosat  for  ensure*  that  the  probability 
of  direct  Vibrational  excitation  is  small 5^ 

The  LX30HS  energy  equation^  for  the  free  electrons  in  the  presence 
of  nitrogen  colt ''ilea  can  be  written  in  the  form 


It  <»,..)  -  wiiE2i  ♦  -  wi?**  '“»i  -  wV 

*\ <*21  -  *«1>VV«  *I(Eai  '  Ea.>V«E 
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e  l,  n  nc  cn  e 
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N  «,  -  II  tfH 


*  ".3ZKcnEcn  ‘  ‘  ».2<‘W>»  “A 

n>l  n^l  u>v 

■  uw 

u>w 


(1) 


where  (j)  represents  the  mean  translation  energy  and  the  density  of  the 
free  electron.  Note  ^his  equation  represents  the  sodium  resonance  state 

population  density  and  (<sm“^)  represents  the  population  density  of  the  n 
level  In  solium,  and  ny  represent  the  respective  nitrogen  population 
densities  in  the  u  and  w  vibrational  states  of  the  ground  electronic  state. 

The  first  term  on  the  RHS  of  equation  (1)  represents  the  superelastic 
heating  term  for  the  free  electrons  arising  from  electron  collision  quench¬ 
ing  of  the  laser  sustained  resonance  level  population.  E^  (j)  represents 
the  energy  separation  between  levels  m  and  n  of  the  latter  pitted  speexss, 

US*  ®co  represents  the  ionization  energy  of  level  n  of  DPS  and 
represents  the  three-body  recocfcicatloa  rate  coefficient  for  level  n. 

The  next  three  term*  of  equation  (l)  represent  the  contribution  to  the 
free  electron  energy  fraa  the  three  most  likely  electron  seed  creation 
processes.  In  the  case  of  sodium; 

(l)  Resonance  State  Two-Photon  Ionisation 
Ra(3p)  *  Pfcv  -*Ha+  ♦  •’  +  KE  (a.  1,17  eV) 


(ii)  Associative  Ionization 

!*a(3p)  ♦  Na(3p)  -*N^  +«'♦!*(•.  0.<A  eV) 

(lit)  laser-induced  Penning  Ionization 

Sa(3p)  +  R*(3p)  +  bv  -*!Ta*  ♦  Ra(3s)  +  e*  +  XS  (•»  1.17  tV) 

Ka(3p)  and  3«)  represent  a  resonance  and  ground  state  atom,  respectively, 
bv  represents  the  laser  photon  energy  and  corresponds  to  the  energy  differ¬ 
ence  between  the  resonance  and  grouni  states,  E^.  (cn*  sec)  repre¬ 

sents  the  two-photon  ionization  rate  coefficient  (sometimes  referred  to  as 
the  generalized  cross  section)  for  the  resonance  level,  and  represent, 
respectively,  the  cross  section  sad  ionization  energy  corresponding  to 
associative  ionization,  represents  the  cross  section  for  laser-induced 
Penning  ionization,  v  represents  the  mean  velocity  of  the  DPS, 

(cm  )  represents  the  single  photon  ionization  cross  section  for 
level  n,  and  the  sum  extends  over  all  levels  n  >  2  for  which  ionization  by  a 
■ingle  laser  photon  can  be  achieved,  F  (photons  cm*2  sec*1)  represents  the 
appropriate  photon  flux  per  unit  area  of  the  laser  beam,  i.e,, 

hvT  »  j  rVidv. 

C  represents  the  ion  of  energy  per  electron  due  to  the  net  colli si anally 
Induced  upward  movement  Qf  bouni  electrons  (exclusive  of  the  resonance 


c  ’  I  I  (V»  -  VU«U  *  Wa 

o*l  ton 

where  represent*  the  Ism r- induced  super® la* tic  hsstiz«  term  for 

the  free  electron*.  H#i  represents  the  rate  of  elastic  energy  transfer  to 
ion*  through  Coulcefc  Mattering  colli  si  CBS  and  if  the  electron*  ml  loo* 
h*v»  JitxveUlan  velocity  distribution  take*  the  form  (lS) 


«  5^-) 

where  T#  ( *E)  and  Ti  (*K)  represent  the  respective  electron  and  Joo  teagiera- 
ture*  (viz.,  «,  •  f  Iff.)  and  n#  (g>  ami  »t  (g)  their  respective  masMs.  *  1* 
the  electronic  charge  (esu) .  R  (csT5)  represents  the  total  neutral  atom 
density,  and  represents  the  rate  of  elastic  energy  coupling  between  the 
free  electrons  and  the  neutral  stone.  We  can  approximate  R#^  by 

/  8^,  \l/2  4l» 

H  «•  (  - -  )  (T  -  T  )ff  - - 

»*  '  *  t  ut  ea  m 


where  represents  the  electron  atom  elastic  scattering  cross  section. 

P(o)  represents  the  radiative  recombination  r%te  coefficient  for  level  n 
of  the  LFS. 

The  last  two  terms  of  equation  (l)  take  account  of  the  energy  gained 
and  lost  to  the  free  electrons  through  rib  rat  ion  si  de-excitation  and  excita¬ 
tion  respectively.  represents  the  electron  eolUslonel  excitation 

rate  coefficient  for  excitation  from  w  to  u. 

The  electron  energy  equation  can  be  divided  into  two  equations  -  one 
taking  account  of  the  rate  of  change  of  the  electron  density  and  the 
other  gives  the  rate  of  change  of  the  r»n  electron  translation  energy. 

If  we  assume  that  the  free  electrons  are  described  by  a  Maxwellian  velocity 
distribution  corresponding  to  an  electron  t«per.ture  T  ,  then  we  cm  write.  Tar 
e  2  e* 


\  k  JT  *  LIB0RS  l™  ^  *UW  (2) 

U>w  U>w 

and 

dHe 

-  LOOSS  terms  (3) 

We  see  that  the  influence  of  the  nitrogen  molecule*  i*  felt  through 
their  effect  on  T.,  and  subsequently  on  the  various  electron  colli sloosl 
rate  coefficients.  In  Order  to  determine  the  extent  of  this  Influence  we 
need  to  solve  not  only  the  atomic  population  rate  equations  and  the  energy 
equations,  but  also  the  appropriate  set  of  nitrogen  vibrational  state 
population  rate  equations, 

HT  •  \  {JT  \  } 

♦ <tTv>E  u* v*  u*  *» 

*  ViV1’1'  u+1»  wn»  *)  *  nw-i\p^u*  ®“l»  *-*»  *))]  (*•) 

The  first  two  suBnatioa  tern*  in  equation  (It)  take  account  of  the 
electron  induced  rates  of  excitation  and  de-excitation  of  the  ground  elec¬ 
tronic  state  haring  a  vibrational  quantum  tnmoer  "u".  The  sit) sequent  four 
sWstlon  terns  are  concerned  with  the  raU  of  creation  end  the  rate  of 
destruction  of  the  population  in  this  u- vibrational  state  through  vibration- 
vibration  energy  exchange  collisions  between  nitrogen  molecules .  Thus 
<enr>*  represent*  the  elastic  (gas  kinetic)  collision  rate  coefficient 
between  nitrogen  molecule*.  P(u+1,  u,  v,  v+l)  represents  the  prcbebVUty 
that  in  a  collision  between  two  nitrogen  molecules  there  will  be  a  vibra¬ 
tional  exchange  of  energy  equal  to  one  quanta  of  vibration  energy,  l.e.  the 
molecule  that  was  In  atnte-(u+l),  drops  to  u,  while  the  molecule  that  was 
in  ctgte-v  1*  raised  to  v+1, 

for  a  harmonic  oscillator  it  can  be  shown  that 

p(u+l,  u,  w,  w+l)  -  p(u,  u*l,  w*l,  w)  ^  («n)(w+l)P(l,  0,  0,  i)  (5) 
and 

P(w,  u-l»  w-1,  w)  ■  P(u-1,  u,  w,  w-i)  *  uwp(i,  o,  0,  l)  (6) 
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If  we  further  introduce  the  total  nitrogen  density, 


\mY\  (7> 

U 

end  the  average  vibrational  quan;*, 

’  ■ZV°W/nT  <*> 

w 

then  we  can  write 

Hr ' {}  \  «'»™  •  } 

vjhi  wfii 

♦  -J  [(u+l)(5+l)nu+1  -  (u  ♦  (2u+l)tf)nu  +  (9) 

where 

tV  *  [<ov>E  rTP(l4  0,  0,  l))*1  (10) 

represents  the  effective  relaxation  tine  for  vibrational  exchange, 

A  aet  of  equations  of  the  form  given  by  equation  {9)  need  to  be  aolved 
In  conjunction  with  the  normal  aet  of  UB0R3  equations  and  the  amended  elec¬ 
tron  energy  equation  -  equation  (2).  In  general  the  rate  of  electron  excitation 
from  the  lowest  vibrational  state  18  larger,  or  of  the  fame  order  of  magnitude, 
as  that  from  any  of  the  higher  vibrational  states.  In  addition  a  molecule  in 
an  excited  vibrational  state  can  be  colllslonally  de- excited  and  heat  the  free 
electrons.  Consequently,  If  we  assume  that  essentially  all  of  the  nitrogen 
molecules  are  in  their  lowest  vibrational  state,  the  rate  of  electron  cooling 
will  be  as  large  as  possible  and  Invariably  we  will  underestimate  Te-  TM*  ln 
turn  will  lead  to  a  ccmsernwtiuv  estimate  of  the  lengthening  of  the  Ionization 
time  associated  with  the  molecular  background  gas.  With  the  assumption  that  all  the 
nitrogen  molecules  are  in  the  lowest  vibrational  state  we  avoid  having  to  solve 
the  set  of  vibrational  rate  equations  and  can  write  for  the  electron  energy 
equation:  ^ 

I k  5 r  •  LIB0*S  '  %  Y  <rv>^<  (u) 

WXJ 

where  n  Is  set  equal  to  the  total  nitrogen  density. 

It  should  also  be  noted  that  <dV>^  decreases  as  V  increases  and  so  in  our 
preliminary  work  we  have  sunned  to  v  ■  U  in  equation  (ll).  These  rates  were 
numerically  evaluated  from  the  cross-sectional  data  available  ln  the  litera¬ 
ture/20^  The  subsequent  vibrational  cooling  term  in  equation  (ll)  Is  deperxient 
upon  the  electron  temperature  and  is  fitted  by  an  expression  of  the  form 
10 


In  general  we  can  write 

v  "  *21  +  "e  {*81  *  *12  +£*2a  }+£°u£<*V>£,  (13) 

n>2  u  w 

where  <5ry>^w  Is  the  nitrogen  collision  quenching  rate  Involving  vibrational 
transitions  from  u  to  v  and  A^l  represents  the  Xinsteln  spontaneous  transit loo 
probability  for  the  resonance  to  ground  transition,  if  we  as  sums  as  before  that 
the  majority  of  the  nitrogen  molecules  are  ln  the  lowest  vibrational  level  of 
ground  electronic  state  then  we  can  write 

a>? 

•here 

<«) 

w 

and 

Q 

ffow  la  the  cross  section  for  vibrational  excitation  from 

o  to  w  through  quenching  with  Ma{  3p) , 

*-{*r 

T  is  the  gas  temperature,  ar.l 

U  is  the  reduced  mass  for  the  sodi un-nlt rogen  collialon. 

(22  23) 

Ih  table  1  we  present  estimates  of  the  Ma(3p)-Hg  quenching  cro** 

section,  sad  the  corresponding  erergy  defects,  i0  the  absence 

of  nitrogen  the  electron  superelastlc  quenching  leads  to  a  minima  resonance 
state  lifetime  of  about  12.5  nsec  (compared  to  t*"  -  15 .9  nsec)  fqr  «  sodium 
atom  density  of  lO1*  cm  ^  corresponding  to  a  vapour  te^erature  of  613 
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n^AI^expCCTs)  in  our  LIBORS  code,  where  A  «  1 .3l5xlO'35,  B  »  7.3^1  and 
C  =-?.32xlO-U  and  T#  is  in  1C. 

Quenching  Through  Vibrational  Exaltation  of  Hitixxien 
In  addition  to  cooling  the  electrons  the  nitrogen  molecules  will  also 
quench  the  resonance  state  sodium  atoms  through  the  following  interaction 


«a(3p)  ♦  *2(w)  -  5a(3s)  ♦  Jlg(u)  ♦  (E?1  - 


Here  again  H?{w)  and  N^(u)  represent  nitrogen  molecule*  in  the  w  and  u 
vibrational  states,  respectively,  in  such  *  quench! rg  collision  the 
vibration  energy  of  the  nitrogen  molecule  is  increased  and  the  difference 
between  the  atom  resonance  energy  and  the  gain  in  vibrational  energy  will 
appear  ln  the  form  of  Increased  translational  energy  of  the  collision  part¬ 
ners.  This  fora  of  heating  will  be  considered  later. 

Each  tit*  a  sodiimi  resonance  ate*  la  returned  to  its  ground  level  by  a 
quenching  collision  an  additional  laser  photon  is  required  to  elevate  the 

ate*  back  to  its  resonance  state.  To  effect  this  raises  the  value  of  the 

(21l 

saturated!  speotral  irradutno*,  I^v^,  As  was  first  shown  by  the  author, 
the  resonance  and  ground  level  populations  of  the  laser  exalted  atoms  are 
only  locked  in  the  ratio  of  their  degeneracies,  g,  if  r*{  v)  »  I#(v),  where 
rV)  la  the  laser  epectral  ir reliance  (assumed  broadband  cosparcd  to  the 
atomic  linevldth).  The  *aturat*d  opactral  irradiate# 


I  (v)  ' 


8-hv3 

?(M*) 


02) 


where  hv  is  the  laser  photon  energy,  c  Is  the  velocity  of  light,  v  is  the 
frequency  of  the  laser  radlstion,  r*®  ts  the  radiative  lifetime  of  resomece 
level  and  t  1«  the  actual  lifetime  ln  the  presence  of  quenohli*  collisions. 

In  the  case  of  the  soil  via  resonance  transition. 


I#(v)  -  6.07  x  llfW  Wen"2  he'1 


ln  the  absence  of  colllsional  quenching.  If  the  laser  bandwidth  is  about 
0.1  nm  (or  *56  <lKz)  this  eorrespo-is  to  stating  that  for  saturation 


The  presence  of  2  x  10*7  cm  3  of  nitrogen  molecule,  further  reduces  this 
lifetime  to  about  7.6  nsec  which  ln  turn  increases  the  requirement  for 
laser  saturation,  Trom  equation  (12) 

I1  »  109  Vca*2 

This  poses  no  problem  as  the  values  of  laser  ir riddance  throughout  our 
analysis  considerably  exceeds  this  value. 

The  power  extracted  per  unit  volume  from  the  laser  field  restating 
from  these  nitrogen  quenching  collisions  can  be  expressed  in  the  form 

«io  •  U« 

v 

and  is  included  in  the  LlKflS  code  to  allow  for  the  additional  lasar  power 
density  absorbed  In  the  presence  of  nitrogen.  K?  1%  again  the  sodium  resonance 
state  density. 

for  times  somewhat  lees  than  the  ionization  time  wa  can  approximate 
equation  (16)  by  the  fora 

•  ",'Ai  (i7) 

whera  we  have  assumed  lasar  saturation  which  gives^1* 

*2  "  °*o  (18) 

with 

0  *  rfg  (19) 

For  *o  -  lo15  cm'3,  y?,  1017  cm'3  «  have  *  1.3  ,  l0U  ^  ^ 
tha  valuas  of  given  In  table  I.  This  represents  a  fairly  large  rata  of 
laser  energy  deposition  Into  the  nitrogen  gas.  Indeed,  the  rata  of  transit 
tl cosily  besting  the  nitrogen  molecule,  (^  U  mean  thermal  energy  of  a 
nitrogen  molecule)  i,  given  by  the  expression 

ah  (») 

At  »,  •  10U  o."3  .M  a.  •  J  i  X017  o»*5  l^uioi  (80)  tl,., 

-  1.35  *  to*  *7  eec*1  molecule**. 


I1  »  5?  Wcm*? 


and  lead  a  to  an  Increase  In  te^eratnre  of  15?  K  for  each  100  n**c  of  laser 
Irradiation  (until  appreciable  ionisation  occur*). 

This  gas  heating  would  invariably  lead  to  thermal  depopulation  of  the 
vibrational  ground  state  and  a*  indicated  earlier  would  reduce  the  electron 
cooling  rate.  This  would  noan  that  T#  would  be  higher  and  the  ionisation  tine 
reduced.  This  is  consistent  with  our  statement  that  our  calculations  would  be 
conservative.  That  Is  to  say  our  estimate  of  the  increase  in  laser  cnmrgy 
required  In  the  presence  of  will  be  larger  than  required  in  reality. 

We  have  alao  oaiitted  the  Ion- nitrogen  translation  energy  oc change  tarme 
fro*  the  ion  energy  equation.  Such  a  tern  would,  in  effect,  tend  to  prevent  the 
Ion  tenperatur*  from  rapidly  rising  towards  the  electron  te^eratu re  and 
thereby  Dean  that  the  electron-ion  energy  loss  tera  of  equation  1) 

would  remain  large  for  a  longer  period.  This  should  not  affect  the 
ionisation  tine  and  in  order  to  test  this  we  ran  one  case  where  the  Ion 
temperature  was  forced  to  stay  at  rood  temperature  (an  extreme  exxaple). 

The  results  of  this  run  confirmed  our  Justification  for  avoiding  the  com¬ 
plexity  of  Including  the  ion- nitrogen  energy  exchange  terws. 

UBOHS  Ccfyuttr  Cod*  frwau.'ts  With  •jnd  Without  nitrogen 

In  order  tc  see  which  of  the  teed  ionisation  p rooessss  is  likely  to 
dotjj-ate  at  any  giver,  value  of  last-  irradlance  we  have  evaluated  the  ionisa¬ 
tion  rate  for  each  as  a  function  of  I*  assuming  *  NQ  ■  10*^  cm”3.  The 
results  (for  the  cross  sections  used  in  most  of  our  code  runs)  are  presented 
aa  the  full  lines  in  figure  3.  It  is  evident  that  associative  ionisation  will 
dominate  for  107  Wcm'2,  while  two  photon  Ionisation  of  the  resonance 

level  will  be  the  primary  seed  lori ration  process  for  higher  values  of  laser 
lrradia.ee. 

(P5) 

Recent  experiments  y  have  suggested  that  the  laser  induced  Penning 
cross  section,  could  be  as  large  es  5-6  r  10  ^  cs  see  and  so  we  have 

added  this  result  (broken  curve)  *0  figure  3.  Clearly,  if  this  is  correct 
laser  indeed  penning  ionisatio-  is  likely  to  be  the  dominant  seed  electron 
creation  process  over  the  1^  rang’  of  Interest,  nevertheless,  most  of  our 
ItBOfiS  code  runs  used  the  smaller  value  of  2  x  10  cm  sec  for  It 

la  evident  that  If  we  had  used  the  larger  value  of  the  ionization  time 
would  have  beer  somewhat  shorter  er.i  this  would  have  the  effect  of  reducing 
the  amour.*-  of  laier  energy  reqtil-ei  to  for n  «  plasma  channel  of  a  given 
length.  have  proved  this  to  c?  •.h*  ca*e  for  cine  set  of  conditions  where 
the  i-.*ni  *  ntli  r.  burnout  time  w-.s  ir.ieel  reduced  by  about  a  factor  of  3* 


vhars  Q^—  is  taken  to  bs  tbs  paak  value  of  Q*  corresponding  to  ths 
I*  «  10®  Wcm-2  nxn,  i.e.. 


«L  -  Ov  v  J  -  w®> 

Tbs  ioni»ation  burnout  tint ,  from  squat  lorn  (21), 


(25) 


v  &  l>  + -  0  *  vn°  r 


(26) 


end  fna  equation  (22) 

V  4'  L>  “  — 


atH0.  +b(«e.  »0)(I ;  -  O1 


(27) 


The  valises  of  a,  b  am  x  were  obtained  for  the  two  sets  of  runs  shewn 
in  figure  9»  and 

(  r  1.50  x  10*4  Wcm”3  for  no  •  0 

2.0b  x  10*4  Wcm"3  for  n0  ■  2  x  1017  cm’3 
Equations  (26)  and  (27)  were  used  to  shew  the  variation  of  the  laser 
energy  with  laser  pulse  length  (chosen  to  achieve  at  least  75#  ionization  along 
a  5m  long  channel  of  1  cm2  cross  sectional  area,  see  figure  H.  The  ftoll^urve 
assumes  N  -  lO1^  cm’3  and  n  -  0,  while  the  broken  curve  assumes  N  -  10l  cm”3 

O  I  IT  ^  ^ 

and  n  •  2  x  10  c»"  .  Clearly  the  former  situation  leads  to  a  minimi*  laser 

o  j  ,5  l 

energy  density,  Er ,  •  3.2  J  »  ,  a  corresponding  laser  pulse  duration  r  - 

f  min  72  -2 

300  nsec  and  1^  -  1.07  x  10  Wcm  ,  while  the  latter  case  gives  E^n  •  9.1  J  cm  , 

»  650  naec  and  1^  -  1.**  x  107  Wcm”2. 

We  see  that  although  the  presence  of  13  tarr  of  nitrogen  increases  the 
minimum  laser  energy  density  required  for  creating  s  5m  plasma  channel  from 
3.2  Jem-2  to  9.1  Jem”2,  this  If  not  excessive.  Furthermore  the  actual  laser 
energy  required  to  create  such  a  channel  can  still  he  less  thar  U  If  a  0.1  cm 
channel  cross  sectional  area  is  acceptable  and  consistent  vi  th  laser  beam 
divergence  and  beam  shaping  optics.  It  shoild  also  be  notid  that  the  increase 
In  the  laser  pulse  duration  (300  nsec  to  650  nsec)  la  also  reasonable  as  It  would 
lead  to  an  ionization  time  that  is  within  the  time  required  to  avoid  turbulent 
o.-eak-up  of  the  channel  (>  i  uSvc). 
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Thi*  reaffirms  <xir  assertion  that  our  predictions  will  be  conservative. 

In  tabic  2  wc  present  the  results  of  a  broad  series  of  LIB0R5  code  runs 
used  to  gauge  the  likely  conditions  of  interest.  These  runs  were  also  used 
to  ascertain  the  degree  of  ionization  achievable  within  1  usee  at  a  laser 
irradlance  of  1  Wcm  .  The  results  are  presented  in  figure  U  and  from  these 
we  can  see  that  in  order  to  obtain  an  appreciable  degree  of  ionization  (i.e., 
ti  >  tO'f),  the  sodium  doping  density  J  )  ji  10^  cm  3  for  a  nitrogen  density 
nQ  ~  3  x  10l®  cm”3. 

Since  the  presence  of  nitrogen  Leals  primarily  to  a  cooling  of  the  electrons 
we  decided  to  check  the  variation  of  this  cooling  with  sodium  density  for  a  given 
laser  irradlance  and  nitrogen  density.  The  result*  are  presented  in  figure  5 
where  it  is  evident  that  both  the  plateau  and  mini  aim  values  of  the  electron 
tenperature  decrease  with  decreasing  solium  density  for  a  given  nitrogen  den- 
•  Ity.  The  corresponding  variation  in  ionization  burnout  tl^  can  be  teen  by 
reference  to  figure  f.  Clearly,  th-  Ionization  time  increases  and  the  fined 
electron  density  achieve.!  decreases  with  decreasing  tt0.  It  Is  evident  frem 
these  results  that  plasma  channel  formation  with  an  electron  density  of  about 

1015  cm”3  should  be  attainable  within  1  usee  provided  N  >  cm"3,  T*  >  10^ 

. .  -2  17-3  17  0  “  ~ 

-cm  and  1*  <  2  x  1C  cm  .  The  value  of  2  x  10  f  cm- 3  for  the  nitrogen 

density  is  in  keeping  with  current  thinking  about  the  likely  atmospheric 
conditions  within  a  light  loo  beam  fusion  reactor/*'23 

We  can  gauge  directly  Trom  the  influence  of  a  nitrogen  atmosphere  of  about 
13  to rr  upon  the  Baximua  electron  density,  peak  electron  te^eraturs  and 
ionization  burnout  time  attained  with  I*  -  10^  Wcm-2,  for  N0  -  1015  cm”3  by 
reference  to  figure  7.  The  full  lines  correspond  to  laser  ionization  in  the 

absence  of  nitrogen,  while  the  broken  curves  refer  to  the  case  where  n  - 
17-3  0 

2  x  10  cm  ,  it  can  be  seen  that  the  presence  of  the  nitrogen  reduces  both  ths 

peak  and  the  final  electron  te^i-ratur-  by  shout  a  third  and  this.  In  eTfect, 
doubles  the  ionization  time  and  prevent*  full  ionization  (81*1  as  opposed  to 
99. 9i  for  nQ  •  0)  from  being  attained. 

The  difference  in  the  power  density  extracted  from  the  later  field 
under  these  same  conditions,  can  also  be  seen  from  figure  8,  where  the  tenoral 
Variation  in  is  presented  with  that  of  both  th*  resonance  state  l»2  sad 
electron  densities,  Clearly,  although  reaches  a  peak  where  the  product 
is  a  maximim  for  both  no  -  0  ani  np  •  2  x  1017  cm'3,  inlying  that  electron 
super el ax tic  quenching  always  dominate*  the  peak  absorption.  It  Is  quits  apparent 
that  ths  presence  of  the  nitrogen  significantly  Increases  th*  laser  power  density 
’  V1  • r* !  ;•*■'!  *r  this  peak,  W’  '0;-'  ■<*^iirn  to  this  point  in  the  next  section. 

Tn  fU'.re  we  present  the  vir1*'<ou  nf  the  l*A  loniretion  » lme,  ip, 
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with  laser  irradlance  for  -  lo15  cm"3,  nQ  -  0  (P-points)  and  nQ  -  2  x  lO37 
cm  3  (0-points).  The  full  and  broken  curves  are  respective  fits  to  these  data 
by  a  relation  of  the  form: 


^7? 


(21) 


We  ctn  see  from  figure  9  that  for  values  of  laser  Irradlance  ranging 
from  10^  to  10  Wcm”2,  the  scliicz  vapour  ionisation  time  Is  roughly  doubled 
In  the  presence  of  13  torr  of  iltrogen  as  a  result  of  electron  cooling. 
nitrogen  Effmatt  on  later  ffequirtmtnta  for  Plattmg  Channel  Creation 

We  shall  now  proceed  to  estimate  the  energy  density  of  the  laser  pulse 
needed  to  create  a  plasma  channel  of  length  L  in  sodium  vapour  of  Initial 
atom  density,  HQ  -  10^  cm”3  in  the  absence  of  soy  nitrogen  background  gas. 
We  shall  then  recalculate  the  energy  density  required  In  the  presence  of 
13  torr  of  nitrogen  (I.e.,  nQ  -  2  x  1017  cm-3).  As  previously  shovn^ 
ths  ionization  time  increases  with  depth  for  as  the  laser  pulse  burns  Its 
wsy  through  the  sodium  vapour  Its  irradlance  decreases.  This 
will  aean  that  even  if  the  original  laser  pulse  were  step-like  its  tenoral 
nature  will  change  with  penetration  into  the  vapour  coluxsi .  Although  we 
are  currently  working  on  this  complex  problem  the  analysis  Included  in  this 
report  assumes  that  we  are  always  dealing  with  a  step-like  laser  pulse.  We 
therefore  assume  that  the  laser  energy  density  required  to  create  a  plasma 
channel  can  be  written  In  the  fore 


^  4  <»> 

where  (Wcm”2)  represents  the  Incident  laser  Irradlance  and  (nc) 
represents  the  7H  Ionisation  tire  after  a  penetration  depth  l  (cm).  Itoder 
condition  of  laser  saturation  It  has  been  shown^  that  the  laser  Irradlance 
at  the  end  of  the  channel  length 

A  (23) 


This  relation  assume  1  that  the  laser  power  density  Absorbed  la  Independent 
of  the  laser  Irradlance  (and  therefore  position).  This  Is  not  strictly 
true  for  I*  i  1(7  Wcm”2,  as  can  be  seen  from  figure  10.  Ths  deviation  fro* 
constant  &  It  due  prlmrlly  to  the  growing  l^jortanee  of  tlrgls  photon  Ioniza¬ 
tion  at  the  higher  values  of  laser  irrsdlanea.  Thus  to  be  conservative  v*  shall 
in  fact  writ* 
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Lb  figure  12  we  have  evaluated  the  variation  in  the  winiwes  laser 
energy  density  with  length  required  to  achieve  at  least  T5<  ionisation  in  pure 
sodium  vapour  at  about  feOn  torr  ClQ  -  1015  cm'3)  and  in  the  presence  of  about 
13  torr  (n0  «  2  *  1 O*7  cm"3)  of  nitrogen. 


v")  The  presence  of  aitnagar  not  only  cools  the  free  electrons  but  in 
quench i ng  the  laser  pumped  species  is  Itself  directly  heated. 

This  could  further  throve  the  channel  capabilities  of  UBGCtS 
through  the  creation  of  a  lav  density  care  to  the  channel. 

In  light  of  our  encouraging  theoretical  results  we  would  recoMsnd 
that  an  extensive  experimental  program  be  undertaken  to  teat  our  conclusions 
and  to  see  If  plasma  channels  created  through  laser  resonance  saturation 
offer  any  advantages  over  other  approaches  towards  the  formation  of  discharge 
guide  channels.  In  particular,  it  would  be  worth  cohering  UBQBfi  with 
direct  photolonlxatlan  (either  single  ar  two  photon)  techniques.  The  latter 
only  provides  energy  to  the  free  electrons  an  their  creation  while  in  LIB0R8 
the  laser  continues  to  provide  energy  through  superelastic  anergy  conversion. 
This  could  be  important  when  an  electron  cooling  background  gas,  such  aa 
nitrogen,  la  present.  Another  important  difference  expected  between  direct 
photoionisation  and  LXBORS  la  tha  predicted  gas  beatirg  associated  with  the 
latter.  This  leads  to  tha  fancalicu  of  a  hot,  low  density  core  to  the  channel . 

It  is  expected  that  the  choice  between  the  ca^teting  laser 
techniques  for  creating  the  discharge  guide  channel  will  probably  be  made  on 
the  basis  of  siilch  technique  can  create  a  channel  that  Is  of  high  enough 
qinllty  to  ensure  very  high  transport  efficiency  for  the  ion  beams  at  a  laser 
energy  that  is  snail  enough  to  avoid  appreciable  pellet  preheat.  We  feel 
tnat  LXB0R3  will  prove  to  be  the  nost  suitable  technique  when  all  factors 
are  taken  Into  consideration. 
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3.  CONCTUSIOftS 

An  inertial  confinement  fusion  reactor  based  upon  light  ion  beam 
drivers  will  require  an  efficient  mans  of  transporting  the  Jon  beams 
across  the  reactor  to  the  fuel  pellet.  It  appears  that  multiple  plasma 
discharge  channels  could  accomplish  this  task.  One  of  the  most  attractive 
methods  of  ensuring  that  these  discharge  channels  are  straight,  veil  defined 
and  brought  to  a  corron  focus  Is  to  create  a  system  of  veil  ionized  plasma 
channels  prior  to  striking  the  discharges.  Laser  ionisation  based  on 
resonance  saUax:tian  of  an  alkali  vapour  has  been  proposed  as  ideal  for  the 
task  of  forming  such  highly  ionized  plasma  channels  very  rapidly  with  a 
relatively  low  energy  laser. 

The  purpose  of  the  present  study  was  to  evaluate  whether  the  presence 
of  an  appreciable  background  of  colecular  nltre^ea  was  likely  to  effectively 
degrade  this  means  of  creating  the  discharge  guide  channels  or  to  increase 
the  laser  energy  requirement  to  the  point  where  it  might  become  lnpractical . 
In  an  attempt  to  compensate  for  the  many  poorly  known  cross  sections  used 
In  our  analysis  ve  have  been  conservative  wherever  possible.  Thus,  although 
the  accu.  acy  of  our  analysis  nsy  not  be  better  than  a  factor  of  two  we  feel 
that  our  conclusions  wm  be  reliable. 

Our  conclusions  are: 

(1)  The  p-lmary  effects  of  the  nitrogen  background  gas  upon  LIB CHS  are 
to  cool  the  free  electrons,  thereby  increasing  the  Ionisation  time, 
ani  to  Increase  the  laser  energy  requirement. 

(li)  Near  to  fully  (>  90$)  ionized  channels  can  be  rapidly  created  within 
a  solium  (or  lithium)  se-dei  nitrogen  atmosphere  providing  the  alkali 
Is  present  at  a  cor>cent-&tlon  of  greater  than  0.5$. 

(ill)  If  the  nitrogen  density  within  a  reactor  Is  expected  to  be  about 
2  x  1017  cm'3,  then  an  electron  density  of  close  to  lO1^  cm"3  can 
be  created  over  a  5c  path  with  a  sodium  density  of  103^  cm"3,  pro- 
viiei  the  leser  energy  denvity  Is  at  least  9.1  Jem  and  the  pulse 
duration  is  about  650  nsec, 

(iv)  The  Ionization  time  could  t"  reduced  significantly,  but  only  at  the 
cost  of  much  Ligher  laser  en*'gy  density.  For  example,  full  Ioniza¬ 
tion  could  be  achieved  within  300  nsec  along  a  channel  of  5m  length 
if  a  laser  with  an  output  anergy  density  of  23  Jar®  ware  to  be  used. 


Table  I.  8 odium/ nitrogen  resonance  quenching  cross  sections  am  energy 
defects. 
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Tabic  II.  LIE  OHS  code  results  with  and  without  nitrogen. 


*<2>  , 


1.66  *  10"1*9  cm1  tec 


Sodium 
Density 
no  (cm"3) 

Bitrogen 

Density 

no 

Cooputer 

Run 

Time  (nsec) 

Electron  Density 
Achieved 

<“  (« 3) 

Electron  Temperature 
At  end  of  Bun 

Te  (K) 

IxlO15 

0 

60 6 

9.988x10^ 

6275 

Sxio15 

0 

220 

4.982X1015 

6089 

IdO16 

0 

120 

9.95?xl015 

7581 

lxlO15 

2xl017 

1015 

8.397X1014 

3935 

5xl015 

2xl017 

412 

4 .7O2xl015 

4554 

IxlO16 

2xl017 

177 

9.5i8xlo15 

5090 

a:oli 

3x10 13 

1004 

5 .744xlOX1 

3288 

5xl315 

3xio18 

825 

3.439XI015 

3615 

lxi  a16 

3xl018 

337 

7.48Qxlo15 

4006 
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FIGURE  CAPTIONS 

1.  Variation  of  99(  Ionization  burnout  times,  with  laser  irradia^-’ 
r*  as  predicted  by  simple  model  for  lithium,  sodium, potassium  and 
nabidium,  corresponding  to  initial  density  !»Q  «  1015  of*. 

Z.  Cos*>ari8oa  of  the  temporal  variation  of  the  free  electron  line  density 
cM~z  of  *  aodium  column  evaluated  to  a  radius  of  0.6  cm  for  the  three 
cases:  (i)  Steplike,  flat  laser  pulse  for  I*  -  lo6  V  cm"2.  (li) 

Laser  pulse  Is  assumed  to  have  realistic  tenoral  variation  hut 
uniform  radial  distribution  out  to  0.6  cm.  Energy  in  the  pulse  Is 
set  to  he  0.800  J  corresponding  to  s  peak  lrradlance  of  10”  V  cm"2. 

(ill)  Laser  pulse  Is  assumed  to  hsve  both  a  realistic  tenoral 
variation  and  a  Gauaalan  radial  distribution  with  an  exponential 
radius  of  0.20  cm  and  an  energy  of  O.O91  J,  giving  a  peak  lrradlance 
of  106  W  cm”2. 

3.  Variation  of  three  primary  teed  ionisation  rates  with  laser  lrradlance 
fear  sodium  at  lo7^  cm"'*, 

<» 

3F~'  ”  t,,ocla^iv'5  ionization  rate,  |  -  laser  induced  pennlt^ 

*  L 

*! 

ionization  rate  and  |  -  two  photon  resonance  state  Ionization. 

T 

4.  Percenter  degree  of  Ionization  achieved  within  one  microsecond  for 
several  sodi’ja/nltrogen  a tro spheres  as turning  a  laser  lrradlance  of 
10°  Worn'2. 

5.  Dependence  of  electron  ar.d  ion  temperature  time  histories  on  aodium 
density  assuming  10  Wet"2  laser  li-adlance  sod  a  nitrogen  density  of 
2  x  IQ17  cm"3. 

6.  Dependence  of  electron  density  time  histories  00  sodium  density  assuming 
a  laser  lrradlance  of  lo6  Vcn"2  «nd  a  hitn*en  density  of  2  *  1017  cm"3. 

7.  A  c<mg>  arisen  of  the  t«*po-ai  variation  m  electron  density  R  ,  electron 
temperature  and  ion  ter^erature,  Tj,  with  {?,  lo17  cm"3) 'and 
without  nitrogen,  asstimirg  a  laser  lrradlance  of  106  Wcm"2  and  a 
sodium  density  of  101^  cm"3. 
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LIBORS  SIMPLE  (CONTINUITY)  MODEL 


8.  A  coe-t-iTlson  or  the  tenpO'el  variation  In  the  electron  density,^, 

resos-ivicc  state  density  and  the  absorbed  laser  density,  « 

(Wce*3),  vlth  (2  *  1017  cm*3)  and  without  nitrogen^as string  a  laser 
Irradiance  of  10^  Wcm"?  and  a  sodium  density  of  10  cm  . 

9.  Variation  of  the  Tji  ionilaMon  tt»  with  laser  irradiance  Tor  •odium 

density  1015  cm"3  assualt^  a  zero  and  ?  a  1017  c*'3  denslty^of  nitrogen. 
LDiCKS  code  points,  (•)  for  nQ  -  0  and  (o)  for  %  -  2  *  10  CM  .  ^ 

Curves  are  analytical  fits  to  the  points  of  the  for*  “  (a+totl  i  3 

10.  Dependence  or  laser  absorbed  pq«r  density  Q*  tine  history  upon  laser 
irradiance  for  .odium  density  of  10 15  c*'3  assuaing  tero  (<UL1  curves) 
and  2  *  1017  cn”3  nitrogen  (broken  curves). 

u<  Variation  of  l««r  energy  density,  K*  (Jen  2)  vith  laser  pulse  dura¬ 
tion  necessary  to  achieve  T%  ionization  at  the  end  of  a  5m  channel 
In  sodium  vapour  at  lO1^  cm”3  and  In  sodium/ nitrogen  (10  cm  /  2  x  10 
er.”")  atiiosphere.  ^  ^ 

12.  Dependence  of  the  minium  laser  energy  density,  (Jem  ),  necessary 

to  achieve  75*  ionization  burnout  at  the  end  of  channel  -  upon  the 
channel  length,  assuming  sodium/ nitrogen  atmospheres  of  (10  cm  /xero) 
and  (1015  cm”3/2  x  1017  cm”3). 
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Abstract— l  aser  energy  absorption  measurements  have  been  undertaken  in  an  experiment  involving  the 
transmission  of  a  pu'se  of  laser  radiation  through  sodium  vapour  The  wavelength  of  the  laser  was  tuned  to 
overlap  the  589  nm  resonance  transition  of  the  sodium  atoms.  Although  simple  radiation  transport  appears 
to  account  for  the  attenuation  of  the  laser  beam  at  low  values  of  incident  laser  irradiance.  anomalously 
large  absorption  has  been  observed  at  high  values  of  incident  laser  energy.  We  suggest  that  this  unom<j/ou.f 
absorption  of  laser  energy  can  be  regarded  as  evidence  of  superelastic  electron  heating  and  subsequent 
ionization. 

INTRODUCTION 

in  recent  years  several  experiments  have  provided  evidence  of  anomalous  behaviour  associated 
with  the  interaction  of  intense  laser  radiation  tuned  to  saturate  a  resonance  transition  within  a 
metal  vapour.  Measures'  was  the  first  to  suggest  the  concept  of  laser  selective  saturation  of 
atomic  transitions  for  the  purpose  of  undertaking  both  plasma  diagnostics  and  the  measure¬ 
ments  of  atomic  parameters. 

However,  several  attempts  (Sharp  and  Goldwasser,2  Burgess  and  Eckart,'  Driver  and 
Snider4)  at  observing  saturation  of  laser-induced  fluorescence  failed.  Salters  has  attributed  this 
failure  to  radiation  trapping,  while  Rodrigo  and  Measures6  and  Daily7  proposed  an  alternative 
explanation  in  terms  of  the  laser  beam  spatial  profile  and  spatial  averaging  of  the  observed 
signal.  Recently,  Salter  et  a(.H  also  failed  to  observe  saturation  of  the  laser-induced  emission  but 
they  were  able  to  confirm,  using  an  interferometric  technique,  that  the  laser  locked  the  atomic 
level  populations  in  the  ratio  of  their  respective  degeneracies,  as  predicted  by  Measures.1 

In  another  group  of  experiments  involving  laser  saturation  of  metal-vapour  resonance 
transitions.  Lucatorto  and  Mcllrath,4  Mcllrath  and  Lucatorto,10  Skinner,"  and  Young  et  aL12 
observed  substantial  ionization,  while  Leslie  et  AHegrini  et  aL ,4  and  Bearman  and 
Leventhal1*  recorded  strong  emissions  from  levels  that  were  well  above  the  laser-pumped 
transition.  These  results  can  be  understood  in  terms  of  the  superelastic  heating  and  reduced 
ionization  potential  ideas  enunciated  by  Measures.1617  The  initial  pool  of  free 
electrons  necessary  for  these  interactions  could  easily  be  created  by  the  various  collision- 
radiative  processes  discussed  by  Geltman.1*  Further  theoretical  work  on  laser  ionization  based 
on  resonance  saturation  (LIBORS)  has  been  recently  undertaken  by  Measures  et  d/.,g,2°  and 
Measures  and  Cardinal.21 

In  this  paper,  we  present  the  results  of  an  experiment  in  which  the  attenuation  of  a  laser 
pulse,  tuned  to  the  589  nm  resonance  transition  of  sodium,  was  measured  as  a  function  of  atom 
density  and  laser  energy.  We  compare  these  results  with  the  attenuation  predicted  on  the  basis 
of  simple  radiation  transport  theory. 


EXPERIMENTAL  FACILITY  AND  RESULTS 

Our  experimental  facility  includes  a  Phase-R  Model  DL-2100A  flashiamp  pumped  dye  laser, 
a  specially  designed  and  built  crossed  heat  pipe  (Heat-cross),  an  energy  meter,  two  pho¬ 
todiodes.  two  cameras  and  two  Heath  monochromators  with  accompanying  photomultipliers.  A 
schematic  of  the  experimental  arrangement  is  presented  as  Fig,  I.  The  Heat-cross  oven  used  to 
contain  the  sodium  vapour  is  shown  in  Fig.  2  and  has  four  windows.  These  provided  access  for 
the  laser  beam  and  enabled  us  to  evaluate  the  absorption  of  the  laser  pulse  after  transmission 
through  the  column  of  sodium  vapour,  the  length  of  which  was  about  40  cm  and  was  determined 
by  the  position  of  the  water  cooling  coils  and  the  argon  buffer  gas  pressure  (sec  Refs.  22  and  23). 
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Fig.  2.  Heat-cross  sodium  oven  showing  heaters  and  water  cooling  coils 


Observations  of  any  laser-induced  emission  or  laser-scattered  radiation  can  be  made  through 
the  other  pair  of  windows. 

Rhodamine  6G  in  ethanol  was  used  in  the  dye  laser  to  provide  energy  output  pulses  that 
ranged  between  100  and  400  mJ.  Tuning  of  the  output  to  the  589.0  nm  O  resonance  line  of 
sodium  was  accomplished  through  use  of  a  316  line/mm  grating  that  was  blazed  at  63°3S\  The 
spectral  width  of  the  dye  laser  output  was  further  narrowed  to  about  0.06  nm  by  the  inclusion  of 
a  one  to  two  intracavity  beam  expander.24 

The  pair  of  photodiodes  was  calibrated  against  the  energy  meter  so  that  the  shape  and 
energy  of  both  the  incident  and  transmitted  laser  pulses  could  be  monitored.  Figure  3  presents  a 
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Fig.  3.  Fig-  4- 


Fig.  3.  Photodiode  signals:  the  lower  trace  corresponds  to  an  incident  laser  pulse  of  about  300  mJ.  while  the 
upper  displays  the  inverted  signal  for  the  same  laser  pulse  after  transmission  through  the  heat-cross  in  the 

absence  of  sodium  vapour. 

Fig.  4.  Photodiode  signals:  the  lower  trace  corresponds  to  an  incident  laser  pulse  of  about  400  mJ,  while  the 
upper  trace  displays  the  inverted  singal  corresponding  to  the  distorted  laser  pulse  after  transmission 
through  a  40  cm  column  vapour  at  a  density  of  about  2.9  *  |0,5cm  \ 


representative  pair  of  photodiode  signals  corresponding  to  an  absence  of  sodium  vapour  within 
the  Heat-cross.  The  lower  trace  corresponds  to  an  incident  laser  pulse  of  about  300  mJ,  while 
the  upper  trace  displays  the  inverted  signal  that  corresponds  to  the  transmitted  laser  pulse.  The 
FWHM  can  be  seen  to  be  about  700  nsec.  The  distortion  in  the  laser  pulse  shape  arising  from 
transmission  through  40  cm  of  sodium  at  about  2.9  x  101*  cm  *  can  be  seen  by  reference  to  Fig. 
4.  When  the  sodium  vapour  pressure  is  increased  to  correspond  to  a  density  of  about 
1.5xl016em  \  almost  complete  absorption  of  the  laser  pulse  is  observed  (see  Fig.  5).  The 
variation  in  the  amount  of  laser  energy  absorbed  as  a  function  of  incident  laser  energy  for  three 
values  of  sodium  atom  density  (2.9  x  fOJ\  6.3  x  10‘\  and  1.5  x  10,hcm  *)  is  presented  as  Fig.  6. 

The  incident  dye  laser  pulse  temporal  history  can  be  reasonably  well  fitted  by  the  expression 
/'(f)  =  /«(f/r,)(  I  -  (f/n)l  exp(  -  tlr2l  where  we  have  found  that  n  ^  160.  n  =  1154,  and  = 
8 10  nsec.  Also  I{)~  I .  f  x  107  where  I1  and  /<r  are  in  Wcm  :  if  E1  is  in  joules.  In  Fig.  7  we 
present  a  representative  comparison  of  the  incident  laser  time  history  curve,  as  given  by 
the  above  analytical  expression,  and  the  experimentally  observed  points  taken  from  one  of  our 
laser  pulse  oscillograms.  The  laser  beam  has  been  assumed  to  have  a  gaussian  radial  profile  with 
an  exponential  radius  of  0.2  cm. 


THEORETIC Al.  ANALYSIS 

In  an  attempt  to  account  theoretically  for  the  observed  absorption  of  laser  energy,  we  shall 
represent  the  sodium  atom  by  the  ground  state  |1),  the  resonance  state  |2),  and  the  continuum 
| c).  We  shall  consider  a  one-dimensional  situation  and  assume  that  a  rate-equation  analysis  is 


Fig.  5.  Photodiode  signals  the  lower  trace  corresponds  to  an  incident  laser  pulse  of  about  270  mJ.  while  the 
upper  trace  displays  the  inverted  signal  corresponding  to  the  residue  of  laser  energy  detected  after 
transmission  through  a  40  cm  column  of  sodium  vapour  at  a  density  of  about  1.5  x  10’*cm  \ 
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LASER  PULSE  ENERGY  I  mJ ) 

Fig  b.  Experimentally  observed  variation  of  absorbed  laser  energy  with  incident  laser  pulse  energy  for 
three  values  of  sodium  atom  density. 


TIME  t  n  sec  } 


Fig.  7.  Comparison  of  the  laser  pulse  temporal  history  with  the  theoretical  pivfile  given  bv  Hq.  (70).  Both 
curves  are  normalized  by  (heir  respective  maximum  values. 


justified  because  of  the  broad  bandwidth  of  the  laser  and  the  fast  dephasing  time  associated 
with  the  densities  of  interest.  Under  these  circumstances,  the  appropriate  population  rate 
equation  for  the  resonance  state  can  be  written  in  the  form 

jN*:.  t)  =  N,U.  DRr. C.  0 -  N*z. t){R^:.  0+  Ai,  +  /)  +  ^N,<-\  /)(/?,  +  R, (:.  f))J.  (0 

where  /V,(:. /)  and  N:(:,t)  represent  the  ground  and  resonance  state  population  densities, 
respectively,  at  a  distance  :  from  the  boundary  of  the  sodium  vapour  and  at  a  time  t  after  the 
laser  pulse  is  sw  itched  on.  Also,  A»({  *  represents  the  radiation  trapped  spontaneous 

emission  probability;  y:\  is  the  relevant  photon  escape  factor  and  A:,  is  the  Einstein  transition 
probability:  /??,(:,  f)(  <r^F2(r,  f)|  represents  the  two  photon  ionization  rate  for  the  resonance 

state;  <rV'f  is  the  appropriate  two  photon  ionization  rate  coefficient  and  F(r, /)  corresponds  to 
the  laser  photon  flux  density  at  (:,  M,  i.e.  FU,  t)  =  /'(c,  t)ihv\  V[:*t)  represents  the  laser 
irradiance  and  lu>{  -  F?1)  the  laser  photon  energy;  RA(  -  Un)A)  represents  the  effective  collision 
rate  coefficient  for  associative  ionization  of  two  resonance  state  atoms;  R(  (:,  t )|  -  (or),  F(r.  /)| 
represents  the  laser  induced  Penning  ionization  rate  coefficient  for  two  resonant  state  atoms; 

/)|  HjJ/'tr.  t)/\{iv)dvl4ir]  represents  the  stimulated  emission  rate,  where  is  the 
Milne  stimulated  emission  coefficient,  represents  the  spectral  laser  irradiance,  and 

y?,U')  represents  the  resonance  to  ground  state  line  profile  function;  Ri;C,  Ml  =  M| 

represents  the  absorption  rate;  #  -  is  the  ratio  of  degeneracies  for  the  resonance  and 
ground  levels,  respectively. 
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If  we  invoke  continuity. 


NAz,t)  +  NAzj)  =  N» 


where  N0  represents  the  original  atom  density  (all  atoms  are  assumed  to  be  in  the  ground  state 
prior  to  laser  irradiation).  If  we  also  assume  that  the  laser  spectral  width  is  very  much 
greater  than  the  line  width  of  the  atomic  transition,  then 

^rWi.FO.O.  (3) 


(T2i-hvB2\l47rhvl.  (4) 

If  we  use  Eqs.  (2)  and  (3)  in  Eq.  (I),  we  can  write 

t)  =  ffAUi, F(:,  t)  -  N2(:,  /){( 1  +  Ft:,  t)  +  a?,  +  ogFHs.  0 

+  ^N2{:,  +  ( (tv)lF{z ,  f)]  J.  (5) 

The  corresponding  one-dimensional  radiative  transfer  equation  can  be  expressed  in  the  form 


£{F(:.  t)lc\  +  yj(:,  t)  =  {N:(r.  t)  -  f>N> <;. /)(<r',,F(;.  /)  -  2N,(:.  t)<r$FH l) 

oi  <fZ 

-~N:2(2,  t)((rv),F(:.t). 


For  the  conditions  of  our  experiment,  the  speed  of  light  may  be  regarded  as  infinite.  Under 
these  circumstances,  the  pair  of  coupled  partial  differential  equations  (5)  and  (6)  can  be 
nondimensionalized  and  expressed  in  the  form 

r)  =  <t>(x ,  t)  -  n(x%  r){  I  +  4>(x ,  r)  +  ax}>\x,  t)  +  n(.r,  r)|u  +  p4>{x ,  r))},  (7) 

oT 

r)  =  {n(.v.  r)|  I  -  bit>(x.  r)  -  n(x,  r)/>]  -  1}<^(x,  r),  (8) 

P  X 

where  r  =  A'2]t  and  x  =  GN^a1:  are  the  nondimensional  time  and  space  variables:  G  =  gl{ I  +  g) 
and  (r‘  s  (r2i(\  +  g).  The  nondimensional  excited  atom  density  and  photon  flux  density  are, 
respectively,  n  =  N2jGN0  and  4>  =  FIF The  saturated  photon  flux  density  is 

F%  =  A2\I(t‘  while  co  =  <rgFf2Mii , 

b  =  Itr’f.'FJa1.  a  ^^GN„(m  )AlA'2h  and  p  =  -GN^Ja'. 

Equations  (7)  and  (8)  do  not  lend  themselves  to  analytical  solution  and  would  have  to  be 
solved  by  numerical  techniques.  For  the  conditions  that  apply  to  our  experimental  work,  we  can 
make  a  number  of  simplifying  assumptions  that  are  reasonably  well  justified.  In  Table  1,  we 
have  listed  the  values  of  the  parameters  that  correspond  to  our  experimental  conditions.  If  we 

Table  )  The  sodium  parameters  corresponding  to  experimental  conditions. 

4ji  =  b.29xf07  y2tsec  *:  o' =  8.0*10  ’'cm2;  1.67*10  49  cm4  sec;  <<rr)4  =  2.4  x  10  *’  cm'  sec  (<n)t  = 

1.6*10  ,9cm';  G  =  0.75;  F,  =  7.9x|019  y7)  photons  cm  2 sec'1;  «  =  1.66x10  17  y2l;  h  =  3.3*10  17  y2);  Fm„ 
(corresponding  to  400 mi)  1.8*  l(F4  photons  cm  2 sec  *  1.79*  fO'/rn ;  A/(T**  =  1.5 x  I0,fc cm  \ 
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set  4>  =  0max,  N0-Noax  and  evaluate  each  of  the  terms  within  the  braces  on  the  r.h.s. 
of  Eq.  (7),  we  have 

f  1  +—[1.79  x  I05  +  5.33  x  10'7  +  n(x,  r)(2.15  x  ]0'5  +  2.02  x  10'6)]}. 

(  y2\  J 

Since  n(jc,  r)<  1.0,  it  is  clear  that,  for  the  conditions  that  prevail  in  our  experiment,  we  can 
neglect  the  last  three  terms  on  the  r.h.s.  of  Eq.  (7)  and  write 

—nix,  r)**<t>(x,  r)-  n(x,  r)[l  +  4>(x,  r)].  (9) 

dT 

In  a  like  manner,  it  can  be  easily  demonstrated  that  the  second  and  third  terms  on  the  r.h.s.  of 
Eq.  (8)  are  unimportant  for  our  experiment  and,  consequently, 

T)  =  [n(jr,T)- l]<Mjr,T).  (10) 

If  we  further  assume  that  the  laser  irradiance  changes  slowly  during  the  times  of  interest, 
then  we  can  approximately  set  (dldT)n(x,  r)  =  0,  which  leads  to  a  pseudo-steady-state  non- 
dimensional  resonance  level  density  given  by 


n(x,  r)  *= 


jK*,  t) 

1  +  <«jc.t)’ 


(11) 


Under  saturation  conditions,  t)>  1  and  n(x,  r)=  l,i.e.  N2  =  GNq .  If  Eq.(M)is  used  with  the 
nondimensionalized  radiative  transfer  equation,  we  can  write 


<£(*,  t) 

I  +  <t>(x,  r)‘ 


(12) 


which  has  the  solution 


4>(x.  t )  +  ln{<M.r,  t)}  =  4>(0,  t)  +  ln(<M0,  t)}  -  *.  (13) 

Although,  Eq.  (13)  is  a  transcendental  equation  and  has  to  be  solved  iteratively,  there  exists 
two  limiting  analytical  forms  of  solution  to  Eq.  (12).  Under  saturating  conditions,  <fi(x,  t)?>  1 
and  Eq.  (12)  yields  the  linear  attenuation  equation 

<j>(x,  t)  =  4>(Q,  t)  -  x  (14) 

first  proposed  by  Measures.16  Alternatively,  if  d>(x,  t)«1,  Eq.  (12)  yields  the  Beer-Lambert 
exponential  attentuation  equation 


<Hx,  t)  =  <£(0,  r)exp(-x).  (15) 

Using  our  earlier  definitions,  we  may  write  for  the  saturation  photon  flux  density 

r  _  8rrcAA'y2i 

F‘-W+ 1)  (16) 

and  for  the  nondimensional  scale  length  (corresponding  to  the  length  of  the  sodium  vapour 
column,  L) 

xL  -  GN0 AM2l(l  +  g)LlZircbA\  (17) 

where  A  is  the  laser  wavelength  and  AAl  corresponds  to  the  laser  line-width. 
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The  energy  in  the  transmitted  laser  pulse  is  given  by  the  relation 

Et  =  hvF,  f  <f>{xL,  t)  dr.  ( 18) 

Jo 

The  laser  energy  absorbed  in  transmission  through  the  vapour  is 

Ea  =  E‘-Et.  (19) 

COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 
In  order  to  determine  theoretically  both  the  energy  absorbed  and  the  distortion  in  the  laser 
pulse  time  history  after  transmission  through  the  40  cm  column  of  sodium  vapour  within  the 
Heat-cross,  Eq.  (13)  is  solved  iteratively  for  x  =  xL  and  the  resulting  solutions  are  then 
numerically  integrated,  as  indicated  in  Eq.  (18).  In  order  to  solve  Eq.  (13),  the  nondimensional 
incident  photon  flux  density,  0(0,  r),  has  to  be  ascertained.  This  can  be  achieved  by  employing 
the  analytical  approximation  to  the  incident  laser-pulse  time  history  described  in  the  experi¬ 
mental  section,  viz. 


0(0,  t)  = 


4.14xI05E7  r  \[ 
y2i  \\0yjl 


(20) 


where  E1  represents  the  total  energy  ( J )  in  the  incident  laser-pulse.  The  radiation  trapping 
factor,  y2u  is  treated  as  a  fitting  parameter.  In  general  the  value  of  y2\  was  selected  such  that 
the  theoretically  predicted  value  of  the  laser  energy  absorbed,  EA ,  closely  approximated  the 
observed  value. 

A  representative  example  of  our  results  is  presented  in  Fig.  8.  where  the  theoretical 
laser-pulse  shapes,  before  and  after  transmission  through  the  sodium  column  (N0  =  2.9x 
10l5cm'3),  are  compared  with  the  corresponding  experimentally  determined  laser-pulse  shapes, 
The  incident  laser  energy  for  Fig.  8  is  400  mJ  and  the  values  of  EA  are  225  mJ  (experimental 
value)  and  227  mJ  (theoretical  value  corresponding  to  y2l  =  0.97).  Although  some  aspects  of  this 
comparison  look  reasonable,  it  is  clear  that  our  theoretical  analysis  leads  to  too  small  a  degree 
of  absorption  with  regard  to  the  leading  edge  of  the  laser  pulse  and  tends  to  overestimate  the 
amount  of  absorption  in  the  tail  of  the  laser  pulse.  The  degree  of  overprediction  is  made  even 
worse  when  allowance  is  made  for  the  overprediction  of  the  laser  irradiance  available  in  the  tail 
of  the  incident  laser  pulse,  as  described  by  Eq.  (20). 

We  interpret  the  difference  between  our  theoretical  predictions  and  the  experimental  results 
(as  illustrated  in  Fig.  8)  as  evidence  of  anomalous  absorption  of  laser  energy  by  the  sodium 
vapour.  Even  stronger  evidence  of  this  anomalous  absorption  is  provided  in  Fig.  9,  where  the 


Fig.  8.  Comparison  of  the  laser  pulse  temporal  histories  (for  both  incident  and  transmitted  laser  pulses) 
with  the  corresponding  theoretical  profiles. 


544 


P.  G.  Cardinal  el  at. 


INCIDENT  LASER  PEAK  IRRADIANCE  (MWcm**) 


Fig.  9.  Comparison  of  theoretical  and  experimental  variations  of  the  absorbed  laser  energy  with  the  incident 
laser  pulse  energy  for  a  sodium  atom  density  of  about  2.9  x  to*5  cm  \  using  the  photon  escape  factor  y2i  as 

a  fitting  parameter. 


laser  energy  absorbed  by  (he  sodium  vapour  is  plotted  as  a  function  of  the  incident  laser  energy 
for  an  atom  density  of  about  6.3  x  10t5cm  \  Several  theoretical  curves  are  shown  correspond- 
ing  to  values  for  y2i  of  0.40,  0.45  and  0.55.  It  is  clear  that  a  choice  of  y2i  =  0.45  gives  a 
reasonable  fit  to  the  absorbed  energy  for  incident  laser  energy  values  up  to  about  250  mJ. 
However,  for  larger  values  of  incident  laser  energy  the  theoretical  curve  underpredicts  the 
observed  absorption. 

A  similar  result  is  found  at  other  sodium  vapour  pressures.  At  each  atom  density,  it  is 
possible  to  find  a  value  of  y2i  which  provides  a  reasonable  fit  of  the  theoretically  predicted 
energy  loss  with  incident  laser  energy  at  low  values  of  laser  irradiance  but  fails  to  account  for 
the  degree  of  absorption  at  high  values  of  incident  laser  energy.  This  can  be  seen  by  reference 
to  Fig.  10.  The  values  of  y21  that  allow  a  theoretical  fit  to  the  observed  absorption  data  appear 
to  be  somewhat  larger  than  would  be  estimated  using  the  radiation  trapping  theory  of 
Holstein.25  This  result  is  not  surprising  in  light  of  the  strong  radial  distribution  of  excited  states 
created  by  the  laser  beam. 


LASER  PULSE  ENERGY  ( m  J ) 

{PEAK  LASER  IRRADIANCE  (MWcm-*)} 

Fig.  10.  Comparison  of  theoretical  and  experimental  variations  of  absorbed  laser  energy  with  the  incident 
laser  pulse  energy  for  three  sodium  atom  densities:  2.9x10',  6.3xiO'\  and  I.5xl0,fccm  V  The  ap¬ 
propriate  value  of  the  photon  escape  factor  y:i  was  determined  by  requiring  a  best  fit  to  the  experimental 
data  at  the  low  end  of  the  incident  laser  pulse  energy  range. 
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DISCUSSION  AND  CONCLUSIONS 
We  have  been  able  to  demonstrate  that  the  energy  absorbed  from  a  laser-pulse,  tuned  to  one 
of  the  resonance  transitions  of  the  attenuating  medium,  can  be  understood  in  terms  of  simple 
radiation  transport  for  low  values  of  laser  irradiance.  However,  anomalous  absorption  is 
observed  to  arise  for  increasing  values  of  incident  laser  energy.  This  additional  absorption  of 
energy  does  not  have  its  origin  in  either  two-photon  ionization  or  laser-induced  Penning 
ionization  of  the  large  population  of  resonance  state  atoms. 

We  believe  that  this  anomalous  absorption  of  laser  energy  results  from  the  development  of 
an  appreciable  density  of  free  electrons,  which  are  then  able  to  extract  a  significant  amount  of 
energy  from  the  laser  beam  through  superelastic  collision  quenching  of  the  laser  maintained 
resonance  state  population.  This  interaction,  in  effect,  heats  the  free  electrons  and  enables  them 
to  produce  further  ionization  and  excitation.16 19~21  The  initial  pool  of  free  electrons  necessary 
for  this  mode  of  energy  extraction  is  expected  to  be  produced  by  associative  ionization, 
two-photon  resonance  ionization,  and  laser-induced  Penning  ionization.  Preliminary  measure¬ 
ments  on  the  lifetimes  of  line  and  continuum  emission  from  the  Heat  cross  substantiates  this 
possibility.  This  radiation  is  found  to  persist  for  several  microseconds  beyond  the  laser  pulse,  a 
time  consistent  with  a  decaying  dense  plasma. 

ff  appreciable  ionization  does  occur,  then  the  theoretical  underprediction  of  transmitted 
laser  irradiance  in  the  tail  of  the  laser-pulse  (see  Fig.  8)  can  be  understood  in  terms  of  a  reduced 
neutral  atom  density  along  the  path  of  the  laser  beam. 
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